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PREFACE 


This course is intended for students who have taken 
a first course in calculus and presupposes some familiarity 
with the processes of algebra and analytic geometry. It 
aims to give a concise and clear account of the most useful 
methods of solution of differential equations and to pro- ‘ 
vide ample material for practice in these methods. At the 
same time the student is kept in close touch with-ihose 
applications of the subject to which it owes much of its 
development. In these applications the problems chosen 
have been confined to those in which a student with an 
elementary knowledge of mechanics can follow the formu- 
lation of the equations. 

In the brief treatment of partial differential equations 
the topics have been selected with the object of indicating 
the possibilities of the subject and stimulating interest in 
its further study. 

I desire to express my indebtedness to my colleagues, 
Dean J. Matheson and Professor C. F. Gummer, from 
whose suggestions I have profited, and to Mr. H. S. 
Pollock, M.Sc., who drew the figures. 

N. M. 

queen’s UNIVEESrrY 
June, 1035. 
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CHAPTER I 


DEFINITIONS. FORMATION AND GEOMETRICAL 
MEANING OP DIFFERENTIAL EQUATIONS 

1. Introduction and definitions. Equations involving two 
variables a; and y have been studied in analytic geometry where 
they are seen to represent plane curves. It may be possible by 
solution of such an equation, /(»,y) = 0, to express y in terms 
oi X by means of elementary functions.f Whether or not this 
is possible, the given equation, in case f(x, y) is jsubject to certain 
very general conditions, determines y as a function of x, which 
from .the nature of the relation is called an implicit function, 
In the study of this function the locus represented by the equa- 
tion f{x, y) = 0 is of great use, as from it many properties of the 
function may be read off. 

Equations frequently arise in which the relation among the 
variables (two or more in number) involves derivatives or dif- 
ferentials of these variables. These are differential equations. 
An ordinary differential equation is one' in which only' one 
variable is independent so that the derivatives which enter are 
all total. A partial differential equation is one in which two or 
more variables are independent so that the derivatives which 
enter are all partial. 

Both ordinary and partial differential equations arise in many 
problems in the study of natural phenomena. In fact, in the 
majority of cases in which a scientific problem admits a mathe- 
matical statement, that statement is a differential equation. 
Many examples will appear in succeeding chapters. 

The elementary study of differential equations is concerned 
with solving the simpler types, that is, with finding either in 
explicit or implicit form the functions which are defined by the 
equations. A solution is then any relation among the variables 
free from derivatives which, together with the derivatives ob- 
tained from it, reduces the equation to an identity. The term 

t ^0 so-called elementary functions include rational functions, irrational 
functions involving radicals, logarithms, exponential functions, and trigono- 
mclrio functions direct and inverse. 

"‘513 



2 FORMATIOV OF DIFFERENTIAL EQUATIONS 
solution 13 also applied to the locus (curre or surface) which this 
relation represents 

The types of differential equations for which solutions can be 
found in closed form in terms of the elementary functions are 
comparatively few in niunbn^ This raises the larger question 
as to the existence of functions satisfying differential equations 
of certain forms and the properties of such functions when they 
ate not elementary The treatment of existence theorems is 
beyond the scope of this boede Our point of view la the more 
elementary one with wluch we approach the study of algebraic 
equations, viz assume in the first place that the equation has 
a solution and find a relation among the vanahles which satisfies 
it, the fact that a differential equation is satisfied of course 
proves the existence of a solution for that indiTidnal case 

The order of a di/fercntia) equation is the order of the highest 
denvative which it contains 

The de^« of a difleKUtut equation of the nth order la the 
degree to which the nth denvative is raised when the equatwa 
H rational and integral in the denvatives 

For the pnsent we shall be concerned only with ordinary 
differentia! equations 

An equation of the first mder and first degree may be written 
in either of the forms j,. 

or Jf(x,y)dr+iV(r,y)«?y =5 0 

ExiMPLZ 1 The equation ^ | or zdtf s= o, of 

the first order aod first degree is satisfied by y » 1/z, also by 
y ~clxav xy =-c for any constant value of c 

ExawtleS ~ of the first order and second 

degree has the solutions p = 2* and y = Sz also the more 
general solutions y = 2*-K and y — 3z+c where c is any con- 
stant 

Examvls 3 ^ — y = 0 ofthc second order and first degree 
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INTRODUCTION AND DEFINITIONS 
has the solutions y ■= and y = e“® and the more general solu- 
tion y = AeF^+Be-^ where A and B are arbitrary constants. 
This may be verified by direct substitution. 


2. Geometrical meaning of differential equations of the 
first order. When x and y are coimected by a relation 
f{x,y) = 0, the derivative dyjdx gives at any point {x,y) the 
slope of the curve / = 0. Thus in order that a locus satisfy 
the equation dyjdx = 1 it must be subject to the single condi- 
tion that at every point its slope is 1. It is obvious that all 
straight lines of slope 1 are solutions of this equation and that 
no other locus satisfies the equation at all its points. Again, the 

equation ^ = - states that at any point the slope of the locus 

is the slope of a line joining that point to the origin. The solu- 
tion then clearly consists of the pencil of all straight lines 
through the origin. Also a locus which at any point satisfies 

the equation ^ must cross at right angles the locus 

which belongs to the previous equation, from which it is evident 
geometrically that this^last equation represents a set of con- 
centric circles about the origin. 

Consider now the general equation of the first order and first 
degree dyjdx = /(x, y). If values be assigned arbitrarily to x and 
y a unique value will be determined by the equation for dyjdx. 
There is thus determined in the Cartesian plane a point (x, y) 
and a definite slope at that point. This combination of point 
and slope is called a line element. It may be represented on 
a figure by marking a point and drawing through it a short line 
segment with the proper slope. By drawing a sufficient number 
of these line elements in the plane we have a geometrical clue to 

the solution of the equation. Thus for the equation x^4-y = o 

dx 


we obtain the adjoined figure. This figure strongly suggests the 
existence of a set of curves such that every point on each curve 
and the tangent to the curve at this point constitute a line 
element satisfying the differential equation. The Cartesian 



4 FORUATIQN OF DIFFERENTIAL EQUATTOVS 
eq»witionofthesecTirreaiaeeentob©*y = c (ex 1 art Ijwhere 
every conataat value of c gives Ofie curve of the set 
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A coofigurstion of hoe elemenla m a plaoe u mhied by the 
faTBikar expenmeat of epnoUusg son fihoge oo a sheet of paper 
vtbieb la theo held over a magnet The iron fihogs heeome 
magnetued and airaoge themselves so that each one is a hse 
elecaeat of the owe of foroe vtuch passes through it 
I'or a differential equation of the first order aod higher degree 
there vrould correspond to each point two or moTO hue elements 
and we infer that there would pass through ever7 point two or 
more curves satisfpng the equation 
These geometncal con&ideratioDS funuah strong evidence to 
the effect that a differential eq;aation of the first order cannot 
represent onlj one curve but does in general represent a family 
of curves 

3 The derivation of first order differential equations 
Further evidence of the fact just stated is funushed on the ana* 
Ijtical side by the manner in which differential equations may 
be derived Suppose we have given an equation F(x y e) » 0 
where c is an arbitrary constaut or parameter that is for a 
definite value of e the equation represents a particular curve 
and different values assigned to e give different curves The 



DERIVATION OF FIRST-ORDER EQUATIONS 6 

curves of such a set constitute a me-paramder family. By 

differentiation we obtain ^+-:r- — 0. Between these two 

dx By ax 

equations it is in general possible to eliminate the parameter c, 
the result of this elimination being a differential equation of the 
first order. It follows that every differential equation which 
may be obtained in this way (which certainly includes all first- 
order equations which have solutions in closed form in terms 
of elementary functions) represents a one-parameter family of 
curves. 

Eor a first-order differential equation a solution which con- 
tains an arbitrary constant is called the general solution or 
primitive. A particular solution is any solution obtained from 
the general solution by giving to the arbitrary constant a parti- 
cular value. It will later appear that there exist equations 
which have solutions that are not contained in the general solu- 
tion. Such solutions are called singular. 

Example 1. The equation a^+y^—cx = 0 with c arbitrary 
represents the one-parameter family of circles which touch the 

y-axis at the origin. By differentiation we get 2z+2y^~c — 0 

ax 

and, by eliminating c between the two equations, 
x-+y-~x^x+2y^ = 0 

or x-~y^+2a^^ = 0, 

dx 

which is the differential equation of the family. 

Example 2. Bind the differential equation for each of the 
following one-ptirameter families of curves: 

(i) the family of straight lines of slope m, 

(ii) the family of concurrent lines through the point (0, 6), 

(iii) the family of equal circles of radius r with centres on the 

axis of X, 

(iv) the family of equal circles of radius r with centres on the 

line y = x. 



6 FORiUTION OF HUTERENTIAL EQUATIONS 
4. Differential equations of bigher orders. A differential 
eqnabon of tfee second order may be symbolized by 
F{x,y,dyjdx,d^ldt?) = 0 We elull take the eqoation to be 
of the first d^ree and leave to the student the vanation in the 
argument for equations of higher degrees If values be arbi- 
trarily assigned to *,y and dy/dx. the equation viU determine 
a uiuqne value for <ffy/<£c* ^us v^ue of together with 

the value assigned to dy/dx, fixes the value of j*' 

which gives the curvature of the curve described by the point 
(x,y) The conaequence is Uiat at an arbitrary point and with 
an arbitrary slope at that point there exists a curve having a 
definite curvature If the fiinction F is continuous m its four 
arguments, then a slight change Sn the position of the point 
(x.y) and m the slope dyfdx wili produce only a alight change 
m the curvature Geometncaliy it «a plausible therefore that 
the equation represents a famOy of curves of which an infinite 
BUpbn pass tbiongh every pomt All these corves will to 
general cut a given enrve, say the x ass, and so the curves 
sstufyiog the equation ace aa niunetoiu as all the duectiona 
which can be associated with aD the points on a hne To include 
all the solutions in a single pniaitire equation would necessitate 
twoconstanta and the curves constitnteatwo parameter family 
Conversely, it is clear that from a primitive equation con 
taming two arbitrary coostanta a differential equation of the 
second order may be denved Thus from /(x.y.s.h) =*= 0, in 
which a and b are constants, we obtain by two differentiations 


Sz^dydx 


0 , 


3+2i2.S , = 0, 

ftc* SxSy dx^ dy^\dz) ~ey<ir* 


and from the three equations the two constants may be eli- 
minated 

From the results obtained m the cases of first and second 
order equations the correqiondu^ facts regarding wth-order 
equations may be inferred by anahigy and wiU be merely stated 
(1) From a pmmtive equation containing n arbitrary constants 
an equivalent differential equation of the nth order may be 
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derived. (2) The general solution of a differential equation of 
the nth order involves n arbitrary constants, or the equation 
represents an n-parameter family of curves. 

Example 1. The equation (x — o)"+(y— h)* = r® with a and 
b arbitrary represents the family of all circles of radius r. To 
obtain the differential equation of this family we have 

x-a+(ij-b)^ = 0 


and 

from whiofr 1,-6 = - *-«. = Sni- 

‘ dhjlda^ d-ij[dx~ dx 

stitution of these values in the original equation gives 


1 + 


ldy\^\%ld^^ 


\dx, 


dx^ 


= r; 


wliich is seen to be a statement of the fact th^ at every point 
on every curve which it represents the curvature is constant 
and equal to 1/r. 

Example 2. Find the differential equation corresponding to 
the equation of ex. 1 in which a and r are arbitrary constants 
and b is fixed. 

Example 3. Find the differential equation of all circles in 
the plane. 


5. Boundary conditions. It has been seen that a differential 
equation has not a unique solution, that, in fact, the general 
solution contains a number of arbitrary constants equal to the 
order of the equation. When a differential equation is used in 
the solution of a problem which has a definite answer, the solu- 
tion must satisfy one or more conditions besides the differential 
equation. These are called boundary conditions. Thus a first-- 
order equation has in general a solution which satisfies the 
condition y = a when x = Xq; a second-oi-der equation has 
a solution which satisfies the two conditions y = a and 
dyjdx ~ ftj, -when x ~ x^, etc., the boundary values a, a^, etc., 
being fixed constants. 
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Exonu jClie eqnabon ^ has been seen 

to have the general solabon (x~a)*+(t/~b)* = r* Let ns im 
pose the two bonodazy coodittons p 0 when c =: 0 and 
iyjdx = 0 •when x «= 0 Then *= f* Also aince 
X— o+(y— h)dy/<ir =» 0 

It foUovra thate ^ 0 &om which h = Hence the solution 
of the equation snbject to the two boundary conditions consists 
of the two circles x*+y*±2fy = 0 Further fflostmtions of 
boundary conditions will occur in the chapters on apphcations 
of differential equations 

EZiuUFLES OW CBaTTZB I 

Find the diSmotial eqcatwiis for the foB«nn&g husihn of mm* 
t AB straight hnea m the pUoe 
S All straight lues taagrot to the circle et+y* 1 
a Ihe family of spirals r « 

4 The family of parabolas each with vertex at the ongm and axis 
along the z-axis. 

5 AUcncleainthe&retaodthirdqriadrastetoiicbingtheeoordirtato 

8 The famSy of hues each of which mekea intm^ita on the co 
ordinate axes wboeo sum is coostaot 

1 The family of lines each of which rnakee with the cooidinatoaxee 
a tnangte of constant area 

8 A family of equal parabolas snth axra along the z-axis. 

9 All ell pees of constant area with axes along the coordinate azea 

10 All circles coaxial with = 1 and » 8c 

11 All conics coofoeal with the eUtpaa ~ -f ^ ° t 

12 All circles of umt radiua pasong through the onguu 

13 All circles with centres on the line y w 2c 

14 All cILpsaa and hyperbolas with axes along the coordmate axes. 



CHAPTER II 


EQUATIONS OF THE FIRST ORDER AND 
FIRST DEGREE 

The general differential equation of the first order and first 
degree is Mdx+Ndy = 0 where M and N are functions of x 
and y. Some special cases will now be considered in which the 
general solution is readily found. The fact that these cases make 
up a comparatively small portion of all possible equations of 
this form is due to the nature of the subject. The functions 
defined by differential equations are in many cases not ele- 
mentary functions and more advanced methods must be used 
in studying them. 


6. Exact equations. If « is a function of x and y possessing 

first partial defivatives, then du = ^dx -\-^dy is called an 

ox oy 


exact differential. Thus, in order that M dx -f N dy be an exact 
differential it is necessary and suflScient that a function u exist 


such that ^dx -{-^dy ~ 3f dx -\-N dy or, in other words, such 
dx dy 

that ^ = M and ^ = N. It follows that = ~ = — — • 
dx dy dy dx dxdy 

provided that this last derivative exists and is continuous. The 

condition = ~ is therefore necessary in order that the 
dy dx 

differential be exact. That this condition is also sufiScient is 
not so obvious. The following proof makes use of a definite 
integral with a variable upper limit. Assuming that the con- 


X 

dition holds, let us form the function « = j” M{x,y) dx -f ^(y), 

where .Tq is any constant, y is indopendent of a: in the inte- 
grand, and 4>(y) is a differentiable function of y whose value 
will be assigned later. Now 
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10 EQUATIONS or FIRST ORDER ANEW FIRST DEGREE 

S” I w’’” “ / s * +■*'**’ 

daisy wiU be equal to S provided 

^ (j/) = y) oc 4{v) = ^S{x^y)dy, 

v> 

•where y# la BOtae constant Hence Mdx-\-Ndy is the exact 
diiferentiftl of 

M = J U{x s) dr + J y) dy 

If JUdz+J/dy 13 an exact differential, the equation 
Jlf ix -i-V dy a 0 13 called an exact differential equation The 
solution of such an equation may usually be wntten down by 
inspection 

Example 1 (Sx*— 2y*)<i*— <4ry+2j>)dy*A. 

This equation u found to be exact Inspection of tbe first 
term shows that the function whose differential is the left 
member of the equation must contam tbe terms :e>—ixy* and 
can hare no other term contauung x The only term in the 
equation unaccounted for is (heiefbre — 2ydjf which is the 
differential of — y* The solution of tbe equation is therefore 
ji*— 2xy*— y* = c 

Example 2 (2t*— 3x*y+4*y*)da:+(y+fte*y*— x’ldy = 0 
Examples (e*4-e>'-l-yoo8x)dr+(einx— smy+xe'ldy = 0 

7 Separation of variables A simple case of exact differen 
tial equations is ^(x)dx+^)dy = 0 which has the solution 
/ ^{xj dx + / ^(y) dy=se It follows that any process by which 
an equation may be thrown into thia form lea^ to its solution 
The process is called eepu&Uon of the variables 
Example 1 tanysec^xdx— tanxdy == 0 
Here the variables are separated by dividing by tan x tan y 
*hichg.ve> ^ 

J— 
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whence Idgtana; — logsiny = c, or tana; = 0 sin y where C = ef^. 
Example 2. ‘^ j { l — y-)dz +‘ j ( l —3 ^)dy — 0. 

Example 3. {l+x){l+y-)dx+2xydy = 0. 

Example 4. (y~+2y—3)dx +2{z'^—4:){y+2)dy = 0, 


8. Homogeneous equations. This name is given to equations 
which have or may be made to take the form dy/dx = F{y/x). 
This form is assumed by the equation M dx -{-N dy = 0 when- 
ever M and N are homogeneous of the same degree since then' 
M[N becomes a function of the single variable y/x. This Sug- 

dv dv 

gests the substitution yjx = vory = vx, from wliich ^ 

and the equation becomes v-\-x^ = F(v), or — = — j in 

dx F{v)—v X 

which the variables are separate. 

Example 1 . (x+yf dx —2x-dy = 0. 

Set y = vx, dy = vdx +xdv. This gives 

{l-\-v)-dx—2{vdx-\-xdv) = 0 

or = 0. 

dx 2dv 


Hence 

and 


0 


Z 

log* — 2tan-^p = c 


Exasiple 2. 
Example 3. 
Example 4. 


logx— 2tan-^- = c. 

z 

{z+y)dx~{x—y)dy = 0. 
(x-+y-)dx+xydy = 0. 
{x^—2y^)dx +3xy"dy = 0. 


9. Equations reducible to the homogeneous form. The 
equation ^ is homogeneous. The more general 

equation ^ i^educes to this form when the 

origin is translated to the point of intersection of the lines 
ax-}-by-i-c = 0 and OjX-f = 0. If, however, these lines 
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are parallel, that u, if a,3r+6|y = I(p+6y), then the substitu- 
tion aj+6y=s from which ttanaform the 

equation into in which the variables are 

separable 

In particular the equation 

{a*-b6y-i-c)dr-Kaj»+hjy-fe,)ify = 0 
maj always be solved by the methods of this section 
Exastple 1 (a;-f-y— l)*dr— 2(i— 2)*dy = 0 
The Imes x+y— I = 0 and x—2 = 0 intereect at (2, —1) 
Hence set ar * jr-|-2, y = r—1 This gives 
(X-{-r)*rfX-2J*dr*0 
which ta ex 1, art 8 

ESAiCPU. 2 <i—y-f l)d*— (2*— 2y--3)dy ** 0 
Set z— y a 3, dy ss dx-^ds Then 

( 2 +l)dz— { 22 — 3Miz— dr) ■« 0 

or (4— 2 )dr+( 22 — 3)<i2 = 0, which is easily integrated after 
dividing by 4—2 

ExssrPLE 3 {at+4y— 8)<fr-Kv+y— l)(/y M 0 
Exaxfls 4 (3Lr— y-f-4)dr-}-(2y— 6x+3)dy = 0 

Exssifle 5 (z-f-y-{'2>*d!r-f-(x-f-y)*i/y t= 0 

10. lotegrftbng factors. If the eqoation 
^(z,y){J/dr-{-A'rfyl= 0 

is exact, the fonctian is called an integrating factor of the 
equation Jf dx -f-iV rfy = 0, which is now supposed not exact 
Thus, for example, any fonctioD the mulhphcstion by which 
separates the variables in an equation is an integrating factor 
of that equation In other equations an integrating factor may 
Bometunes be found by inspection The equation y dr —zdyssO 
13 rendered exact by moltijdying it by 1/xy, which separates 
the variables, or by 1 /y* since =g rf^.orby l/x* since 
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INTEGRATING FACTORS 
ydx—xdtj _ _^y jjj if/fV -where / is any integrable 

X' « y- W 

function, is also an integrating factor since it tiansforms the 
equation into f{xjy)d(xji/) = 0. This iUustrates the general fact 
that if an equation has an integrating factor it has an unlimited 
number of such factors. For certain special types of equations 
rules for finding integrating factors may be given. In these 
cases, however, it will generally be found that the factors can 
be discovered without using the rules. Except in the case of the 
linear equation (art. 11) the finding of integrating factors will 
be left to the student’s ingenuity with the suggestions given in 
the following examples. 

Example 1. {c(?+y)dx—xdy = (i. 

The term a^dx is an exact differential and will remain so on 
multiplication by a function of x. Hence a function of x which 
makes ydx—xdy exact is an integrating factor of the equa- 
tion. Such a function is obviously 1/a^, which gives the equation 

xdx+ ^^^- = 0. The solutionis^-- = ’^oxx?—2y = cx. 

x® 2x2^ 

Example 2. {x^+y-+x)dx—(2a^+2y'^~y)dy = 0. 

The terms may be grouped thus, 

{a^+y-){dx—2dy)-{-xdx-\-ydy = 0 , 

from which it appears that l/{a^-{-y") is an integrating factor 
giving the solution x— 2y-f ilog(x--f y-) = c. 

Example 3. mydx +nxdy = 0. 

The variables here are separable. It is useful, however, to find 
a general integrating factor, as follows. Since 

= ia:^'^-^y^'"-^{7ny dx +nxdy), 

it follows that is an integrating factor of the given 

equation for any value of k. The fact that this constant k is at 
our disposal enables us to find an integrating factor for the 
equation 

x<^yP{mydx+iixdy)+x°‘'y^7niydx-l-v^xdy) — 0. 

For integrating factors for the first and second groups of terms 
are respectively 
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Bame integrating factor i3 for both groups of terms if 

It and ij are chosen to satisfy the equations 
tm— 1— « = fiOTi— 1— 
in— I— ^ = ti»i— t— ft 

F xastp le 4 **(y«iir+2af<fy)— y*(yiZr+4jrfyJ =r 0 

EXAiirLS 6 yiiB+a:(ay+l)rfy = 0 

Exsmfle 6 (*y— y)«ir+(2*V+*)<fy — 0 

ExiJmx 7 <7y4-3rj/*)dc— (8x+2x*y)(fy = 0 

11 Linear equations A linear differential equation of snj 
order is one which is linear in the dependent vanable and its 
detivatiTcs The linear equation of the first order may be 
WBtten^+jPy* Q wherePand Q ite functioxis of x Con 
eider first the special case in which Q u zero Py -b 0 By 
separatuig the ranables we obtain the solution logy^ f Pd!z e e 
osye/^'t** C whewf?"* ^ Differentiation of this result gv^es 
« 0 which shows that is an integrating 
factor of the equation just solved But aince it involves only 
X efi’^ is also an integrating factor ®fg + Py “ ^ Making 
use of this factor we find for the solution of the general equation 

EXAStPLS 1 !^+ycot3: = emx 

(tx 

Here JPdr = J cotrdr= logsin® 

and e*®**>v s= 6inx 

Hence ysinx = j em"xdx = i(i~Bin*cosx)+e 

t Since eny veliwol tbe>cdefii> temtegnl jFdxznakea e^^^SnintegTSt 
mg feetOT there la so need to «odb a constant of mtegration. Moreover 
i£m thjaformnla J Pit be repiseed bf J P dv it la readily »Ben that the 
■ohitios IS no mote and bo leas geneni 
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LINEAR EQUATIONS 

Example 2. = ap , 

ax X 

Example 3. (l+rc2)^+a^ = -^(1+*=). 
Example 4. ~-\-xy = x. 


12. Equations reducible to the linear form. Ifin the equa- 
tion ^ -\-Py Q, y be replaced by a function of y, the solution 
dx 

of art. 11 Tvill give the value of this function of y. Thus 

2y^-^Py- = Q has the solution = I dx+c. 

dx 

Tins is a particular case of a more general form knovm as 

Bernoulli’s equation, viz. ^-\-Py = Qy", which takes the form 

dx 


y-n^^Pyi-n _ Q or (- (1— = (l—n)Q which is 

dx ax 

linear in y^~" as dependent variable. 

In the equation M dx -\-N dy = 0 either x or y may be 

regarded as the independent variable. It is sufficient for the 

purpose of solution if the equation takes or may be made to 

take the linear form either when x or when y is regarded as 

the independent variable. 


Example 1. ^4-?!^ = x ^ y *. 
dx X 


Dividing by y* gives = x® or 

dx ~ 


a: 


dy-® 

dx 



Now 


r 3 

dx = — 31ogx = logar® and = x-®. 

•* X 


Hence 


x-^yS — j dx -fc 


= — 3x-j-c. 
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ExAMPtE 2 (*y*-y*-y‘)^+l+jr* = 0 

If 1 / be taken as independent TtmablSj ibu etjnation is 

Since J j^i = log(l+s®)* the integrating factor la (1+^)* 
and the solution is therefore 


Eeamtle 3 
Ezampu 4 


=* / (t+y*)‘y*<fy +e 
« l(l+S»*)’+e 


-z’y as y‘e' 




^a+y*}-y 


13 Change of Tariables The student mil recall from his 
study of the integral caletdos that one of the most eifeetire 
means of eralaating mdefiiute mtegrals la to make a suitable 
change of the ranable of integration -^This principle eiren 
wider application m the solution of d^erential equ^ons The 
soluiron of The homogeneous equation (art 6) was effected by 
the change of the dependent variable from y to o where y = vx, 
while the method of reducing Bernoulli s equation (art 12) to 
the linear form amounts to the substitution = r It m^ 
happen that the solution of a differential equation is {acihtated 
by a change of either one or of both vanables If the original 
variables are x and y and one of them, say z, is replaced by it, 
then u will be some functum of x and y or of x alone When 
the student has ascertained that a differential equation does not 
conform to any of the simple types his next effort should be to 
simplify its form by a change of one or of both variables For 
this no precise rules can be given dose attention to the maimer 
m which the variables occur m the equation together with a 
bttle expenence will often suggest a useful substitution 
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CHANGE OF VAKIABLES 

Example 1. (l-\-xy)dx—23?{l—2xy)dy = 0. 

K this equation be divided by it takes the form 


which suggests the substitution u — 1/a;. Then 
{u+y)du-\-2{u—2y)dy — 0 
which is homogeneous. 

Set X — rcos0, y — rsin0. Then 

a; 2 -f-y 2 _ ^2 anij -^ydy = rdr. 

Also 


xdy—ydx = rcos0(rcos0d0+8in6dr)— 

— rsin^(— rsin'^df^-fcos^rfr) = r^dO, 
The equation now becomes 

rdr jJ{a^--r-) 

WdB~ r 


or 





from which sin~^- = 0+c or r =_c^{fl+c) 


EXAMPLES ON C3HAPTEB H 

1 . 1+y* 

dx l+a;3- 

2. (i^+2y’z+a:’)tte+x’(y+3a:)rfy = 0. 

4. (2xi/+a:’)dy = 0. 

^ — *tonS =,!e’sm0. 
no 

6. x(3xyH2y5+3x)dx+y(2a:»+3x*y-3y)dy = 0. . 

7. (®^)'^+3(x4-a)i/ = 6,^ 

8. (2x-Sy4.2)(fx+(10y-4x)tfy = 0. 

9. xdy •‘rydx — (x*j/*+ai/)[ir. --tV" 

10. r8m0(?r+(r»+l)cos0d0 = 0. 

<212 _ 
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11 (ysec*J — Un*z)i* + t»Bad'y = f‘ Xa^*" 

12 — {w^^»+y*)rfy*»0 

13 i2ua«Bcc*«— — ^ 
li y*4t— (zy+l)dy «= t> 

IS (3j*y+Vl'i*-(**-*»»^<*y'=0 

18 («*»-3x)<Jy = «I* 

-17 fdy— (y-(-xUoy/2]£rs 0 

19 

19 iJog*^ + 2y = 4Vylog* 

20 *y*<lr-Ki:*+y*+y»)<fy -=0 

ai Vf**-*"^*^*^^ ■"V*') “ 

22 3y*g-**(y*+l) 

23 (»-f3y+l)iB+(2*+y-4)ify»0 

24 *^+Vlo8y ^ 

25 *'y^+2’-** 

27 (tr-yi<£*+(l+*Ky-*)rfy-o 

28 yl(9yd;s+12£tfy)4-*'(l2y^49i<Iy} • 0 

29 2(y+J)«4r-(x+y-3)»dy- 0 

90 ((>(an;r+2y(«n'x)dy— c*«c^d» K 0 



CHAPTER III 


EQUATIONS OE THE FIRST ORDER BUT NOT 
OF THE FIRST DEGREE 

If a differential equation of the first order is not of the first 
degree this circumstance in general renders the equation more 
diScult of solution. Our effort will be directed towards making 
the solution depend upon that of one or more equations of the 
first order and first degree. Certain devices wiU now be given 
which in this way effect the solution of a limited number of 
types of equations. It will be convenient to denote the first 
derivative which occurs in the equation by p. 

14, Equations solvable for p. If an equation can be solved 
for p in terms of functions of x and y the equation is replaced 
by one or more equations of the first degree. The equation 
p"+PiP"-i-fP2P’‘-2+...+P„_jp4-P„ = 0, where the P’s are 
functions of x and y, is of the first order and ntii degree. It 
may be feasible to express this equation in the factored form 
(p— 6i)(p-~Q2)—(p— Cn) — where the Q’s are functions of x 
and y. Then every solution of each of the equations p—Qi = 0 
is a solution of the given equation and every solution of the 
given equation is a solution of one of these first-degree equa- 
tions. If now it is possible to solve each of the equations 
p—Qi = 0, the solutions being P^(a:, y, c) = 0 ({ = 1, 2,..,,n), 
then the general solution of the given equation is 

Fi(x,y,c)F2[x,y,c)F3{x,y.c)...F„{x,y,c) = O.f 

Example 1, p‘—(3x^—2xy)p—G3^y = 0. 

This equation is equivalent to {p—3a?){p+2xy) = 0. 

^_3a? = 0 has the solution y—7?—c = 0 and ~-\-2xy = 0 

has the solution y~ce-^' = 0. Hence the general solution of 
the given equation is {y—s?—c){y—ce-^) = 0. 

t Two pnrticular solutions mny bo P,(x,!/,Ci) = 0 and F,{x,y,c,) = 0 ■nilh 
diftoront values of c, but the totality of all particular solutions is ovidrntly 
given by using the same nrbitraiy” constant c in each factor. 
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ExAJllFLB 2 p* 5= 9**+fi** 

This equation resolve into 

{p-arV(i+**)){i»+arV{i+a^)) * 0, 

and gives the general soiatton 

{y+c-(l+3t»)iK!/+e+(l+**)‘} * 0 

or {y+c)» « (I+3^)» 

Eiampij: 3 jif~6^+8p = 0 ^ 

Exautlb 4 p* = l—y* 

Example S p*— (»+y)p+ay *= 0 
ExAMPLS 6 xyj>*^(3C*+i/*)p+xy = 0 


15 Equations solvable for y Consider an equation which 
may be thrown uto the Ibnn 

In order to dense a method of solnng it, let a solution be 
assumed m the form y ^(z) If this be granted we get by 
.ubstKatonmd) ,S, 

whence the solution is y It follows that the solu 

tioa IS obtained by the eubsUtution of a euitsble function ^'(z) 
for p m (1} It remauiB to detemine this function ^'(z) From 
(2) It must satisfy the relation 


4' 


ex'^e^ dx' 


(3) 


and this equation, provided it can be Solved, suffices to deter 
mmo^' 

K in (3) ^'(x) be replaced by p, the equation is the result of 
differentiating (1) with regard to x, that is. 


P 


if 4-^^ 

ae'*'0pdz 


(4) 


Hence the method may be condensed into the following rule 
Differentiate the equation with regard to x Solve the resulting 
equation for p in terms of x and substitute this value of p in 
the given equation If the solotion of (4) is found m the implicit 
j'(.,r,c)=.o, (5) 



equations solvable for y 
the required solution is found by eliminating f between (1^) and 
(6). If this elimination is not easily made it is preferable to 
regard the ’two equations as parametric equations of the solu- 
tion, p acting in this case as a parameter.!* 

ExAiffiiE 1. y = aj-fp®— 2p. 

By differentiation with respect to a: we get 

whence p = 1 or 2p = x+c. p = 1 gives y = x~l which is 
a solution but, since it contains no arbitrary constant, is not 
the general solution, p = |{a;-}-c) gives the general solution, 
viz. 2 / = x+l{x+cf—x—c or (ar-fc)^ = 4(y-|-c), 

Example 2. 2y - i ~ x ^ = 2p", ^ 

By differentiation we get 

2p-t-2a: = or {p-l-*)da; = 2pdp. 

The substitution p — vx leads to the solution 
(p— x)’*(2p-f-a;) = c. 

Now p may be eliminated between this and the given equation 
or the two equations together may be regarded as parametric 
equations of the general solution. 

Example 3. 2y+x^ = xp. 

Example 4. y+zp = x*p^. 

16, Equations solvable for x. The method here is analogous 
to that used in the preceding case. If the equation may be 
made to take the form 

^=Mp), ( 1 ) 

t Blndent should distinguish carofully the three meanings oi p in this 
^iclo. In (I) p stands for dy/di; in (4) p is the dependent variable, i.e. a 
function of x which satisfies this differential equation; finally (5) and (I) con- 
B 1 uto the general solution of (1) when p is a parameter in the sense of that 
term m the study of parametric equations of curves. 
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flssame that x = ^(y) is a solution Then 

to) -4^) w 

and the Eolution 19 

Ncrw replace the functiMi 1/lf (y) p(jr) or p for abort Then 
by difierentiatmg (2) we gel * 

J (41 

1»{») ey'^’apdy' 

which 18 a differentud equation to determine p(y) If this equa 
tion can be solved forp(y), the aolotion of (1) is x ^/(yiPCy}} 
It will now be observed that (4) u the equation resulting from 
differentiatuig (1) with regard to y and wntmg dr/dy « 1/p 
Bence the rule if a diHerential equation la solved for z, dif 
ferentiate with regard to y and solve the resulting differential 
equation for p, this solution, together with the gives equation, 
form parametne equations of the general solution, from which 
the parameter p may be etimmated if this is feasible 
Exampib 1 X‘s yd-Sp— p* 

We have here a cboice of differentiating with regard td either 
f or y Differentiating with regard to y we get 


The latter equation girea y4-o Bence the general eolu 
tion results by Buhstitutingp = ±^{c-(-y)m the given equation 
Tlus gives 

* = ±2V<y+c)— c or (r+e)* = 4(y+c) 

Example 2 * = yp*+p* 

Example 3 p(r— y*) = y 


17 Clalraut’s equation Among the equations to which the 
method of art 15 apphes, the equation y = ;^-i-/(p) known as 
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Clairaut’s form, deserves special notice. Differentiation of the 


equation gives 




.dp 


or 

which is satisfied if either factor vanishes, dpjdx = 0 gives 
p — c, whence y = ca:+/(c) is the general solution. A solution 
may also be obtained by eliminating p between x-\-j'{p) = 0 
and the given equation. The nature of this solution will be 
explained in the following chapter. 

It is sometimes possible by a change of one variable or of 
both to transform an equation into Clairaut’s form. The general 
solution is then at once written down. 

Exajiple 1. y — 2px-\-yp^. 

This equation is not in Clairaut’s form as it stands. Multi- 
ply it by y. Then y- = 2ypx-{-y'p-. Now set = v. Then 

dv 1 fdv\^ , , . 

^ ~ solution is y"^ = ca:-f Jc®, or 

V- = 20x4- CP where d = Ir,. 



CHAPTER IV 


SINGULAR SOLUnOire or FIRST ORDER 
EQUATIONS 

18 Ulustrstlon. I«( us consider agam tbe equation of ex 1, 

*« *5, (1) 

The general solution of this cqnatioQ was found to be 

<x+e)* = i(yi-c). (2) 

which denotes a famOy of equal parabolas with axes parallel to 
the y axis and vertices on the tine y x The method of aolu' 
lion sure nlso (jj 

and we proceed to examine the relation of this locus to the 

equations (1) and (2) In the first place the line (8) is a solution 
of (1) since along it p » 1 and the substitution of z-'l for y 
and 1 for p reduces (1) to an identity Secondly the line (3) 
touches each of the parabolas (2), for the substitution x «« y+l 
in (2) gives (y'^l+ej* 4^+c)or (y+c— !)• 0 and the two 
values of y are equal It is sow clear why the line (3) is a solu 
tion of the differential equetion although it is not one of the 
curves given by the general eolation For every point (x,y) on 
it together with its elope pel constitute a line element be- 
longuig to one of the parabolas (2) The line is aa envelope of 
the family of curves which rnahe up the general scJution 

A solution of a differential equation which cannot be obtained 
&om the general solution by assigning a particulat value to the 
arbitrary constant is called a singular solution If the family of 
curves which coostitate the general eolation has an envelope, it 
13 clear that the envelope is a singular solution It will now be 
assumed without further mrestif^tion that all singular solutions 
are of this nature 

From the fact that the singular solution is a single curve as 
opposed to an infinite set of curves of the general solution it is 
not to be inferred that it is of little importance In many 
problems the greatest mteiest attschea to the singular solution 
and no equation which has s singular solution is considered 
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completely solved until all the solutions, singular as well as 
general, are found. 

We tirni now to a consideration of methods of determining 
the singular solution if one exists. For this purpose we may 
start either from the differential equation or from the general 
solution. The purpose of the following argument is to show 
what the facts are in the cases which are ordinarily met with. 
A complete theory of singular solutions presents difficulties 



wluch can be appreciated only by one who has a grasp of the 
salient facts.t 

19. Thep-discriminant. Returning to the example of art. 18, 
consider a point P on the upper side of the line y = a;— 1. 
Through tliis point pass two parabolas having two different 
slopes or two different values of p. As the point P approaches 
the envelope the inclinations of the two tangents through P 
become more and more nearly equal. Finally, when.P is on 
the envelope the two slopes coincide with that of the envelope. 

t Sto E. L. Ince, Ordinary Diffmntial Equations, London, 1927, p. 87, where 
dotnilcd references to the literature are given. 



28 SINQUIAB SOLUTIONS OF JTOST ORDER EQUATIONS 
The envelope is therefore distingwiabed b 7 the fact that at each 
point on it two values of p are equal which at non speciahzed 
points are different The locus will therefore be determined 
from the differential equation by expressing the condition which 
must hold between * and y m order that the equation may give 
equal values for p In the example under discussion the con 
dition IS easily found by wnting the equation as a quadratic 
in p, p*~2p+x^}/ = 0 The condition for equal roots of this 
quadratic IS \~~x-\-y « 0 which is the singular solution already 
found 

From the argument used m the case of the foregoing example 
we conclude that if the equation /(z.y.p) ss 0 haa a suigular 
solution that solution will be found by expressing the condition 
on X and y which must hold in order that this equation should 
have equal roots for p It does not follow conversely, how 
ever, that every locus which is determined m this way is 
a eingulat solutioa Other poesibihties will be comidered fit 
art 21 

Tor an equation quadratic u p, Aj^-^Bp-^C s* 0 the eon 
dition necessary and euScient for equal roots is » 0 

Generally, for an algebraic equation of any degree, a function 
qf the coefiicients is defined whose vanishing determines that 
the equation has two or more equal coots Tbia function is called 
the diacTuninant Aa applied to the equation m p, /(z, y p) « 0, 
we shall call it the p discrumnant and the equation obtained 
by equating it to zero the p-discrumnant relation From the 
theoiy of equations we leam that in order that /(z,y,p) = 0 
have two or more equal roots in p, this root must also satisfy 
the equation BflSp » O whence tbe p-discrumnant relation is 
obtained by elitninaUng p between these two equationa 

In the solution of Clammt’s equation (art 17), y = p*4-/(p) 
it appeared that a solution was obtained by eliminating p 
between the differential equation and x+f{p) = 0 The latter 
Is now seen to be the result of differentiating the given equation 

t If It should happen that the e n eriope u • bne parallel to the p axis, then 
along thu lociu dzjdy 1 /p e> 0 Tbep ducnnunant relation should give all 
the loo along which p or 1/p have equal finite value*. 
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THE p-DISCRBIINANT 
partially with regard to p. Hence according to our present 
argument the solution obtained is singular. 

Prom the above discussion it follows that one method of 
searclung for singular solutions is to find first the 'p-discriminant 
relation. This relation may resolve itself into a number of com- 
ponent equations. Each of these equations should then be 
tested to see whether it is a solution of the differential equation. 

Examples. Find the p-discriminant relation for each of the 
following equations and test whether it is a solution. 

1 , 2 / = 0 . 

2. yp"—2yp-\-x = 0. 

3. y‘P^—^yp = x'^—2y'^. 

4. p^—Zpx+Zy = 0. 

20. The c-discriminant. The singular solution of a differen- 
tial equation may also be found fi-om the general solution. The 
problem is here that of determining the envelope of a one- 
parameter family of curves. This problem is discussed in books 
on calculus, to which the student should refer to refresh his 
memory. The family of cmves F{x, y,c) = 0 with parameter c 
being given, the method consists in differentiating the equation 
with regard to c and eliminating c between the two equations 
F[x, y,c) = Q and 9F/Sc = 0. If the family of curves has an 
envelope it will be found in this W’ay, but again it does not 
follow that every locus which results from the elimination will 
be an envelope. Each locus found should be tested by substitu- 
tion in the differential equation. 

The method of finding the singular solution from the general 
solution is the same as from the differential equation. The pro- 
cess gives us the c-discriminant and the c-di^criminant relation. 
The latter may be interpreted as the condition that the equation 
F{x, 1 /, c) = 0 have equal roots in c. That this property pertains 
to the envelope is evident from the figure on p. 25. Through 
ever^' point above the line y — x — I there pass two curves given 
by two different values of c in the general solution. For a point 
on the envelope, however, the two curves through it coincide. 
It follows that the singular solution may be thought of as 
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a locus charactcnzed by equal Taluca of in the differential 
equation or by equal Taluea of c m the general solution 
Exauiu 1 Find the singular solution of 
{g—pxy = l+P* 

either from the differential equation or the general solution 
Example 2 Solve examples 3 and 4 p 27 and find the 
singular solutions &om the general solutions 

21 Tac cusp and node lod It has been pointed out that 
the arguments used in the last tiro articles do not ensure that 

|Y 



lY 

every locus found from (he p and c-discnminant relations is a 
singular solution It is interestuig to observe certain other types 
of Joel which may occur 

As a first example the equation has the 

general solution (x— = a* which denotes a set of cireles 
of radius a with centres on the axia of z The p-discnminant 
relation la y*(y*— o’) = 0 which is composed of three loci y = 0 
y a and y ^ ~a The last two satisfy the differential equa 
tion and are envelopes of the set of circles p = 0 is not a 
solution of the differential equation Through any point on this 
locus pass two cuclesof the family whose tangents at that point 
comcide But the tangent to the circles is not tangent to the 
locus This explains why the locus occurs in the p-discnnunant 
relation without beu^ a solubon of the equation Such a locus 
along which two or more curves of the family touch each other 
without touching the locus is called a tae locita Since curves 
which pass through any point on the tac locus are all distinct 


TAG, CUSP. AND NODE LOCI 29 

it is not a locus along which values of c are equal and hence is 
not given hy the c-discriminant. In the present example the 
c-disciiminant relation is f—a? == 0 which gives only the en- 



At each point on a tac locus two distinct curves of the fainilj' 
touch each other. A locus may exist at everj’ point of which 
a curve of the family has contact with itself, which hajipcns 
at a cusp on the curve. Consider the example 

The substitutions x-\-y = 7, x—y = X lead readily to the solu- 
tion 4t{x—yf = 0(x-{-y-f c)-. The p-discriininant is x—y and 
the c-discriminant {x—yf. The Une x—y — 0, however, does 
not satisfy the differential equation. It is a locus of cusps of 
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the curves of the f«m3y Hmn^h any point below the line 
y = I there pass two curves with two distinct values of p As 
the point approaches the enveh^ the curves through it ap- 
proach coincidence and the two values of p become equal This 
example makes clear that a locus of cusps if one exists will 
occur in both the p and c discnmuiant relations but will not 
ordinarily be a singular aolutioo In an exceptional case the 
curves at the cusps may be tangent to the cusp locus in which 
case this locus is also an envelope 
In the study of singular points of curves a type of singularity 
other than a cusp is a node that is a pomt where a curve 
crosses itself If each curve of a one parameter family should 
have a node will the locus of nodes be a solution of the dif- 
ferential equation t WtU its equation be contained in the p or 
t discrunioant relatioruT The answers to these questions can 
be inferred tcom a study of Uie figure on p 31 m which 
the axis of x u a node locus Through each point near the 
«-axis pass two curves which approach coincidence as the pout 
approaches the x axis Bence every point on the node locus is 
chsractented by the fact that conesponding to it » only one 
curve of the family instead of two in other words two values 
of c are equal in the general solution The equation of the node 
locus will therefore occur in the c-chscruninant relation, On the 
other hand the two values ofp at a point on the node locus are 
distinct from each other as they are at other neighbourmg 
points The p discnnuuaat lelatioa accordn^Iy will not contain 
the equation of the node locus and the latter is not a solution 
of the difierential equation Ad exception to this statement 
arises in the rare case in which the node locus is tangent at each 
node to one brunch of the curve passing through it 

22 Summary and examples Our result is that in the cases 
where such loci exist the p discnminant relation contains the 
equations of 
the envelope 
the tac locus 
the cusp locus 
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and the c-discriminant relation contains the equations of 
the envelope, 
the cusp locus, 
the node locus. 

Two further remarks may be added. A differential equation of 
the first degree in f containing only single-valued functions of 
x and \j denotes a family' of curves of which only one passes 
through each point of the plane. Such a family can have no 
envelope nor any other locus of the kinds we have considered. 



An equation of Clairaut’s form, y = px-{-f{p), has a singular 
solution which may in particular reduce to a point. Since the 
equation denotes a family of straight lines, it can have no tac, 
cusp, or node locus. 

Each of the following examples illustrates three types of loci. 

Example 1. 9p\2—yf = 

Tf this equation is solved for p the solution is found by 
integration to be (x— cp = p^(3— y). The p-discriminant rela- 
tion is {2~y)-{3—y) ~ Q and the c-discriminant relation is 
!/®(3-y) = 0. Of the three loci y = 0, y = 2, y = 3, the last 
is the only one which satisfies the differential equation. It is 
the envelope of the family of curves, y = 2, which occurs in the 
p-discriminant but not in the c-discriminant relation, is a tac 
Icwus, while y = 0, in the c-discriminant but not in the p- 
dificrimmant relation, is a node locus. 
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Exampie 2 (*— y)*(I+J>*)* =* 

Wntmg thia equation m the form 


*-y = 




and diiferentiating mth regard to x we get after a sLght 
„duct,an 


from wluch 1— p = 0 or (£r= — i-r^mSp 
{i+p*y 

the latter equation a 

x~c = a(l+P*) » 




at 


(.■-(r-oWI 

x—c 


The solution of 


These ralues substituted u the given equation give the general 
.olution (r_o,.+(,_o,l _ „l 

which denotes a set of equal 4 cas{>ed hypocyclnds with axes 
parallel to the coordinate axes and centres on the hne y « z 
The suhstitubon of p =■ 1 which was obtained in the course 
of the work gives the singular edation i— y *= 
us examine the p-discnnuoant relation DiSerentiating the 
original equation with regard top we get 

(»-yni+j^ = o*(J+P*)P* 

By multiplying the left member of the given equation by the 
right member of this and nee versa we ohtam 

(r~y)S(l+p*)»(I+p»)p* =» (*-y}*(l+p*)*(l+i>’)*P 
from which either x~y = 0ps=0p=— lorps=l p=:l 
has already been considered, p = 0 gives on substitution 
x—y s= ±0 and p = — 1 ^vea *— y = 0 There are thus five 
loci contained in the p-discnmmant relation vu 

x—y => ±iaV2 i— y =* 0 x~y = ±a 

Of these the first two are the mily ones that satisfy the dif 



33 


STOEMARY AND EXAMPLES 
lerennal equation. These lines are envelopes of the family of 
curves. It is evident from the figure that the line x — y — 0 is 
a tac locus and the lines x—y = ±o are cusp loci. 

Y 


X 


EXAMPLES OK CHAPTER IV 

For each of the follo-sving equations find the general and singular 
solutions and any tac, cusp, or node loci that exist; also sketch the 
curv’es. 

1. y = xp-Jrajp. 2. {y—px)'‘+a‘p = 0. 

3. y—px ~ V(l— p*)— pcos->p. 4. x—y = logp^p. 

6. 27(y— px)* = 4p>. 6. ity = 2pa:*— 2p*. 

7. 4p*(2p-3) = 27(y-*). 8. (*-y)»(p»+l) = (i+yp)». 

0. y-px=^j. 10. 27y-8p5 = 0. 

11. p’-Apyp+Sy- = 0. 12. 5y»p*+2aOT+a;’+4y*-4 = 0. 

13. x’yp’—(2x^y^—l)p+xy^ = 0. 

14. xyp*—(x'-j-y-~l)p+xy = 0. (Hint; Set a:* = j,3 _ 

16. (y-pa:)«(l+p>) = atp>. 16. (16p’-27)x = 24p*y, 
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APPLICATIONS OF DIFFEBENTIAL EQUATIONS 
OF THE FIRST ORDER 

It ba5 already been remarked tbat differential equations arise 
in the solution of many problems in geometry mechanics and 
other epphcationa of mathematics The present chapter is con 
cemed with some examples of problems for which the mathe 
matical formulation leads to differential equations of the first 
order In each problem the solution depends on a law which 
belongs to the snhject matter of the problem considered In 
most cases this law is one which is readily granted requiring 
little techmeal knowledge It la important to recognize how 
erer that we are not here concenied primarily with the eonect 
ness of the law the latter is merely an assumption from which 
the mathematical eolution proceeds 

33 Laws of growth I Suppose that a variable increases at 
a rate which is proportional to the vanable itself Such a law 
IS strongly sngg^ire of the loUing of a snowball It represents 
approximately the rate of growth under proper conditions of 
certain oigamc substances If the variable be called z we have 
the diffcTential equation dtjdi ss or which gives at once the 
solution X In a specific problem the constants a and c 

may visually be determined from the initial conditions It will 
be observed that the successive values of x at the ends of equal 
jnteivala of time are m geometnc progression If the constant 
0 in the differential equation is negative the equation gives the 
law of decrease or decay of the variable x Physical examples 
of this are provided in the cases of the evaporation m the open 
air of moisture from a porous eubstance the solution of salt in 
an abunihint supply of water etc 
H A biological law of growth which agrees closely with 
observation m the case of many plants and young animals is 
that by which the rate of growth is proportional jointly to the 
vansbl" itself and to its defect from a certam maximum value 
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If this maximum bo m, the law is expressed by the differential 
equation dzjdt = az[m — a:) in which the variables are separable. 

Exajeple 1. If in a culture of yeast the active ferment doubles 
in 3 hours, by what ratio vill it increase in 15 hours, assuming 
the first law of growth. 

If X(, is the amount of ferment at the beginning and x the 
amount after t hours, then dxjdt = ax and x = ce°‘. Now when 
t — 0, X = Xq, hence c = Xq. Also when t = 3, x — 2xq, hence 
2*0 = * 0 ^^“ = 2. Then when t = 15, 

X = *o6^*“ = iCo^® = 32*0. 

Example 2. If a sum of A dollars bears interest at 6 per cent, 
continuously compounded, find the amount at the end of t 
years. 

Let * be the amount at time t. If simple interest were reckoned 
on this sum for a time At this interest would be O-OBxAt. The 
compound interest continuously compounded for the time At 
differs from this by an infinitesimal of higher order than At (this 
is the law on which the solution depends; it furnishes the defini- 
tion of the term continuously compounded), Ax — 0-05*Af-f-i2 

where lim ^ = 0. From this equation, by dividing by At and 

letting At approach zero, we find dxjdt = 0-05*, whence 
a; = ce° When t — 0, x = A, hence c — A and x — Ae° 
The law which is expressed by the differential equation or by 
its solution is sometimes called the compound interest law. 

Example 3. Assiuning the temperature of the atmosphere to 
be constant for all altitudes, find an expression for the atmo- 
spheric pressure at any height h above the earth. 

Lot the pressure at height h be p and the density w. Let pg 
and Wg be the pressure and density at the surface of the earth. 

By Boyle’s law ^ ^ or w = ^ n. An increase of A7i in the 

to Wg pg-^ 

altitude produces a decrease in pressure of approximately wAh 
(approximately since this formula neglects the change in density 
throughout the increment Ah). Then Ap ~ —wAh -f where 
ij is an infinitesimal of higher order than Ah. This leads, as in 
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r ORDER EQUATIOKS 


ex 2 to the equation ^ — w s= —^p wtencep =» ee p * 

OA Pa 

VHien h <=s 0 p as Pa Hence c*=j), and p =i)oe p 


Exaufle 4 If the natural late of increase of the population 
of & coontc; IS piopoitional to the number of people bring there 
and the population incie&sea from SCO 000 to 600 000 ui ten 
years find the population as a liinctjon of the time 


Exautle fi Assume that a wet porous substance m the open 
air loses its moiatuxe at a rate piopotUonal to the moisture 
content If in a giren case half the mcutnre dries out in an 
hour when will 95 per cent of the moisture hare disappeared? 


ExAKPtx 6 If in the eqnatioa of law II m u ten tunes the 
initial value of x and x attains half ila maximum m tune tj 


show that z 


jn9'« I 

l+y* » 


34 htiscelUoeous problems Imolviog rates of change 
We shall not attempt to classify all the diverse prohlems which 
UTolre rates of change of pbyaca} quantities Problems con 
ccnung velocities pressure and temperature gradients rates of 
chemical reactions and rates of flow of liquids and gases are 
among those which lead to differential equations of the first 
order 

Exaufle 1 A tank contains I(K> gallons of hnne in which 
are £0 lb of dissolved salt f^h water is let into the tank at 
the rate of 2 gaUons a minute The solution is kept nearly 
homogeneous by stirring and is led oS through a tap at the rate 
of 1 gallon a minute How much salt does the tank contam 
after one half hour ? 

Suppose that after t minutes the tank contains z Ib of salt 
This IS dissolved m 100-t-^ gaUons In an interval of time At 
the amount of solution which is drawn off through the tap is 
At gallons If the concentrafion of the solution remained con 
stant during this interval the amoont of salt m At gallons would 
be As the concentration of the solution changes very 
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slightly during the time Af, the actual amount of salt which 

•C 

leaves the tank during At minutes differs from — ^ At by 

an infinitesimal of higher order than At, (This is the physical 


law which we assume.) Hence Ax = 


X 


100+t 


At-\-rj where 


= 0. Dividing by At and letting At approach zero we get 
At 

dt 


dx 

dt 

t = 


X , dx 
* — "r" I whence — 
lOO+t X 


lOO+t 


and X — 


100+t 


, When 


0, a; = 60, hence c = 5,000 and x = At the end of 


one half-hour t = 30 and x 


6,000 

130 


100+t 


== 38-6. 


ExAiiPLE 2. If the tank in ex. 1 contains initially 100 gallons 
of firesh water and a brine solution containing i lb. of salt per 
gallon is entering at 2 gallons a minute, the solution being led 
off through the tap at the same rate, find the amount of salt 
in the tank after one half-hour. 

Exampia: 3. Assume that when two reagents unite chemically 
the rate at which the compound is formed is at any moment pro- 
portional to the product of the amounts of the original substances 
untransformed. (This is called the law of mass action.) If a parts 
by weight of one substance combine with b parts by weight of 
another, find the amoimt of the compound formed at any time. 

If A and B are the amoimts of the two substances at the 
beginning and x the amount of compound formed at time t from 

the beginning, then ~ Solve this 

equation with the values A = B — 5, a — 1, 6 = 2 and the 
boundary conditions a: = 0 when t = 0, z = 5 when t — 10. 

Example 4. Suppose that a hot body cools in air at a rate 
proportional to the difference between the temperature of the 
body and that of the air (Newton’s law). If the air is kept at 
20° C, and the body cools from 100° to 76° in 10 minutes, when 
will its temperature become 25° ? What will be its temperature 
in half an hour ? 



38 APPLICATIONS OP FIRST ORDER EQUATIONS 
25 Determination of curves from geometrical properties 
(a) Cartavin coordinala If a property of a plane curve can be 
expressed in terms of an arbitraiy point on it and the slope at 
that point the statement of this property is a differential equa 
tion of the first order The sdotion of the latter is the ordmaiy 
equation of a set of curves having the given property If a 
singular solution exists it also baa the same property and may 
be the most significant solution of the problem It is convenient 



in stating certain properties of curves to make use of the line 
segments shown in the figure If P (z y) is an arbitrary point 
on the curve with tangent nonoal and ordinate drawn to the 
X axis the names of these segments and their values are 

subnormal =tyt8aT = y^, 
ax 


normal PN = y j, 
if = ycotv = y~. 


tangent 

Also we recall that the differential of arc on the curve is 
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DETERMINATION OF CURVES 


Example 1. Find the curve which has a constant tangent 
length. 

The differential equation is 




^ — 1 or -i-dx — 


y 


Integration of this equation gives the result 

db(a;+c)'= V(°^-y-)+alog °~ - ^^^ - ~ ^ ^ 


The curve obtained by giving c a fixed value is known as the 
tractrix. 

Example 2. Find the curve with a constant subnormal. 
Example 3. Find the curve with a constant subtangent. 
Example 4. Find the curve in which the subnormal is equal 
to the abscissa. 


(b) Polar coordinates. When the equation of a curve is given 
in polar coordinates the segments known as the polar tangent, 
normal, etc., are defined by drawing through the pole a line at 



right angles to the radius vector. Then in the figure tan ^ = 
and the values of the segments are: 



Bubtangent OT = r tan ^ = r-— , 

dr 


tangent 
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subnormal OJf *= fcot^ = 

£xAHn.s 5 Petenunetbe cnrve tvboBe polar tAngeot at any 
point IS k tunea the radius vector to that point [k > 1) 
ExauTLS 6 Deduce the polar equation of the family of 
curves in which the length of arc is proportional to the vectonal 
angle 


Also the diHereniial of arc la 




26 Trajectories The term trajectory may be applied to a 
curve which cuts a given family of curves in any stated manner 
An important case u that m which the trajectory cuts every 
member of the family at a constant angle If this angle is a 
light angle the trajectmy is called orthogonal, otherwise it la 
called oblique 

(a) Curiestun coord\MU$ 8uppoee that a given family of 
curves has the diifereniial equation /(x,p,p) <a o If [x,y,p) is 
a tnple of values satisfying this equation, then through the 
point (z,p) there passes a curve of the family with slope p If 
p’ be the slope of the orthogonal trajectory through the same 
point, thenp = ~~ andp' must therefore satisfy the equation 
= 0 In this equation p may now be replaced by 
p = dt/ldx and we get the differential equation of the orthogonal 
trajectories, y — = 0 

If a trajectory makes a constant a^igle a with the curves of 
a family given by f{z p,p) =» 0, and if p' is the slope of the 

trajectory at the point (* s) we have tana = ^ or 
, ^+PP 

p =! — - — i By the same reasoning as in the previous case 
1 -f-p tan a 
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the differential equation of the oblique trajectories is 



j)-tana: \ ^ ^ 
’ 1 +-_p tan ot/ 


(b) Polar coordinates. Orthogonal trajectories. Consider a 
family of curves whose differential equation in polar coordinates 
is /(r, B, dB/dr) = 0. At the point (r, B) on one of these curves 
let tj) be the angle between the tangent to the curve and the 
radius vector. For the orthogonal trajectory through this point 
let the corresponding angle be tfj' so that ip' = Then 

taniA' = —cot ili = — - If now {dBjdr)' be the value of the 
r dB 

derivative found from the equation of the orthogonal trajcctorjs 

we have ^ value 

\drj r dB dB \dr) 

of [dBjdr)' at the point (r, 6) must satisfy the equation 



0 . 


If wo drop the prime, which has now served its purpose, we 
have the differential equation of the orthogonal trajectories, 



Exasitle 1. Find the orthogonal trajectories of a family of 
circles which touch a given line at a given point. 

Take the given line as axis of y and the given point as origin. 
Tlie family of circles has the equation (a;— c)“+y® = ®". The dif- 
ferential equation is found by eliminating a. Thus x—a-{-yp = 0 
or a = .r-f 2/p; hence (p--t- l)y- = [x+yp)” or a?—y--{-2xyp = 0. 
The differential equation of the orthogonal trajectories is 

= 0 or [x-—y")dy—2xydx == 0 or 

y 

from which — ]-y — b or — by, which denotes a set of 

y 

circles touching the x-axis at the origin. 

Example 2. Find the orthogonal trajectories of the family 
of ellipses x--f 2y= = a~. 



♦2 APPLICATIONS OF ITOST ORDER EQUATIONS 

Example 3 Find the orthogonal trajectories of the family 
of logarithmic spirals r 

Exasitle 4 Find an eqnationfora set of ohliqne trajectonea 
to o famil} of coDcentno arcles the angle of intersection being 
40° 

27 Determination of cun ea from mechanical or physical 
properties A wide ranety of problems may be proposed in 
each of which a ctirre is required to satisfy some niecbaoical 
or physical property the mathematical statement of which is 
a diSerential equation 



Example 1 Find the shape of a plane corre such that light 
striking It from a fixed point sonrcein the same plane is reflected 
m a parallel beam Take the ongm at the pomt-Bonice and the 
z axis in the duect on of the parallel beam At {z y) ary pomt 
on the curre let the slope of the corre be p Since the inc dent 
and refiected raj’s make eqoal angles with the tangent at (z y) 

we havep = tana = or y Soar+pV 0*^ 

i+(y/^)p 9+py 

mnltipljing this eqoation by y and setting = c we get 
V =; This equation 13 m Clairaot 8 form and gives 

the soiut on y' — cz-J-Jc* whndi denotea a family of parabolas 
with axes along the x axis and common focus at the ongm 
Example 2 A dnek swuna across a nrer heading always for 
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DETERMIXATrON OP CUR'^'ES 
the point on the hank directly opposite her starting-point. 
Assuming the speed of the rivet to bo a constant v and the 
duck’s speed in still water a constant -w, find her path. 



Take the origin at the point on tlic bank towards which the 
duck swims and the axes as shown in the figure. The com- 
ponents parallel to the axes, of the duck’s velocity, me 

— ^—ucosO and ^ = v—usmB.. 
at at 


Hence 


v—ue\n9 


dx —ucosB 


. Now sin0 : 




—.y^osB-. 


X 


Hence if ^ , On substituting!/ = Kcwe 

irfl "x — = 0 from which 


^oglz+^/(l+z-)}+klogx = logo, 
or z-f-./(l+22) = ca;-fc. When a: = a, y = 0 and z = 0, hence 
c = and Z'f.y/(I-f-2-) = (a;/o}-* Solving for z we get 



It is interesting to examine this solution for different values of 
L If i < 1 the curve passes through the origin. That is, if the 
speed of the river is less than the duck’s speed in still water 
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she reachefl her destination JlJbss 1 y =» (z/a)*} wh ch 

cuts the y axis at (0 Ja) That u srhen u s= v the duck reaches 
a point on the bank distant )a from her destination If i > 1 
the y axis is asymptotic to the carve Hence when v > tt the 
duck is earned down eticam without teaching the opposite ban*. 

examples OtCEATTEB V 

1 Deternaue the curve m wh ch the teogih ef the teagent is equal 
to the mtemept made by the tangent no. the z-asis 

2 Detenmae tl e cune in vhwb the perpendicular distance from the 
ongm on any tangent is equal to the abscissa of the point of contact 

3 Cetenmne the curve whose tangent cuta oS from the coordinate 
axes mtercepta whose sum is constant 

4 Zlesennme the curve whose tangent malces with tlie coordioate 
axes a trutngle of constant area 

5 Determine the curve whose potar s^motmal la n times the polar 
fejbtangenx 

6 Find the orthogonal injectonea of the cotucs z'+oy*'— I and 
lUiutrate by a figure 

7 Fmd tbeorthogonal trajectoresofthe aet ofpanbolM S'* 

8 Find the orthogonal tnjectonee of the hyperholaa ly £* 

9 Find the orthogonal trajectories of the circUe z*i-(y 
and illustrate by a figure 

to Detemuoe the curse for which the sum of the perpeodumlatt 
from two fixed pomes to any tangent is constant. 

It Determine the cune for whicb the rectangle contamed by the 
perpendiculars from two fixed points to any tangent is constant 

12 Determine the curve whose normal at any pomt makes equal 
angles with the radius vector and the m till line 

13 Fmd the orthogonal tra)ecton*9 of the spirals r^ « o 

14 Find the orthogonal trajectones of the cacdioids r a(l— cos^l 

] 5 Find the curve for which the length of the tangent intercepted 

between the coordinate axes ts constant 

16. Show that the family of parabolas y* 4a(x+a)i9 self .orthogonal 

17 Show that a frumly of confocal ellipses and hyperbolas is self 
orthogonal 

13 Suppose that a solid sphere (d salt dissolves m running water at 
a rate proportionsl to the s<^ace area, of the sphere If hs^ the salt 
dissolves in 15 minutes in what time will it be all dissolvedt 

19 Water is poured slowly into a cask containing vinegar and the 
m xture (assumed homogeneous) is drawn oS through a tap at the same 
rate tVhal percentage of vm^ar does the cask contain when three 
tunes the volume of the cask has passed through the tap T 

20 A tank coctains 100 ga^iKisof btme with 251li ofdisso) ed salt 
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EXAJIPLES 

A solution containing J lb. of salt per gallon enters the tank at the rate 
of 2 gallons a minute and the resulting solution (assumed homogeneous) 
is drawn off at the rate of 1 gallon a minute. Find the amount of salt 
in the tank at the end of one half-hour. ■ . 

21. A 12-qunrt pail full of maple sap is exposed to a rainfall of 2 inches, 
the diluted solution overflowing. Assuming that the sugar solution is 
always homogeneous, find the ratio of the amounts of sugar in the pail 
before and after the rain. Take the diameter of the top of the pail to 
bo 1 foot. 

22. Each of two vessels contains 10 gallons of brine with 6 lb. of 
dissolved salt. Fresh water runs into the first at the rate of 1 gallon 
a minute. The resulting solution is drawn off from the first and runs 
into the second vessel at the same rate. Again the solution runs out of 
the second vessel at 1 gallon a minute. Assuming that each solution is 
kept homogeneous by stirring, find how much salt the second vessel 
contains after 10 minutes. 

23. The nir in a room is changed slowly by ventilation, warm air at 
a temperature of 40° C. being admitted. If the room air was originally 
at 10° C., what will be its temperature when a volume of warm air 
sulTioient to fill the room has been admitted t Assume the air throughout 
the room to bo of uniform temperature at any time. 

24. A mass of 10 kg. of brass (sp. heat 0-09), at a temperature of 
100° C., is placed in a vat containing 10 litres of water at 0° C. In one 
minute the temperature of the brass has fallen to 60° C. Assuming that 
no heat escapes to the walls of the containing vessel or to the surrounding 
air, find the temperature of the brass as a function of the time, 

26. Find the path of a projectile which has nn initial velocity w at an 
angle a with the horizontal. 

20. Find the shape of a curve such that rays of light striking it from 
a fi.scd point source in the plane are all reflected to a second point of 
the plane. 

27. A man swims across a river, always heading for the nearest point 
of the opposite bank. Find his path, supposing the velocity of the river 
proportional to the distance from the nearer bonk. 

28. Assuming that the temperature (in absolute units) of the atmo- 

sphere falls off imiformly from ST), at the smdaco of the earth to 0 at 
height H, at which the atmosphere disappears, show that the pressure 
at height h is given by p = where Po pressure at 

the earth’s surface and Ji is constant. 
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LINEAR EQUATIONS WITH CONSTANT 
COEFFICIENTS 

28 Tbe general linear eqoatioD Composition of solu< 
tions The general Imeax equation of the nth order is 

m which the P a and Q are functions of x Both on account of 
the extent to which their theory has been developed and of the 
frequency of their occurrence m mechanical and physical pro 
blema, linear equations constitute the meat important single 
class of ordinary differential equations No method is known 
which will solve all linear equations Solvable cases am, how 
ever, obtained by epeciaUung the coefficients P and the tight 
member Q 

The Lnear equation in which Q is replaced by sero mey be 
called the rtdu<«d equaiton of (I) The Ibllowisg principles, 
which may be vended by direct substitution, serve as a goide 
m searching for general solutions of Lnear equations 

(а) If y a/(x) IS a solution of the reduced equation, then 
y => ej{x) where c is constant is also a solution 

(б) If y >» /^(x) and y = f^{x) are two solutions of the reduced 

equation then y is a solution It follows that 

y “ IS ^ solution Moreover, if by any method 

n linearly independent 8oIutionst/|(>} /*(*) ./*(*) have been 
obtamed then y = c,/,(x)+c^j(»)+ +c„/„(z) is a solution 
containing n arbitrary constants and is therefore the general 
solution of tbe reduced equation 

(c) If y = </>{x) IS a solution of (1) containing no arbitrary 
constants and y = Ci/,(z)+c^i(x)-f- +c«/„(*) is the general 
solution of its reduced equatioo, then 

y = s5{*)+Ci/iW+ +c,A(z) (2) 

t Tha functions/, /, are sud la be Uneu-ty dependent if b set of 
Constantsa, o, not allsen^ existsstich that a,/,-i-a,/,-f -f-a,/, =± 0 
If no swh set of rensltints exists the fanctioiis are linearly indcpcnd nt 
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is the general solution of ( 1 ). For the substitution of ^(a;) in 
the left member of ( 1 ) reduces it to Q and the substitution of 
Cj/j(a:)+...+c„/„(x) reduces it to zero. Hence (2) is a solution 
of ( 1 ), and since it contains n arbitrary constants it is the general 
solution. In this solution of (1) the function <f>{x) is called a 
•parlicular integral and Ci/i(ar)+...+c„/„(a;) the complementary 
function. 


29. Linear equations with constant coeflScients and right 
member zero. The general equation of this typo is 

( 1 ) 




-+an-i^+a„2/ = 0. 


The substitution of y = e'“ in the left member of this equa- 
tion gives Conse- 

quently y = c”*® is a solution of ( 1 ) provided m is a root of the 
algebraic equation = 0 , which is 

called the auxiliary equation. If the n roots of this equation are 
TTii, m„ and are all distinct, then the general solution of 
the differential equation is 

ExAjrPLE 1, 6y = 0. 'i-> 

dx^ dx ^ 


The auxiliary equation is m-— m — 6 = 0 which has the roots 
3 and — 2 . The general solution is therefore y = 026 “^. 

Example 2. l3^-f-12y = 0. 

dx? dx 

Examples, 2^— 4^4-w = 0. 

dxr dx 


30. The case of equal roots of the auxiliary equation. If 
two roots mi and m^ of the auxiliary equation are equal, then 
the solution of the differential equation corresponding to these 
equal roots is = (cj-l-Cjle’"*® which contains 

essentially only one arbitrary constant Ci-f-Co. It follows that 
in this case the solution obtained in the preceding paragraph 
contains fewer than n arbitrary constants and is therefore not 
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the general solnbon The result obtained in the following 
example indicates the nature of the solahon m such a case 
A second method is giesa m art 33 

Example So^-{-o*sf =* 0 


This equation has the anxQiaiy equation {m— o)* and 

hence y = cs" u one solution Let ns replace the constant c bj 
a function r(z) and inquire whether v can be determined other 
than a constant so that y » satisfies the dilTeiential eqna 
tion t If y =<= i«“ we have 




equation u satisfied 

iPe/d^ *a Q or i£ n a» ti+t,* TVas give* y * ** 

the general solution of the given equation 


31 The case of Imaginary roots of the auxUlaiT equation 
It remains to consider the form of solution corresponding to 
imaginary toots of the anxihary equation As this equation has 
real coefficients imaginary roots if such exist most occur in 
pairs of the form a± where * = V— 1 In this caae the expres 
sons and e<“ formaUy satisfy the differential equa 

tion but we are faced with the problem of giving a meaning 
to these exponential functions with imaginaiy exponents We 
shall avoid this problem by solving directly the second-order 
equation whose auxiliary equation has the roots a±i^ The 
extension to equations of higher order will be made m art 32 

Consider first the differential equation 

which the auxiliary equation = 0 has the imaginary 

roots This equation beecooes an identity on the substitu 

t Thu u • simple ezsmple of the method of T&nst on of parsmeten Cf 
art 51 
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tion of cos^x or of siapx for y. Hence the general solution is 
y = Cicos^x+Cgsin^. Take next the equation 


g_2.|+(„.+«, = 0 


( 1 ) 


for ■which the atoiliary equation is m^—2am-\-o?+pr = 0 with 
roots a±tj9. This equation is transformed to the former case 
by the substitution of c“*z for y. For then 


y=e' 


whence 


cPz 


and (1) transforms into = 0. Since this equation has 

the general solution z — A cos /9a:+ B sin Px, it follows that ( 1 ) has 
the general solution y = e^{A cosjSx+BsinjSx). An equivalent 
form for the solution is y = A:e'“cos(j9x+y) "''dtii the arbitrary 
constants k and y. It is verified by expansion of cos(j3x+y) 
that the constants of the two solutions are connected by the 
relations A = fccosy, B — — Asiny. 


Exasiple 1. 


^?'_2^4-2y 


Tlie auxiliary equation is m"— 277i+2 = 0 which has the roots 
1±»- Hence the solution may be •written in either of the forms 
y = e*(A cosx+H sinx) or y = ^•e'co5{x+y). 

Example 2. ^-f-3^+4y = 0. 
dx~ dx 


Example 3. 


g-I0|+6Iy=0. 


32. The operators D and /(E). In the further study of linear 
equations great advantage is found to result from using the 
symbols JD, D-, JC^, etc., to stand for the differential operators 
d/dr, d-fds?, dPfdx^, etc. The linear equation of the nth order 
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mth comst&nt coefficients and ngtit member a (anctioa of x may 
then be written 

>+ -Hi, x2>+a,)y« « 

Theoperatori)*-fax2)" •+ +a, has the /om of a polynomial 
and may be denoted bj/[Z)) Tbe differential equation is then 

The nsefnln^ of the symbohc operaton D and /(D) r^nlts 
&om tbe cucnn^tanee that la many operatiom D bebavea as 
if it were an algebraic qnantity Thu u due to the fact that the 
symbol obeys the following laws of algebra 
(I) it IS dutnbutiTO over the terms of a stun 
D(tt-|-c) a= Dtt'i'Dv 

(3) it u commutatiTe with constants D(ett) » cPu 
(3) it obeys the index law « D^**u- 

It 18 to be noted that tbe {ooperty which differentiates it from 
an algebraic quantity u that it u not commntatiTe with van 
ab\es that is £>uc u not equal to uDc It is for this reason 
that tbe chief nsefiilness of the symbol arises in the stnt^ of 
equations with constant coeffictents 
A fint result of these properties u that if the polynomial 
j{D) be separated into factora (D— a,)(D— a,) (D— a,) which 
may be wniten m any order the two symbols regarded as 
operators hare tbe same meamsg For the change from 
(D— ai)(D— a,) (D— to /(D)y requires the use of only the 
three laws enumerated above 

Suppose now that f{D) is resolved into two factors so that 
/(D)y — /i(D)^(D)y If y u such a function of x that^(D)y s 0 
it follows that /(D)y = 0 Thus a solution of /(D}y ^ 0 u 
obtained from any factor of /(D) If such a factor is 
D*— 2aD4-«^+^ 

the corresponding solutimi of /(D)y = 0 u 
e^ ( A 008^+ fi ran 

The result of art 31 fcv second-order equations u therefore 
extended to an equation of any order in which the ansihaiy 
equation has a pair of imaginary roots 
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33. The device of factoring the operator f{D) furnishes an 
alternative method of treating the equation J{D)y = 0 in case 
the auxiliary equation has repeated roots. The method will be 


illustrated by solving the equation 




da? 


2a — -‘rOry = 0 which 
.dx 


was solved by a different method in art. 30. In symbolic form 
the equation is {D—a)~y = 0 or (D—a){D—a)y = 0. Denote 
{D~a)y by v. The equation is then (J?—a)v = 0 firomwWch 
V = ce°^. Hence y satisfies the equation (Z> — a)y = 
equation of the first order which gives y = e“(CiX+C2). By a 
repetition of this process the student should show that\ the equa- 
tion {D—ayy — 0 has the solution y = e'^{Cia?-yc2X-\-c^) and 
that {D—aYy = 0 has the solution 


Example 1. {D—b){D—a)-y—0. vX 

Setting {D—a)y = v vfe have (D—b)(I)—a)v 
V = CjS®*. Then {D—a)y = Cie'’^-fC2e‘“ gives 

y = e'“|J e'’*-t-C2 c“*) dx -J-Caj = C7ie*'*-}-(c2X-l-C3)e 


0 whence 

u 


where Cj — ejib—a]. 


Example 2. {D—b)-[D~ayy = 0. 

Examples. XD-~{-iyy 

The auxiliary equation has the roots each repeated. By 
art. 31 y = cos x and y = sin x are solutions. The above results 
suggest the general solution y = (c2-j-C2x)cosx-j-{c3-}-C4x)8inx, 
which may bo verified. 


34. The inverse operator It aheady been pointed 

out that the general solution of the equation f{D)y = Q con- 
sists of two parts, a particular integral which contains no 
arbitrary constants and the complementary function which is 
the general solution of the reduced equation. Methods of finding 
the complementary function having now been considered, it 
remains to devise methods of finding a particular integral. A 

particular integral may be symbolized by y = ™ Q. The 

J i 
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gymbol IS thea a new operator defined by the fact that 
Q IS such A fimctiOQ P that 0 = /(D) P This suggests the 
analogy of other mrerso operations Thns if © = sin P then 
P = em^Q ifQ = e** then P — log and if G =* dP/dr 
then P = J G dx This last lUustxation is indeed a particulai 
case of the operations nnder discussion, for if G DP, then 
P = ^G^&^theeymbol^orD ' denotes integration Further 


more, rrbile the result of the operation /(D) P is to giTe a defimte 
function G> the result of the operation ^ indefinite in 

much the same way as ordinary integration la indefimte We 
have by definition other hand, 

^^/(D)G 13 such a function P that /(D)P -=/(D)G<*ai3 two 
different functions P and Q may satisfy this relation Thus, 
for example, D’** » D*(**+Ci**+C|ar+Cj) where the e’a are 
arbitrary constants When therefore we wnte “ C 

we are to understand that the right member of this equation 
IS one of the many functions denoted by the left member 


35. Tbeparticularlntegralof/(D)ysz G Generalmethods. 

First method The Unear equation of the first order is 
(D— a)y » G fts solution has been found to be 
y = e"J<-"G<tx. 

from which follows the impoTtant formula 

If now the operator /(D) be expressed in factors, we have 
/(D)jr = (D— o,)(D— a,) {D~a^\t/ ~ Q Hence 

(D— <xj) (D— = = 
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Also 7 . 

{D-<x^)...{D-cc^)y = J 

' _ ga.i I J ((?a;)-. 

By successive operations of this kind we obtain finally 
y = e“*® J e'tvi-oAc J,..J e“*‘=^Q (da:)". 

Second method. Suppose that, thinking of D as an algebraic 
quantity, we decompose 1//(D) into its partial fractions. Thus 
1 Aj , A.> , , A 

wr' 


D- 




D- 


D—ol,, 


( 1 ) 


BO that 

1 = Ai(2)— a2)...(f)— a„)+A2{i>— q:i)(A>— a3)...(I>— a„)+ 

+ ..>+An(Z) — <Xi)...{D Oji-i). (2) 

Each member of (1) may now be interpreted as an operator. 
The question is, have they as operators the same meaning or 
is the following equation true: 


This will be the case if 




(3) 


Since f{D) may be expressed in factors which may be written 
in any order, the right member of (3) is 

{Aj{D—ou)...(D—oc„)+A2(D—ai)[D—a3)...{D—cCj,)+ 

+...+A„[D—ai)...{D—(x„_^)}y. 

The polynomial in the bracket reduces to 1 by the use of (2), 
and by virtue of the laws governing the operator D this expres- 
sion has the value y. The particular integral now takes the 
formf 

e-“‘*Qda; -fAjC®** J c““«®9dx-(- 

J e-<=^=Qdx. 

t Tho fonnulao of this section may bo applied formally oven when the roots 
— 0 ore imaginaiy, Thorcsultingsolutions, however, involve functions 
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Examtlz 1 (ZP— 4Z)+3)y *= «** 

This equation may be WBttco (U— 1)(D— 3)y = e** By the 
first method 

<“-’)='■= 2CT'“ “ ''T “ T 

By the second method 

J ^,1/1 Ml. 

^ (Z)-l>(Z)-3)^ 2\2)-3 ’ 

Thicb gires the same result The complementary function is 
found in the nsu^ v&y 
ExaufleS (22)>+6Z)+2)yar6 
Example 3 (J5*— 4)y = * 

The general methods of the present article may be employed 
to determine a particular integnl, so matter what may be the 
function <3 in the nght member of the equation There are, 
havcTet, many cases la which a particnlar integral may be 
found more readily than by performing the integrations in 
reived m the genera] methods These cases will now be con 
sidered 


36 f{.D]y =a e" Smee Dtf^ =» oe", 2Fe“ = b*c", , 
= a"<", hence /(Z>)«" = /(o)e“ Dindmg this equation 
b> the constant /(a) 'which we assume for the present is not 


e" Since /{D) is commutative with 


with imsginarx sixum«ota Floin (bne functions real solutions ms)' be denved 
bj the methods of the lheoi7 of funrlions of « complex vansble solutions 

rnsv however bo obtsined directlr Tbusmlheequstion {(i>— a)*+^}y 
J{D) - 

The solution of this equstion nuij be obtsined from the lomuls 

This foncula tna; be veriSed dmttl^ It «t eoD'vementl; deiwed bj the 
mcthoil of tsmtion of psmmetCTS. Cf ert 61 
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S5 


a constantf this may be -written "whence it fol- 

/{®/ 

-c“. A particular integral of the given 


e“ 


‘o- O)- 

equation in case/{o) ,6 0 is therefore y = e'“//(ffi). 

If /(a) = 0 the polynomial f{D) contains the factor D—u, 
repeated, it may be, m times. Then/(Z?) = {D—a)”'<j)(D) where 
9 ^ 0. In this case <li{D){D—a)”^ = 6“. Hence 

{D—aY'y = = -r— r 

<j>{D) ^(o) 


and 


y = 


1 


1 


1 


{D-a)^ 4>[a) ^(a) (D-o)”* 

which may be evaluated by the first method of -art. 36 or by 
the method of art. 37. 

Example 1. (H-— 4)y = (l+c®J-. 

1 


A particular integral is 
1 


Z)=-4 


1 = 


H"— 4 
1 


(l-}-2e®+e®*). Now 


pOx 


Again 


1 2e* 

2e® = 


Z?=-4 


1—4 


= _ — _4 
J52_4® _4 - « 

— |e®. Also 


jD2-4 


X>-2D+2 


glx _ 




D— 2 4 


J da: = ^e^x. 


The complete solution is y = — ^^—i+Cxe-^-{-Cne-^. 

Example 2. {'iD^+10D—S)y = 

Example 3. {D^—a^)y — e'’®. 

t Tho fact that /(D) is commutative with constants is apparent. That 
is also commutative with constants may not bo so obvious. Suppose 

y *= Then f(D)y — cQ. Honco ^/(D)y = Q and /(L)| = Q whence 


2X 


I 
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37. f{D)y = €^r where F 1* a fuactloa of x We heve 
i)e"r = e~Z>F+oe-»r = e"{2)+a)F 
AlsoZ)*«"r = «**(2)*+2Z)a+o*)r We shall wnte the operator 
jD*+2i?(i-j-a* m the abridged form (f?+e)* and Bimilarly for 
higher po^rs of Z)-ro Then = «^(D+o)*F In the 

aanjewayitisfonnd tbat2>’e“*r = e"(2)+a)’r, and the general 
result C'ef^V = e*'{Z»+o)*]’ is easily established by mathe- 
matical induction Hi adding t<^ether terms of the form 
a^D^ef^T it follows that 


J{D)t?^V = t^J{D+a)V (1) 

The result of the inverse operation on the prodnet e“r 

IS now obtained by the following device Prom (1) 


Since F is an arbitrarv function of x f(D-^o)V is also an arbi 
traiy function Let it be denoted by /(Z)+o)r k ^ Shen 

This relation beu^ true for an arbitrary function we may 
now replace T1 by I Hence jr^e"F “ e” ^ — — -F This 


The exceptional case of art 36 may now be considered 
Suppose that F = 1 and that /(/)) = (C— a)* Then 

which disposes of this exceptional case 
Example l (i>*— 2i>^3)y = e** 


The complementary function is Cje“*4-c*e** For the parti 
cular integral (2)+l)(i?— 3}y = 
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PABTICtJLAB INTEGRAL OF /(D)y = Q 
1 

Hence {D~3)y = 

and y = ^ ^ ^ 

ExAjrPLE 2. (Z)-+4jD+4)y = e-^sina;. 
Exa3Ifle 3. (Zy*— l)y = c®. 

Exasiple 4. (D-—3D-\-2)y — xe^. 


38. SiP)y — sinoajor cosax. Proceeding as in the previous 
cases we find D sin ax = o cos ax, 

D'sinox = — a^sinax. 

It follows that 2>'sincx = ( — a^)®sinax, 
jE)®8inax = {—a'^ysiaax 
and generally Z)-”'smax = (— a-)”'sincx. 

Consequently, if F{D-) is a polynomial in D^, 

F{D-)sin.ax = P(— c^)sinax. 

Similarly F{D-)cosax ~ E(— a“)cosax. 

The results of the inverse operation 1/F{D^) on sin ax and 
cosax are now deduced as follows on the assumption that 


F{—a^) is not zero. 


jrE(Z>=)sinax = sinax. Hence 

F{—a-) 


sin ax . , sin ax 1 . -r 

F{D~) ■ j, - ; = smax and -— , -r = .^Tw^smax. In the 


P(-a-) 


same manner 


-cosax : 


F{-a^) F{D^) 

cosax 


F(jD^) F(—a-)' 

These formulae give a particular integral in cose /(£>) has 
the form F{D‘), i.c. it contains only even powers of D, and 
P(-a2) ^ 0. 

The method of finding a particular integral in case f{D) con- 
tains odd ns well as even powers of D is illustrated in the 
following example; 

(IP-fX)--fi?-}-6)y == cos2x. 


«15 
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A particular integral « 

(the juBtihcatioQ of thia step u left to the student) 

= D»-6) 

= -i(i)»+D_J|)»-6)coa2r= -A(Bam2r-2cosar) 
The foUomog vanation of the method u sometimefl shorter 

7-»(-i+ir-<+«“““ “ =55+5“”“* “ 

In the exceptional case n whl^ f(— a®) ^ 0 the foregoing 
method fails to give a particular integral A method ofhandliDg 
this case u lilustrated m the following example 

^4.oV*oc«B» (1) 

If the right member be replaced by cos(a+h)r, It ^ 0, the egua 
tion u no longer exceptional, 

— «>e(o+^)* (2) 


A particular integral of this cqaation is 

co8(a+Jl)* oos(aH-A)z 

D»+o*'' 


~coa{a+A)r = 


■ -(«+A>*+o» -3oh-h* 

By use of the theorem of mean ralue of the difiecential calculus 
•we find 

co8(a+A)r = cosox— &raui(a-(-0A}x, 0<^<1 

Hence the particular Intend of (2) becomes 
cos ox— AzBiD(a+$A)x 
— 2aA— A* 
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The term — is included in the term A cos ax which is 
^2ah—h- 

part of the^ complementary function of (2). It may therefore 
be omitted from the particular integral and we obtain the new 
particular integral 

~7ixmn{a+6h)x —x^{aA-07i)x 

— 2a/!— h- ~ ~2a—h 


This is a particular integral of (2) for all values of h other than 
zero. If we assume that the solution is a continuous function 
of h when h, = 0, we find a particular integral of (1) to be 

O* RlTl 

— . This result is embodied in the formula 

2o 


1 


J)-+a^ 


COBOX = 


a: sin a* 
2a 


wliich may be verified directly. The student should develop for 
himself the corresponding formula 


Exajiple 1. 
Exasiple 2. 
Example 3. 


rsmaa: = 


xcosaa: 


D-+a- "" ■ 2a 

{IA—5D--{-4)y = sin^x. 

{D--\-D-\-l)y — cosx. 

(X) 2 -|- 4 )(i) 2 _|_ 9 )y = 20sin2x— 30cos3x. 


39. f(D)y = x"*, m a positive integer. We consider first the 
simple case (D—a)y — x”'. A particular integral is 

1 r 

y = a:"* = e"* I e~°*x’" dx. 

D — a J 

Integration by parts gives the formula 




oxaf clx 


3fC~'’ 


-- f 

o J 


■ dx. 


By successive applications of this formula we find 

n-a^ 

m{m — 1) x”'~"e~“ 


a- 


w(ot— 1)...2. ] e-°^ 




a 
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It IS important to observe that this same result is obtained by 
expanding the operator by the binomial theorem and operating 
on i" by the successive povrete of D Thus 




, w(ff>— 1) 2 r 


Consider noir the general cose m vhicb y a 


m 


*“ If the 


operator be separated uto its partm) frwtions the fore 
going method mil apply to each faction The reo'dt is eqtu 
valent to that obtained by expanding in a senes of 

pouers of D and operating on z" by each term of the senes 
Since i>z* » 0 vheR r > m the result mil be a finite zrumber 
of terms 

Exawxe 1 (D*— 2)+l)y * 

Tbe ccTOplementary function is e*'(j4coslV324-B8in{V3z) 
A particular integral is 


1 




= )r»»=a:»+3*»-6 

Example 2 (ZI*+4Jy = 3x* 

Example 3 (Z)»-3Z)*+2D)y = 12(r*-2z+4) 


EXAMPLES ON CHIlPTEB VI 

1 (l>‘+x>’+l)y “ lb 2 (I>*+2tJ+I)jr= ITcoeAr 

3 (3Z)»+P-U)jr= 13f» 4 (0‘-X>»-D+l)y = !r 

6 (I)*-®I>+. 10 ^ = a»<*» e (D»-Sr>+Jjy = ** 

7 (I>»+3D*+31>+S)y = 

8 l£>-oyy = «" 

9 (D‘+26X>»+S5)y = 24ginS*eo»3i 



EXAMPLK, fil 

10. (D-a)^y = a^'c". II. 'd+2)'(D+3)i/ = c-«. 

12. {D'+5D’+4)y == 128in2ar. ’ 

13. (6DHn2?-10)y = 19(el*+c-t*) 

14. {D^-2D+l)y = xc^-sinar. IS. [D^-9)y = 3— 9x*+27x‘. 

16. {D^+5)y = sinx+26in*x. 

17. Show that if P is a function of a? and the operator /(D) has the 
form of a polynomial, 

(i) D^ixY) = xD'>r+TiD"-iF. 

(ii) /(D)(xF) = x/{D)F+/'(D)F. 

18. (D>+4)y = xcosx. 19. (D— l)y = 2x-sinx. 

20. (D>+l)y = 4a;cosx. 21. (D*+2D»+l)y = 8cosx. 

22. (D’+2D5+D)y = 4x’-CiA 

23. (D*+4D+8)y=: 12c”“sinxsin3x. 

24. {4DH8D+3)y = c-^Cx’+sinix). 

26. (D*-+D‘)i/ = 4+12xH4sinH:c. 

20. (3D»+30)y = 10co3*x. 

27. The differential equation for .simple harmonic motion is 

^ + a'^ = 0. 

Show that if the particle starts from its equilibrium position with speed 
V, its displacement ( seconds later is a = -sinal. 

28. The differential equation for harmonic motion with a damping 

d^8 ds 

proportional to the velocity is ^ + 2J5^ + a's = 0. Solve the equation, 

and discuss the character of the motion for different values of a and p, 

29. Show that a particle moving in a straight line according to the 

. d'3 

equation ^ + cos at oscillates with a constant period and ever- 

increasing amplitude. 

30. The motion of a particle on a straight line is given by 

"S + 2|+40a = 0, 


where s is measured in inches and ( in seconds. Show that the particle 
executes oscillations which die down to about of their original ampli- 
tude in 15 seconds. 



CHAPTER VII 


METHODS APPLICABLE TO CERTAIN LINEAR 
EQUATIONS 171111 VARIABLE COEFFICIENIS 
Wnzv the coeScjento u> a Lnear equation are not conalant no 
Bucli simple and eflective methods tor its solution am available 
as in the case of constant coefficienta The first type of equation 
to be considered in the present chapter transforms by a change 
of variable into an equation with constant coclficKnta and may 
therefore be solved by the methods of the preceding chapter 
Next follows a discussion of the rare case in which the linear 
equation IS exact Thirdly a description is given with examples, 
of the important method of obtaining solutions in tlie form of 
infinite senes 


40 The bomogeneous linear equatloo The so-called homo 
geneous equation has the fimn 






in which the a s arc constant and Q is a function of x We 
shall change the independent variable ar to e where s se e* or 
i-log, Tten 

dz dzdx xdz 


dx^ 


dr z^lds* dzj 


Also 
and 

<fy 

If the symbol D be used for ijdz these results may be wntten 


ilz}'^x\d3? djfds d-j 


.S^y 

dr* 

<Py 


Dy 

= D(0-l)s 
= i)(0-l)(D-21y 
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HOMOGENEOUS LINEAR EQUATIONS 
These results are readily extended to higher derivatives. The 
general formula, which may be established by mathematical 
induction, is 

By these substitutions the given differential equation is trans- 
formed into an equation with constant coefficients and the 
methods of Chapter VI are then available for its solution. In 
particular it follows that the complementary function is found 
by solving the auxiliary' equation, 

jn(m— l)...(m— l)...(}n— n)-l- 

+-+«fn-l”»+«n = 0. 

If Wj, are roots of this equation and are all distinct 

then the solution of the reduced equation is 

or, in terms of the original variable x, 

...-{-c,, a:’”". 

If is a double root of the auxiliary equation the corresponding 
terms in the complementary function are e’^'^ic^+c'^z) or 
a:"'(c,.-fc'loga;). 

Example 1. — — 0. ^ ^ 

ax- dx 

Example 2. x -^^;-\- Gx^-\-\Gy = 0. 

(vX** uX 

Example 3. a ?^^— x ^-\- 2y -— 0. 
ax^ ax 

41. The particular Integral for the homogeneous form. 
When the homogeneous equation is transformed into an equa- 
tion with constant coefficients we may, as in Chapter VI, use 
the symbol D for djdz. It is possible, however, to avoid explicit 
reference to the variable z by using instead of D a new symbol 
6, so that 
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The first-order equation — Q solution 

s= 2*] ar~* from wlach follows the formula 

If then the general equation expressed m terms of the operator 
9 has the form J(0)y = <? » particular integral may be found 
by factoring F(6) and employing either of the two methods of 
art 35 

Each of the rules which were worked out la Chapter VI for 
wntmg particular integrals corresponding to special forms of 
the ngfat member Q has its coimteipart in the case of the homo- 
geneous equation The only case of sufficient importance to 
dfseire mention is the follofwing The equation /(6)y ** i* 
becomes with the independent vanahle 2 = logar f(D)y e** 

t*** x* 

A particular utegra! u therefore y » e> yjjm) 

Etm) ^ b Example \ shows how to deal with the exceptional 
rase in which fim) ss 0 

Ex-u.n»L8 1 **«— *» ** 

ax* ax 

In terms of the operator $ this equation is 

{0(B-l)-.2B+2}y ^ X* or («— lKP-2)y * ** 

When * «= e* this becomes (f)— 1)(Z)— 2)y = e** Hence a par 
ticular integral u 

(Z)-l)(Z)-2) U>-2 H-1/ *«- e- g l 

The complementary functwa IB y = CjX*-f-c,i and the complete 
solution y » a:*(log»+Cj)-l-c,* where Cj = ej—I 

ExiiMFLE 2 2*^-f-6x^-|.6y= e* 

Exasiple 3 *»^-{.6jr*?|+8r^-t-2y = **+ 3a:-4 

<ir* djt* dx ’ 

\ 

12 Exact linear equations If from an equatwa *nf order 
j— I containing an arbittart constant an equation of order n 
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is derived by differentiating and eliminating the constant, the 
former equation is called a first integral of the latter. An equa- 
tion is called exact if it may be derived from its first integral 
by differentiation and no further ojjeration. Exact equations 
are of rare occurrence in the applications of the subject. In case 
they are Unear they may be recognized and the first integrals 
written down by the rule which follows. If then a first integral 
is exact, a second integral may be obtained in the same way. 

The linear equation to be tested is 






da:"‘ 




( 1 ) 


If the left member is the derivative of a differential expression, 
of order n—l the latter must contain the term — r- This 
term gives on differentiation the first term of (1) and also the 
term To cancel this latter term and to give also tlie 

second term of (1) the first integral must contain (Pi— P o)t — 
This term gives on differentiation the additional term 

fP' 


This term will be cancelled and the third term of (1) provided 

for by a term (P„— P^-fPJ)^j^ in the first integral. By this 

procedure we find that aU the terms except the last in the left 
member of (1) are obtained by differentiating 


p 




d"~^y 


(P-Pi-fPS) 




The derivative of this expression includes also a final term 
Pn- 2 +...(— l)"“^P{,">)y. Then (1) will be exact if and 
only if this final term is identical with P„y, that is, if 

K 


‘ oil 
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In case this condition is eatisfied a fint integral of (1) is 




+ +(!■. .-n 1+ = 

In writing out a first integral for a gircn exact equation it 
IS preferable that the method rather than the formula be used 
The method xaaj in fact be used for other than bnear equa 
tioos though no convenient test for exactness can be given for 
the general case 

For a linear equation whi^h is not exact the question arises 
of whether it can be made exact by multiplying it an 
integratmg factor The following discussion for the second 
order equation might be extended to equations of higher order 
Let the equation 


be ffiultipbed by a function ^(x) Then 

isesutif 


p.&j. 


dx )dx \d; 


Any function fi{x) which satisfies fZ) is an integrating factor of 
(i) As (2) IS generally as difficult to solve as (1) this fact is 
seldom of assistance in the solution of (t) The relation between 
the two equations is however an interesting one The student 
should verify that a solution of the reduced equation of (1) is 
an integrating factor of (2) On account of this reciprocal rela 
tion between (2) and the reduced equation of (1) these equations 
arc called adjoint 


ExAMfLE 1 + = 


The test for exactness is satisfied since 6—18+18—6 = 0 
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EXACrr LINEAR EQUATIONS 
A first integral is 

3x2)^+(18a:— 12x)i/ = sinx+Ci 
da? ax 

or a?^+6z2^+6a:y = smx+Ci. 

OX- ox 

This equation is also exact, giving the second integral 

s?~-{-Z3?y — — cosx+CiX+Ca 
dx 

■which is again exact and gives the final solution 

a?y — — sinx+CiX^+CjX+Cj 

where = ic^. 

Example 2. (ar>+2a^)^+2(5x®+8x)^+4(5x-f'6)y = 0. 
ax~ ux 

This equation is not exact. Its adjoint is found to be 
(i«>+2a?)g-(4*>+te)&+(6*+12),. = 0 


which is of homogeneous type and has the two independent 
solutions /i = X- and n = sp. These are two integrating factors 
either of which makes the given equation exact. 

Example 3. 

«='{^-4)g+12x=(x-3)g+36x(x-2)^+24(x-2)y = 0. ' 


Example 4. 




-(x-"-l)y = 0. 


43. Solutions in series. For such diiferential equations as we 
have been able to solve the solutions have been expressed by 
the aid of the so-called elementary functions. The usefulness 
of such functions as the trigonometric, logarithmic, and expo- 
nential functions is due partly to the fact that their properties 
arc well known but also to the fact that we possess tables of 
their values which make it possible to obtain numerical results 
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readily in problems whose BcJntions inrolre these functions 
These numencal tables are moat easily computed by the aid of 
infinite senes These senes, moreover, have other important 
uses, if we start from them as de&utiooa the properties of the 
functions may be denved For example, 

. ** 

3i+5i- ■ 

and from this senes a table of values of tbs sme may be com 
pated and the science of elementary tngonottetry may be 
built up 

Other functions most which are not met with in elementaiy 
mathematics but which may be represented or defined by in 
finite senes When such a fonebon haa been denoted by a 
symbol a table of its values computed, and a catalogue of its 
properties denved from the eenea definition, it takes its place 
in the family o** transcendental functioiis on a par with suix, 
sin 'z e*, logz, etc In most cases these sew traosceadentol 
functions arise as the solutions of differential equations SIo» 
over, an infinite eenea may be the only practicable way m which 
to express the solution of a diffete&u^ equation even though 
the function represented is not of sufficient importance to de* 
serve a name or a symbol In case we find no combination of 
elementary functions in finite fonn to satisfy a diffierential equa 
tion we therefore endeavour to find an infinite aenea which 
converges for some interval of values of the independent variable 
and satisfies the equation As a fint illustration of the method 
we shail find the solution u senes of the equation ^ ^ 

Assume that this equation can be satisfied by a senes of the form 

It a proved in the theory of infiiute eenea that the denvativea 
of s function represented by a power aenea withm its interval of 
convergence may be obtained by differentiating the se:^ tenn 
by term Consequently 

^ = 2a, 4-2 3a,x+3 4a, ^+4 50,**+ 
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69 


Hence , 


p,+y = (2a2+Oo)+(2.3a3+ai)x+(3.4a4+a2)a^+... . 


In order that the equation be satisfied this last series must be 
identically zero. A power series can vanish identically only 


when each coefficient vanishes. Hence = 


«i 


31’ 


a — a — — 

0“4!’ 4.5 




Similarly, Cg = — 


fli Qq 

= ~7i’ ~ §1’ ’ 


9! 


i , etc. Hence the series which satisfies 


the differential equation is 


'3! 


+ 


\ 

5! ■"}■ 


It is readily established that each series in parentheses con- 
verges for all values of x. Since the equation is satisfied what- 
ever be the values of Og and a^, these are arbitrary constants 
and wo have found the general solution of the equation valid 
for all values of x. This solution, of course, we recognize as 
y — Og cos x-j-Oj sip X, " 

Consider next an equation whose solutions are not elementary 
functions, x^-j-^-{-xy = 0 (Bessel’s equation of the zeroth 
order). As before, we assume the series 


firom which 


y — Co-fOiX-fCza^^-hOsa^^-f ... 


d« 


Ci-j- 202 x-J-SOg x--f 40^ X® . 


and 


dx2 


200-1- 2. 303 x-f 3. 404X® • 


We multiply the first and third series by x, add the three 
resulting series, and equate to zero the coefficient of each power 
of X. This gives Oj = 0, 2 o,-l- 2 o„-)-Oo = 0, 2 . SOg-l-Sag-l-ai = 0... , 
and generally 

n(n-f l)a„+i-f (7i-}-l)o„+i-}-fl„_2 = 0 or (n-f-l)=a,.+i = -o„_i 
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from which o„+, *= — It foUowB that 

s a, s: z= c 

•nd th.t «. - -g, ». = «. - '<» 

j Thesenesinpmn 

theses, which is easilj shown to eonrerge for aH values of x, is 
known as a Bessel funetwn of the seioth order and is denoted 
hj «^o(^) -As y =3 contains only one arbitrary constant 

it u not the general solution of the equation Another mdepen 
dent senes solution may be obtained in a less simple manner 

44 Cootinuatloa. The Indlclal equation. It may happen 
that a solution of a differential equation cannot be expressed 
in a senes of positive mtegral powers of x For example, the 
frmctious log z, •h, wbieh may occur in solutions, have so 
expansions of this form Certain differentia) equations may bd 
satisfied by senes uvolvug negative or fractional powers of sr 
AdiscussionofaUcaseswillDoibeattemptedbeie Weehallcon* 
fine our attention to equations of form ^ (1 

in which P and Q are rational fanctiona of x For such equations 
we shall lUostrate by ezamplea a method of wide generality, 
leaving aside the finding of a second independent solution ifhen 
only one is furnished by this method The method consists in 
assunung a senes of the form y ^ z^(ag-f'(iix+<*iX*+ ) und 
attempting to detenzune the index r and the coefficients so as 
to satisfy the differential equation. 

EsaMPLS 1 {2z+z*>^+^— 2y =/ 0 
Assume y SB o,z'+Ojaf**+Ojaf**+ where o, ^ 0 Then 
^ « ra^x' *+{'-+l)<Haf+<r+2)a,af+‘+ , 

g = {r-J}«i,3f-Hr(r+IJa,x''-H(r+l!(f+2)<i,zr4-’^ 
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COimMATIOK. THE lOTICIAI. EQUATION 
the coefficients of po^ /2fr— lir+rVio = 0- ^ 

2r(r+l)ai+{r-l)rao+(r+l)®i-2«o 0> 


firom which 
Similarly, 



ft r, 

®1 — ■ 

2r+l 


— - d. 

02 ■“ 

2r+5 

03 = • 

L_a. 

2r+6 ■ 


etc. 


UsinK first the value r = 0. we find % = 2ao, 

— fl =0 = ... = 0, and the corresponding solution reduces 

to a polyno'mial, y = Oo(l+2*4-|a^) = «o«. say. Take next 

.4 .1 „ _ 1 


or - - 

the value r = Then 


-3 

«i = 


•0i 


—1 


-3.-1 


4.8 




On - ^2^2 


Similarly, 

- 3 . - 1 . 1. 3 

04 = — 

Hence 


y 


4.8.12.16 

‘h 


06 = 


-3. -1.1 
4.8.12 “ 

-3. -1.1. 3. 5 


4.8.12.16.20 
3. -1.1 


: OnXM 


-3 


4.8 


etc. 


4 8 12 




= OqU, say. 

It is easily found by the ratio test that this series converges for 
values of x numerically less than 2. Hence for these values of 
a: we have two independent solutions. In each case the constant 
Oq is arbitrary. If we call it o in the first case and 6 in the 
second, the general solution is y = ow+6w. 
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Example 2 (1— 6y = 0 (Legendre e eqoa 
tion of Older two) 

Proceeding to satisfy this equation by a senes we msjr set 
down the work as follows 


- 2 r 


1 -x* 


y=eBa'+o,a 5 +i+o, 3 ^**+ 
dx 
dz*" 


^ = ra,ir^+(f+l)a,jr+{r+2)a,ifH+ 


*= (r-l)fa,*' »+r(*-+lK3^ H 

+(r+l)(r+ 2 )a,z'+ 

The indicia} equation is f(r— 1) = 0 whence r = 0 or 1 Next 
r(r+l)ui =* 0 Hence if r = 0 o, is arbitrsiy and if r = 1 
ss 0 Next, 

^ (f+ 2 )(r+ 3 > ‘ 

rjr+i) 


Now if r 


-1 * 

* 23 ' 


(r+SlP'+^J**’ 
i«, o,»= 0 =:o,«a,* 
-1 1 4 e_ 


- 1 X 541 


y = Oa(l- 3 x*)+o,|*+- 


2 4 3 5 
~«j etc 


:U-j-V+ 1 

J 4 3 6 ^ j 


Again ifrs=i ajss<[,Esa, ssQ and a, = 
—114 6 . „ 
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It now appears that this series is already contained in the 
solution for r — 0. Hence for (a:j < 1, which is the interval of 
convergence of the series, the general solution of the equation 
is obtained from one root of the indicial equation. 


Example 3^ 


Obtain the indicial equation for 


da^ 


+ 1 / = 0 , 


Show that the general solution is given by one root and that 
the solution given by the other root merely repeats one series. 

Example 4. Show that for Bessel’s equation of order zero, 
the roots of the indicial equation are equal and only one solution 
is furnished by the present method. 


Example 5. 




o.‘ 


Example B. 2x(l—z)~-}-(l—9z)~—3t/ == 0. 
Example 7. — 0 . 

( l 4 C '‘ a-X 


45. Equations of Bessel, Legendre, and Riccati and the 
hypergeometric equation. In illustrating the method of in- 
tegration in series we have confined attention to equations of 

the form = 0 w'here P and Q are rational func- 

dx- dx 

tions of X . Included under this form and solved by this method 
are tluee equations of importance in various branches of mathe- 
matics. They are 


(i) Bessel’s equation. 




0, 


(ii) Legendre’s equation, 

L 


«12 
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(uj) the hjpet^oznetnc eqaabctn. 

In the case of (i) and (u) the character of the solntions depends 
on the constant n, called the order of the eqnatjon (not to be 
confused with the usual meaomg of order of a differential equa 
lion) The aolutiona of (i) with suitable constant factors are 
called Bessel functions of order n (see also art 77) The indicia) 
equation gives r « <^d leads to solutions denoted bj 

and </-a(z) When a is not an integer these solutions are in 
dependent and the general solution is y = AJ^(x)-\-BJ^{x) 
IThen n is an integer J_a(a) turns out to be a multiple of J^[x) 
and then a senes solution of another sort f called a Bessel 

function of the second kind, u necessary for a complete solution 
The Bolutiona of (u) with auitabSe constant facton are called 
Legendre functions (see also art 77) If n is an uteger, one 
solution reduces to a polynomial (as la ex 2 art 44) This 
polynomial is denoted by when the constant factor is 
chosen so as to satisfy ths condition 1^,(1) » 1 Tbns the first 
four Legendre polynomials are 

|x*-i. iiW «{**-!* 

The usefulness of integration m senes is not confined to 
equations of the form considered above An example of another 
sort IS 

(it) Riccati a equation 

^+6y* = ci" 

a non linear first-order equation whosd'solutions must m general 
be expressed m infinite senes. Tfaesolutioivof Biccati t equation 
may m fact be expressed m terms of Bes^I functions 


EZaUFLES ON CHaiTEB m 


I x> 


i<Py_ , ^ , 


t See B O Peirce A ■SAerf ToUe^/nle^rab, formulM 73S 738 
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3. {i‘-4ar")^+(8a:^-24x=)^ + (12a:»-24a:)i/ ='12x*. 

6. Bj’ the method of series find one solution of 

■ 6y = 0. 


2x^ + 3xJ^. 

ox* dx 


7. Obtain one solution of Bessel’s equation of order n. 

8. Solve in series Legendre’s equation of order one. 

9. 2(ax+b)*^+ 7a(ax+6)^-3aV = a^- 

10. Show that the following equation is exact and solve it: 

(3xV-2x2/)g + (ex--i/-2x)(^)*+(12x2/*-4i/) J + 22/= = 0. 


dx* 


dx 


Vx. 


12. Show that Riccati’s equation is reduced to a linear equation of 

the second order by the transformation « = ^ 

o: dx 

13. Obtain the logarithmic series as a solution of the differential 
equation {l+x)g + J= 0. 

14. Obtain the series for tan~*x os a solution of 





CHAPTER VIII 


SnSCELLAKEOUS iCETHODS FOR THE SOLUTIOX 
OF EQ^7ATIO^S OF THE SECOND AND 
HIGHER ORDERS 

Is this chaptei are toUected some methods of aolvmg dif 
ferential equ&tioas of certain special kinds In each case some 
peculiarity of the equation or some partial knowledge of its 
solution furnishes a clue to the method employed. It may hap 
pen that a particular method u useful only m finding a first 
integral of an equation leaving the further reduction to some 
other method In what foUova x icQl be used for the mdepen 
dent and y for the dependent ranahle unless otherwise stated 


46 Eqoatioas from wblcli y Is missing An equation of 
order fl from which y is missing may be symbolized as 

!')=» 

This 13 equivalent to an equation of order n-»l m the dependent 
variable p =* dyidx If the aolntiou of this equation can bo 
found in the form p » f{x) then a further integration finds y 
Again if y and its first r— 1 denvatires are mixsiBg the sub> 
sUtution d'yjdxf = q depresses the order of the equation by r 


OaBubetitutingpfor^wehaTeaj^+ZMj*— p «= 0 Trhence 
oz az 

q- 2z dr = 0 which has the solution — = c or 
P 

p = — — llus la turn baa the solution y as }Iog(**—e)4-c, 


EhjimJ 
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EQUATIONS FROM WHICH y IS MISSING 
Example 4. a:~+^ = 4a;. 


47. Equations from which x is missing. In this case the 
substitution of p for dyjdx also leads to an equation of lower 
order, hut the higher derivatives must be expressed with y as 

independent variable. Thus g = g = 

etc. The method is chiefly useful in the solution of second-order 
equations. 


= 1 - 

Example 2. 2j,g-(|)'=l. 


Example 3. 



^dxdx^^\dxy ■ 


48. The equation d^jdz^ =f(y)- This equation, which is the 
form assiuned by certain differential equations of mechanics, 
comes under the type discussed in the preceding article and may 
he solved by the method there given. An alternative mfethod 
is to multiply both members of the equation by 2 dyjdx, thus 
making it exact. Then 


2 ^^ 
dx the* 



Integration with respect to a; gives the first integral 
(|)’ = 2 //&)*+«. 

The solution of this equation offers no difficulty apart from the 
actual integration. 

The equation d’^yjdx" =/(d"-*y/da:’*~*) takes the form just 
discussed on the substitution d’^-^yjdx”-- = v which is therefore 
a first step towards Its solution. 

Example I. = o. 

dx- 



18 EQUATIONS OF SECOND AND HIGHER ORDEl^ 
ExAMPM 2 ^ 3a*j^, with the initial condition ^ ® 

when y »= 0 

ExAjm-E 3 2^ =» e*', with the condition ^ = 1 when 

dr* dr 

ysO 

Example 4 =» 2y*— y, with tl» condition ^ ® when 

y = o 


49 The linear equation of the second order. General 
solution found from a known solution of the reduced 
equation In the solution of the linear equation 






( 1 ) 


a method which is aomehmes effeetiTe u to replace y by ^be 
product of two rsnablea, the object then being to take for one 
of these rambles such a function that the resulting equation 
may be solved for the other one Thus, if y (1) becomes 


tion of the reduced equation ( 
tituted for « Hien^+rg 


If now any solution of the reduced equation of ) be knOTTn, 

let this be substituted for « Hien ^ end 

dx* ox 

(2) becomes 


which IS a linear equation of the first order in the ranable 
p = dv/dx If V 18 the generJ solution of (3), then y «= it« is 
the general solution of (1) It la to be noted that in order to 
use this method the complete complementary functiou of ( 1 ) 
need not be known but merely any integral belonging to it 
Such an integral may eometimes be noticed by an inspection 
of the equation 


Example I a;^— (3i+l)^+(2x+l)y = e»» 

The left member of this eolation reduces to zero on the sub 
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LINEAB EQUATION OF SECOND ORDER 
stitution y — eF. (It is clear that this will happen for any linear 
equation in which the sum of the coefficients of the left member 
is zero.) Hence substitute y — d^v. The equation becomes 

cPv 1 

(fo:* \ xjdx X ’ 

whence 

— = r — — Ac = xd^i — -+c] = cxe®— e® 

dx \J xd^.x } \ X ) 

and V — c(xc®— e®)— c^+Cj. The solution of the given equation 
is accordingly y = c-®(cx— c— ll+CjC®. 

Example 2. x-^— (2x34-3x)^-}- (2x^+3)y = x®. 
dx~ dx 

(»/ = X is evidently a solution of the reduced equation.) 

Example 3, x-^— (x-4-6x)^+(3x+12)y = 0 given that 
dz’’ ax 

i/ = ar’ is a solution. 

Exaotle 4. (2a:+3)^+ {2x+2)y = xe-*. 


50, Reduction of the linear equation of the second order 
to normal form. In the preceding article equation (2) was 
reduced by choosing for u a function to make the coefficient of 
V zero. Another way of reducing (2) is to choose for u a function 
to make the coefficient of dvjdx zero. This function will be a 

solution of 2^-f-P« = 0 from which u — only a par- 

ticular value of the integral being necessary. To substitute this 
value of u in (2) we have 

and (2) becomes 


form, which may be abbreviated ^4-7u = 5r, 

Ac® 
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the nonnal form of a linear equation of the second order It 
sometimes happens that the nonnal form of an equation may 
be solved by one of the methods iutherto considered In that 
case the solution of the ongmal equation is ohtamed by multi 
plying by the factor «-*!*“** 

ExamplsI ^+2*^+(»»4-5)y= he*-*** 

Here / = ^ « 4 

and e*!^*** =s e**' The nonnal form of the equation is 
d*e/<fe*+4e s= Oe* •**«•** =s 5e*, from which 
tl 5= <*+«co8ieH-Asin2r 

It follows that y » „ e't**(eV4-acos2z+68m2x} 

EXAS£PtB2 4z^-4(K+2)^+<4+ar)y« 0 
EXiMTLE 3 (**-2r*)^+2i*^-12(-r-2)y 0 


51 Varlatloo of parameters In case the complementary 
fimctioD of a linear equation is known the following interestmg 
method of findit^ a particular integral may be used The 
method which is due to Lagrange, consists in replacing the 
arbitrary constants of tbe complementary function by variables 
and determimig these variables so that the given equation is 
satisfied The solution of example 1 shows the details of the 
method and it is clear that the method apphes to linear equa 
tiona of any order with coewtaut or vanabte coefficients 

ExAJtPLE 1 ^-[-(Z*y =: sec (»-l- CSC oz ( 1 ) 

The complementary function is Mcosoz+.Bsmax Take A 
and B as functions of z to be det«mined and set 


Then 


y = iicosox+fsinax 


iy 

dx 


= — o^smiu-fuBeosax+uosQz 


U 

dx 
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As a first condition for determining A and B we shall place 


dA . . dB 
cos ax— — f-sm ax -r- 
dx dx 


0 . 


( 2 ) 


Then 

Again 


^ _ 

dx 


-a A sinaa;4-a5cosffX, 


. . dA , dB 

-A = — aMcosaar— a-^smax— asmaa:-— d-acosax— . 
dx- dx dx 

Substituting these values of y and d-y/dx? in the given equation, 


we have 


dA , dB , 

— asm ax — +a cos ax— = sec ax-j-csc ax. 
dx dx 


But from (2), ^ = — cotax^, whence 
dx dx 


d_A 

dx 


asi 


cos-ax\ , 

smax— a— : 1 = secax+cscax 

V Binax/ 


dA 

dx 


= — i(tanax+l). Hence A = xV Also 

o a\ a } 

^ = i(14-cotax), and B — A particular in- 

uX Ct €L\^ O' j 

tegral of (1) is therefore 

1 lU * 

y = — (cosax logcosax -[-sin ax logsinax)— -(cosax— sin ax). 


Example 2. 


d^ 


1-fe* 


Example 3. x^^— Ax’^-fCx^ = given the com- 

plemcntary function ax^-\-bxA. 


52. Change of the independent variable. We have seen in 
arts. 49 and 60 how a change of the dependent variable may * 
bo useful in transforming a linear equation of the second order 
into an integrable form. Sometimes a change of the independent 
variable may be determined to achieve the same object. Thus « 

uo may transform the equation = J? by 

dx- dx 


■4212 


M 
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replacing * by z where z u a fosction of i to be chosen. Then 

dx dz dx 


, d*y d‘y/«fc\* rfy i*s 

iz* ’^di <Zi* 

and the gives equation becomes 

It may be that a valoe of z u (enna of x may be chosen so that 
this eqoabon takes a simple form Due possibibty is to choose 
z so that the coeScient of dyfdz ahall disappear Then 

from which 2 u found. 

Example 1 

Cbaegug the variable c to a, ve get 

dy/2 dz . d*2\ . 4 „ ^ 
dt*\dx) '^dz\^dx^dx*r3*^ x* • 


sdx «t** ' ' d* »* 

, ^ , 


The equation becomes, after dividing out z*, ^4-4^ = 


whence ys:z A cos(2z4-oj 


Example 3 4x^4-2^— tt*y » i* 


EXAMPLES ON CHAPTEB VTH 


■ s-ir- 


2 (*^+I)55-5{*+3)2| + rx+5Jy = e* 

3 dhf/dx* = *m 2y with the coitdiiion dy/<fc = 0 when y = 
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4. — (2x»-fx)^+(2x+l)y = i*. 

C. 3^~+3^^—2[x’+l)y — 0, given that y ■ 


K* is a solution. 


6- (®+o)||-3^= (a+a)’- 


dx 


7. ar*^+a*^+(a’— 2)y = 0, given that \jx is a solution. 

ox* ox 

8. 8inx|^— ^^+(sinx-{-co3x)y = c*eosx, given that sinx is an 
integral of the complementary function. 




11. By clianging the independent variable, solve the equations 

(i) = 8x5 sinx’, 

(ii) (l+x*)*g+2x(l+x’)^+4y=0. 

12. Sliow that a^+(2— a)^— y msy be expressed sj’mbolically in 

the factored form {D—l){xD+l)y but that the factors may not be 
viritten in reverse order. Use this method to solve the equation 

ag-f(2-x)|-y = e*. 

13. By factoring the operator, solve 

(i) (x-f 3)g-(2x+7)^ + 2y = x+2, 

(") (*+l)g-(a+3)2+2y = (x-f 1)’. 



CHAPTER rx 


APPZJCATIO%S OP EQt7ATI0^S OF THE 
SECOND ORDER 

Is Chapter V were considered eoine appbcations of differential 
equations of the first order If we allow out equations to contain 
second denvativea as well as first the range and variety of 
appUcatiooa are greatly increased la geometry for example 
we may express by means of second-order equations properties 
of a curve involving eurvalHre In mechanics many problems 
involving motion with known acceleration or the equilibrium of 
a system of forces may be formulated and solved by means 
of second-order equations 

53 Problems Involving the curv'ature o! a curve Tbe 
curvature of a curve is given m Cartesian coordinates by the 



and tbe radius of curvature la l/x Any relation among quanti 
ties which depend on the curvature tbe slope and the co- 
ordinates at an acbitcsry pouit on the curve may be expressed 
as an equation of the second order For a relatively small 
number of such equations tbe solutions may be obtained by tbe 
methods of the precedtug chapters 
ExASiPU 1 Among all tbe curves for which tbe curvature 
IS proportional to the lei^h of the normal determine those 
which are everywhere concave towards the * axis and cross it 
at right angles 

The length of the normal is given by the oamencal value of 
SimJe the curve w to be everywhere concave towards 
tbe X axis it follows that if y is positive d is negative 
and if y IS negative d*y/d*^ m positive Hence the differential 
equation 13 — = I yV(l+P*) Since X is missing from 
this equation we write itm the fona 
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the cur 7 e is to cross the a:-axis at right angles, then, ■when y 
approaches zero, p becomes infinite. Hence c = 0 and 



or p- = 


1-kY- 
khf • 


Tliis equation, with variables separated, becomes 


kiidy 


= dx 


from which — - = x-\-c or (a'+c)'4-i/- = 1/i®, which 

K 

denotes a set of circles of radius l/k with centres on the x-axis. 


ExAjrPLE '2. Determine the equation of the family of circles 
of given radius from the condition that the radius of curvature 
is constant. 


54. The catenary. If the ends of a uniform flexible chain are 
fixed at two points, the curve in which the chain hangs is called 
a catenary. The name is properly applied to an open curve 
extending upwards indefinitely, and such that if any two points 
on it be made the points of 
suspension the lower portion 
of the curve is the form in 
which a chain of proper 
length will hang. Since the 
curve has clearly a vertical 
axis of sjTnmctrj' wo take 
this as the y-axis. The posi- 
tion of the origin ■will be left 
undetermined for the present. 

To find the equation of the 
curve we must translate into 
mathematieal language the statement that anj’^ portion of the 
chain is in equilibrium under aU the forces which act on it. 
Acting on a segment PQ arc three forces, the tension at the 
tension at (?, and the weight of the segment PQ. Denote the 




S6 APPLICATIONS OF SECOND ORDER EQUATIONS 
pojnt P by (x y) the tension at J* by T, and the inclination of 
the tangent by t The eorrespondu^ quantities for Q differ from 
these by increments as shown in the figure If the arc PQ is 
denoted by &3, the weight of the segment is wAs, w being the 
Sleight of a unit length of cham Suiee the horizontal com 
ponentsofthe tensions at Pand^must balance each other, we 
have ycosr = (J’4-A3 ’)cos(t+At) Each of these quantities is, 
in fact equal to the horizontal component of the tension at any 
point of the chain, which in turn is equal to the tension at the 
loll cat point That is. m 

(1) 

ctwr 

un the increase in the vertical component of the tension 
fn 111 P to Q must balance the weight vAs Hence 
A(rainr) =■ tcA# 

which by use of (1) becomes A(3«tanr) tv As Hividmg both 
Bules by Az and letting Az approach zero we get the differential 
equation 7^<^tanr or, in z and y coordinates, 

Set T^fw *a a and denote dyfdx by p Then aS «» ./(l-fp*) or 
™ ^ whence aIog{p+VtI+p*}J =* z+c At the 
lowest point on the curve z == 0 and p=i(t Hence c » 0 and 
log{p+V<l+p')}=‘z/a, as p+^l+^) =a e*'® By inverting 
each member of this equation and rationalizing the denominator 
of the left member we get “ *"**“ 

Hence 2p = 

or ^ €-«(“) 

and y =o 

The position of the ongin on tiie axis of symmetry was left 
undetermined The equation of the catenary will be simplest 
if it IS chosen so that Cj = 0 If ^ » 0 then when z = 0, 
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THE CATENARY 
y = a. Wc therefore take the origin a distance a below the 
vertex of the curve. The z-axis is then known as the directrix 
of the catenary. The equation of the curve is now 

IC 

1 / = or i/ = acosh-. 

cs 

Exajifi-E 1. Show that the tension at any point (z, y) on the 
chain is ity, i.e. it is equal to the weight of chain which would 
hang from that point to the directrix. 

EXAStPLE 2. Eind the form of the cable of a suspension bridge 
assuming that the only mass supported is that of the roadbed 
and that the vertical cables are close enough together that the 
curve may be taken as smooth. 

55. Motion In a line or in a plane under knovra forces. 
Suppose a particle of mass m moving in a straight line is acted 
upon by a force F directed along the line. According to New- 
ton’s second law of motion the distance z of the particle from 

drX 

a fixed point of the line satisfies the equation in-~ = F. 

at- 

Again, if a particle is moving in a plane under the influence 
of a force, the force may be resolved into two components F^ 
and Ey parallel to the coordinate axes; the coordinates of the 

particle then satisfy the equations = ^x> ~ ^v 

If the moving body is a complex system of particles rather 
than a single particle, the foregoing equations apply to the 
motion of its centre of gravity. 

Example 1 . Assuming the Newtonian law of gravitation, find 
the speed -with which a body would strike the earth, falling 
from the distance of a fixed star. 

Wo assume also the result, proved by Newton from the law 
of gravitation, that a sphere attracts an outside bodj’’ as if all 
the matter in the sphere were concentralcd at its centre. Then, 
if in is the mass of the falling body and r its distance from the 

earth’s centre, we have = — -. When r = iJ. the earth’s 

ai~ r- 
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radjus then (Pr/d^* = ~g Hence t «= If we' set 

dr/d/ = ~v so that v a the speed of the fallu^ body, the equa 

Don becomes mvdv s= — i— , whence Jme* == ~4-c tVhen ris 

r* r 

the distance to a fixed star t>= 0, hence c is so small as to be m 
considerable Setting c *= 0 we have v == ^(2i./Bir) = ^(2gS^/r) 
TrVTjen the body strikes the earth r=sJi and v = which 
s 6 9 (This IS called the 
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3960} = 


in miles per second : 

velocity of escape If a projectile left the earth at right angles to 
its surfacq with a veloaty greater than this it would not return ) 
Ezaufui 2 A chain 25 feet long passes through a small 
stationary cased mf pudey and hangs at the beginning mth 
13 feet on one side and 12 feet on the other Neglecting finetion, 
find the tune necessary for the chain to slide off and its speed 
jast as it leaves the pulley 

Let the weight of 1 foot of the chain be ic At 
any tune I meaeured from the beginning of the 
26“$ 11 3 motion let a be the length of the longer segment 
of the chain and let T denote the tension in the 
chain at the pulley Then, applying the equa 
tion of motion to each of the two segments of 
the chain, we get 


dP 


= «tv— r, 


(25-a)w^ = T-(i5-B)K^ 

IVom these equations by addition and dividing out to, 
“ (2«-2% 

or if di/dt = V. S5rda *= (2«— 25)yd* 
from which S5~ = («*— 25«)^+c 


t 31 file pulley la not cased m tba ceotnlugal lorce^rill eventually overcome 
Iho pressure on the pulley due to the tetB on and the chain will fly oS before 
the higher cod has reached tbe jndley Hne chain, might, tor purposes ot the 
probU m, bo supposed to pass tbiuu^ an inverted smooth U tube 



MOTION IN A LINE UNDER KNOM'N FORCES S9 
When s = 13, V = 0, hence 0 = — 156?+c or c = loCy. Then 

when 5 = 25. 25™ == I56flr and v 




in feet per second, is tlie speed of the chain as it leaves the 
pulley. The first integral may now be written 

!(*)■=. (,=- 25 ,+ 156 ), 

“ y(^) V(5-— SSsH-lSG) “ 

whence /+Ci = -|log[5— When f = 0, s = 13 

hence Cj = flogi. It follows that 

t = |log[2s— 25+.^/{(2s— 25)=— 1}]. 

\Wien 5 = 25, < = f log{25+.^(624)} = 2-44, which is the num- 
ber of seconds required for the chain to slide off. 

ExAJtPLE 3. If a body falling slowly through a liquid is 
retarded by a resistance proportional to its velocity, find the 
distance in terms of the time assiuning that the body starts 
from rest. 

Exajeple 4. A chain 12 feet long rests in a smooth horizontal 
tube with 1 foot of its length hanging out of the end. Eind the 
time it takes to slide out and its speed as it leaves the tube. 


56. Harmonic motion. Resonance. If a particle moves in 
a straight line with an acceleration which is proportional to the 
distance of the particle from a fixed point of its path and wliich 
is directed towards this fixed point, then the particle is said to 
execute simple harmonic motion. If s is the distance of the 
particle from the fixed point, the differential equation of the 
d-s 

motion is — -f- a=a = 0 (see problem 27 et scq., p. 61 ). The solu- 
tion of this equation is s = a cos(a<— 0), where a and 6 are 
arbitrary constants which in a specific problem ma}' be deter- 
mined from the initial conditions. When t increases by 2tt/(x 
both 5 and dsjdf return to their former values, hence the motion 
is oscillatory with a period of 2n-/o:. As the value of s oscillates 
between a and —a, 2a is called the amplitude. This motion is 
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approxunately realized bp a vibrabmg tuning fork or bp a simple 
pendulum swinging through a small angle 
If the particle in its hamonic motion encounters a resistance 
which 13 proportional to the Telocitp the acceleration 13 de- 
creased bp an amount which map be denoted bp 2^^ The 
resulting equation la ^ auuliarp algc 

braic equat on is =» 0 vhicn has roots 



If ^ > a these roots are real and the differential equation has 
the solution « « “^'+e,e while if ^ < a 

the roots are imagmaip and the solution of the differential 
equation is 4 = where a and 9 are arbl 

traip constants In the former case each term is of the form 
ee which is a pc^itire decreasing function whose value 
approaches zero as t increases Bence the displacement $ 
gradoallp approaches zero without ever passing through that 
value Such a motion would be realized bp a pendulum bus 
pended in some thick bguid in which the resistance is so great 
that the pendulum sinks to rest without a single oscillation On 
the other hand if the resistance is small so that ^ < a the 
second solution applies This represents an oscDlafion with 
penod 2w/^(af*— ^)t and with conatantlp decreasing amplitude 
t In thn case tbs t«iin penod i* to to undentood an tbe atervn) between 
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It is referred to as damped harmonic motion. This is realized 
approximately by a pendulum swinging in the air in which case 
^ is small and the oscillations persist for a long time, or by 
a tuning-fork which is struck and then immersed in water, in 
which case the vibrations are quickly damped out. 

Consider next the case in which the vibrating system is acted 
upon, in addition to the force of restitution and the damping, 
by an external force which interferes with the free vibrations. 
The most interesting case is that in wliich the applied force is 
periodic. As an example one may think of a vibrating timing- 
fork near one prong of which is placed an electromagnet with 
an alternating current passing through the coil. In attempting 
to represent such an applied force it is natural to use a trigono- 
metric function. We shall therefore consider an applied force 
k cos yt, k and y being constants. The equation of motion is then 

^ + 2^^+“'® — A: cosy/. 

The complementary function has already been found (we shall 
take ^ < a). A particular integral is 


I)-+2§D+ 


1 , , i)2-j-a2-2j5Z) 7 


= k{D"-+a^-2pD)- 


cosy/ 

(<^-Y‘r+ApY 


j (or — ■)r)cosy/-i-2j3ysiny/ 

(«“-r)-+W 

fccos(y/— 

V{(or— y=)2-f4^^y2} 

where tan ^ . The complete solution is then 

a-— 7^ 


«■= ae-^cosy{c^-p^)l-e}i.- 




COs(y/— ^). 


The complementary function gives the free vibrations of the 


system which in time are damped out, and the particular integral 
gives the forced vibrations. The amplitude of the forced vibra- 
tions is 2ky{{<x-—y^)--\-4PY}- If y be thought of as a variable, 
the other quantities being fixed, it is easily shown that the 
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ampbtade is greatest ■when y *= ^(a^— 2^) If ^ is small this 
means that the penod of the applied force S-r/y is nearly equal 
to that of the free vibrations, In this case the 

foreed vibrations which persist afUi the free vibrations are 
damped out are of large ampbtude 
On the other hand if the damping la small and the two 
periods mentioned nearly equal we may get a better idea of 
the motion in another way For the case ^ = 0 and a = y the 
equation haa a slightly different solation (see ex 29 p 61) 
The particular integral is ^ tcoaaI=» > which shows 

that the smphtude of the forced vibrations constantly increases 
■with the tune Such motion may be approximately reahzed in 
vanoos ways Thus for example if soldiers march in step 
across a bridge and if the penod of their step catresponds with 
one of the natural periods of ■nbiation of the bridge the vibra 
tiOQS may become dangerously large For this reason soldiers 
break step while crossing a bridge 
When the penod of the applied force and that of the free 
nbrations are nearly in tnne ao that the vibrations ultimately 
reinforce each other the reeulung phenomenon la known aa 
resonance 


EIAllPLTA OV CHAPTER IX 

1 Detersuae the carve wluch w everywhere concave to the s-axis 
and tor which the ndiua of curvature is twice the length of the normal 

2 IVeteinunetbeairvewhichisrveTywherecoirTextolhex-BXisazid 
for which the radius of c.-urvature is eq^ to the length of the norznah 

3 l>etenmne the curve f<w which the curratore is at any point equal 
to the slope of the normaL 

4 Determine the curve which at the point (0 8) is tangent to the 
line y >3 2 and for which the radius of curvature at any point is equal 
to one-half the square of the ordmato of that point 

5 Detemune the curve which toochee the z-axis at the ongm and 
for which the length of arc from • 6xed point A to a variable pomt P 
la proportional to the slope at P 

d A cable supporta a sreight after the tnumer of a eogpensioa bridge 
Assuming tl»t the density of the mses supported is proportional to the 
horizontal distance from the lowest pomt of the table the curve 
assuroed by the cable 
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7. Suppose that a set of uniform rods are threaded, like needles, on 
a string and that each rod hangs vertically- If the string is suspended 
from its ends and the rods are cut off so that their lower ends are in 
a horizontal lino, find the curve assumed by the string. 

8. Solve the following equation of motion, and find the period of the 
free vibrations, and the period and amplitude of the forced vibrations: 

2S+2|+9s=8cos2t. 

9. A chain 15 feet long hangs through a small frictionless pulley with 
8 feet on one side, of which the end just reaches the floor. In what time 
will the chain slide off the pulley ? 

10. If a body falling in air is retarded by a force which v'aries ns the 
square of the velocity, find the distance and the velocity for a body 
falling from rest for I seconds. 

11. The force of gravity at any point within the earth is that due to 
the sphere of matter nearer to the centre than the point. If a shaft 
were sunk to the centre of the earth, find the time necessary for a body 
falling in it to reach the centre and the velocity with which it would 
arrive. 

12. Set up and solve the equation of motion for a particle which 
moves in a straight line towards a centre of force if the force varies 
inversely as the cube of the distance. Assume that the particle starts 
from rest at a distance a. 
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ORDINARY DIFFERENTIAL EQUATIONS IN 
MORE THAN TWO VARIABLES 
Ir one independent and tm dependent ranablea are involved 
the geometncal interpretation of the differential equations re 
quires a space of three duaensions Moreover it is necessary to 
consider the case of a aingle differential equation and that of 
tuosimultaueous equations A preliminary investigation of the 
geometrical meaning of the equations mil sLov what » to be 
expected m the wav of solutions 


57 Geometric meaning of a pair of simultaneous equa- 
tions ol the first order If a curve is given by the two equa 
lions 1 / r] 8 0 and p z) » 0, then a tangent hue at 
(»« y# *«) on the enrve haa equations 


The direction coemes of the tangent are proportional to 

1^ 

jtV ^ 

ey'i 

Differentiation of the equations of tbc curve gives 

_<ir+-Js, + g*« 


cu cu 

.?tl cu 

ey ct 

r "i 

ec gt 

W fe 

ey « 



— rfzH — dy+— rfr = 
cz eg & 

from which it appears that tbe direction cosines of the tangent 
are proportional to dr dg and d. 

Take now a eet of two equations of the first order 
Pldz+Q^dy-e•Jl^d^ = 0 1 
P.dx-i-QfHg~rSfdz 0 I 


( 2 ) 
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GEOMETRIC MEANING 

where coefficients are continuous, single-valued functions of a-, 
y, and s. These equations determine unique values at every 
point for the ratios dx\dy: dz. This may be interpreted as 
defining at every point a definite direction. We may therefore 
imagine each point of space provided with a line element, and 
it is reasonable to expect that these line elements will determine 
a set of curves of which one passes through every point. Each 
of these ciurves will in general cut a given plane and conse- 
quently they are as numerous as the points of a plane. In other 
words, the curves constitute a two-parameter family and require 
two arbitrary constants in their equations. 


58. The equations ^ ^ Special devices for solu- 

P Q R 

tion. When equations (2), art. 57, are solved for the ratios 
dx ; dy ; dz we obtain 

dx dy dz 

QiRz-QzRi. ~ Rj^Pz-RzPi ~ PiQ^-PzQi 
This form, which we shall abbreviate into 

dx _dy dz 

P~~Q~B’ 


( 1 ) 


is more convenient for the study ofihe equations. We saw that 
the general solution consists of two equations involving two 
arbitrary constants. We shall sometimes for convenience speak 
of these equations separately as solutions. Some devices will 
now be considered which are sometimes useful in obtaining these 
solutions, (i) If one of the equations in (1) is or can be made 
firee from one of the variables, then it may be possible to solve 
this equation. Sometimes it is possible in this way to obtain 
two integral equations, which constitute the general solution. 
Or, if one solution only is obtained, we may make use of this 
together with the given equations to get a second differential 
equation involving two variables, (ii) From (1) we obtain 

^_dy dz Xdx-{-iJ.dy-{-vdz 

P~'Q'~R AP-f/x^-fvi? 

provided XP -\-fiQ-\-vR ^ 0, A, y., and v being any multipliers. 



W EQUATIOMS IN MORE THAN T^VO VARIABLES 
Moreorer li XF+ftQ-f-vJi == 0 it is easjy shown that 
Xrfa: -f-yit = 0 

It IS BOznetimes possjh)c to find multipliers A ^ f which mahe 
XP+nQ+vR = 0 and at the same tune make 
Ad*+/tdy+»'<fes= 0 

an integrable equation Agam if the fourth ratio in (2) can by 
proper choice of X p v be made to take the {brin du/u it may 
be possible to obtain an integrable equation by making use of 
the ranable t* 


a y x+z 

The first equation gives the solution * => log^ or y = ce"* 

Also ^ SR from which y = (itt-fs+o) The integral 

curres a?» the intersections of the exponenhal cylinders y a» «* <* 
with the planes y « c,(ar+2+o) 

Eiaiiplis -iJ . 

TM— ny nx—ls ly— »w 

By using successively the multiphers I m n end x y s we find 
Idt +mdy •\-ndz ^ Q and rdx+ydy+sds » 0 These give 
the general solation Iz-t-my+nr = a **+y'+** ■* 6 The 
curves represented are theuitersectiooa of a set of parallel planes 
with a set of coi centne spheres in other words they are circles 
whose planes ate parallel end whose centres he on a straight 
line perpendicular to these planes 

Exuoiss 

ity y* s* 


EXAtlFLB 4 


dx _ dy dz 
J/+6 x+a *+c 


Examples 

a* y ** 


59 Simultaneous linear equations with constant coeffi- 
cients It 13 generally an aid towards the solution of a pair of 
snnultaneous equations if one of the dependent vanablea can 
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SDIULTANEOUS LtSTEAR EQUATIONS 
be eliminated. This is possible in a -nide variety of cases, the 
resulting equation being frequently of higher order than either 
of the given equations. In the case of linear equations in which 
the coefficients of the dependent variables and their derivatives 
are constants, the elimination can be effected by an algebraic 
manipulation of the operator D, An example will show the 
method. 


Exasiple 1. 







= 3(21-6-*). 


To eliminate z and dzjdx algebraically would require three 
equations instead of two. If we differentiate each of the given 
equations we introduce a second derivative d-zjdx^, but we have 
then foin equations from which to eliminate three unknowns, 
I s, dzjdx, and dhjdx?. The two additional equations are 


dx?' 


dx^ dx dx 


3(2a:+e-*), 


dT? dx? 


= a, 


The elimination is now effected by adding the first, third, and 
fourth of these equations and subtracting the second. This gives 

30+3y = 6e-*+3a;2+6. 

We note now that the foregoing operations may be condensed 
as follows. Write the original equations in symbolic form, 

(Z?+2)y+(J5— l)s = 3x=-3e-* 
{2D—l)y—{D+l)z = 6*- 3c-*. 

Operate on each term of the first equation by 1 and on each 
term of the second by D— 1 and add the results. This gives 

3(2)Hl)y = 6e-*+3r'^6 

which is the equation obtained before. The general solution of 

C12 
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this equation is easily found to be 

y r^+i'-J-Cieosi+ejBina: 

It remama to find z For this we have a choice of two 
methods 

(1) From the given equations ebmuiate I>s (by addition in 
this case) In the result subaUtnte the value }ust found for y 
and solve for z Thus 

(3Z<+l)y— 3;t*+6j;-6e'^. 

•Rhenco 2 = 2 e~*— ^^l^^sinr 

(2) Use the aame method to find z as was used for y This 
gives a=» 2«~*— **+c,co8a;+C4ein* It is to be noted, how 
ever, that the constants e, c, «, are not independent The 
relations among them may be found by substituting the values 
found for y and s m cither of the given equations and equating 
coefRcients of like terms in the resulting identity Thus, sub 
stitutuig in the first equation vc get 
—<'*+ar—Cjami+c,c08»—2<-'—2r—e,8m*+ 

4-C4C083r-h2e'*+2ir*+2(ijCos*+2C|Bin*— 

— 2e-*+**— C|COJar— Cjsmar =» 3z**-3s**, 
whence, on equating to zero coefficients of cosz and of sinz, 

2Ci+c,— Cj+c, -= 0. -~Ct+2ct—e^—e^ •» o, 

or c, =. gi-3ct 

EsaSifle 2 z = 0, y=0 

ax dx 

ExAMTLB 3 ^+^—2? 2cos*— Vsinx, 

dx dx 

^ + 2y= dcosz— Ssmz 
dx dz ' 


Exsmtu 4 
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60. The equation Pda:'-]-Qdy-{-lidz — 0. Condition for 
integrability. From the equation u{x,y,z) = a, involving an 
arbitrary constant a, a differential equation 

^dx+^dy+^dz = 0 ' ( 1 ) 

6x dy ' dz 

may be foimed. This equation (1) is exact and has « — a as 
its solution. In order that the equation 


Pdx-{-Qdy-\-Rdz = a (2) 

be exact there must exist a function u such that dujdx ~ P, 
dujdy — Q, dujdz = R. From these equations result the fol- 
lowing necessary conditions for exactness: dPjdy — BQ/Bx, 
dQfSz = BR/By, BRjBx = BP/Bz. It is sometimes possible to 
render (2) exact by means of an integrating factor. In such 
case (2) is said to be integrable. We proceed now to determine 
a condition on P, Q, Rin order that (2) be integrable. Suppose 
X{x, y, z) is an integrating factor so that 

XP dx -|-AQ dy -j-Ai? dz = 0 

is exact. The above conditions for exactness now give ■ 

BXR BXP 

By Bx ’ 8z By ' Bx Bz ' 


From the first of these equations A— -j-P— = X^ + Q— or 

^ By^ By Bx Bx 



Multiply these three equations in order by R, P, and Q and add 
the results. A is then eliminated and we have the following con- 
dition which is necessary for integrability of (2): 



= 0. 



lOO ZQ^;ATIO^S IS MORE TOA?. TWO VABUBLES 
Tlua condition ts in fact snSorat os well aa necessaij' bottlie 
proof will be omitted here 

61 Geometnclnterprelatlt>Qofi*dr+Qdy+ild =0 We 
hare seen that geometric meaning may be giren (oeqnationaio 
three Tsnables by interpreting dz dy dz as direction nnmbeis 
of a line element in apace Aa a single equation connecting 
dx dt) dz a not sufficient to detennme their ratios it follows 
that at anr point an infinite number of line elements will 
“‘“fy Pd*+Qdy4-«dz=0 {!) 

The geometnc conf^uration suggested by (1) u therefore a set 
of enrres in apace of which an infinite number pass through 
ereirpomt Moreorer the equation states that at any specified 
pomt a curre whcee tangent has direction nnmbert 6x dy dz 
IS perpendicular to a cuire whose tangent has direction nntabers 
P Q R K curve of the tatter sort u that defined by the 

^ ^ ^ 

P “ Q “ P ^ ' 

It follows then that through any point tbe enrres which satisfy 
0) a» at right angles to the amgle ciirre which satisfies (*} 
Agtun if (1) u utegrable we taw that ita soluiioa is of the 
fonn «(7 y s) o Inthiscasetbeequationmaybemterpreted 
as representing either surfaces or carves the latter constitutmg 
the family of all curves whi(^ be on tbe snr&ces It follows 
that m the mtegnble case of (I) tbe surfaces which it representa 
ace the orthogonal trajectories of the carves denoted by (^) 

If (1) u not integrable there is no difference in its lotcrpreta 
tion by means of corres bat tbe interpretation by means of 
surfaces bieahs down Tbe family of curves ("} haa in tbi3 case 
no ortbogonal sor&ce tnjectones 
£xa>rPLB 1 Interpret geometneally the equation 
zdz -hjrdy-f-zds k 0 

and consider its relation to dxjx ^ dy}y « d^jz 

2 The tame problem for zdy+ydx4-zd2 0 m 
relation to dr/y = dy/r = dz^ 
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GEOMETRIC INTERPRETATION 


C?2/ 

Example 3. The same for h— -1 — = 0 in relation to 

X y z 

xdx = ydy — zdz. 


62. Solution of Pdx+Qdy +Rdz — 0, when integrable. 
Por some integrable equations of this type the solution is evident 
or may he readily found by inspection. In particular this is the 
case when the equation is exact or when a proper grouping of 
the terms makes evident an integrating factor. A method of 
general application, however, is the following. Let one of the 
variables, say z, be considered constant and solve the resulting 
equation in x and y. In this solution replace the arbitrary con- 
stant by a function of z and attempt to determine this function 
60 that the given equation is satisfied. Example 2 below illus- 
trates the operation. 

Another method may sometimes be used to advantage in case 
the equation is integrable and the coefficients homogeneous of 
the same degree. This method, analogous to that used for the 
homogeneous equation in two variables, consists in the substitu- 
tion of 2w for X and zv for y. In the resulting equation it may 
happen that z disappears entirely, in which case an equation in 
the two variables u and v is left. Apart from this exceptional 
case, however, the substitution leads to an equation in which 
the variable z is separate from u and v. Since this equation is 
integrable it is easily shown to be exact. 

Example!. 2xzdx—2yz-dy—x^dz = 0. 

By grouping together the first and third terms we observe that 
1/2- is an integrating factor. When it is used the equation takes 

the form^ ^ — ^ 2ydy = 0 and the solution is ^ = c. 

2 - z 

Example 2. 2xyz dx — z(x^-f-z) dy —y{z--\-y) dz = 0. 

Take z constant and hence dz = 0. In the resulting equation 


the variables x and y are separable, giving 

X-4-Z 

whence ^ = c. Now replace c by /(z), 

2 / 


2xdx 

x--\-z 


y 


0 . 



( 1 ) 
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Differentiate this and compare the result ■with thegiven equation 
y[Zxdx (:^+z)<?y _ f^z]d2 

y* 

or 2jy<te — (**+*)dy+{y— = 0 

When thia la multiplied tLnmgh by a the resulting equation 
agrees term for term with the gtren equation provided 
yz-y»-/(r)= -y(**+y) 
or yzf{z) = a+**+y = lif(2)+y 

d{{z) dt 

Hence {1) is a solution provided 2/'<z) =/(2)+ 1 or = — 
from which /(2)+ 1 *= 0* where a is an arbitrary constant Sub 
aUtutu)g/(z) » as~l in (1) we obtam the solution 

I or **+2d'y'»oy2 
Eiuluple 3 By settmg z « zu tf^ tv soke 

(y*+js)dr-{-(zr+z*)dy+^— zy)ife« 0 
EiaaotE 4 (Ziy— 3-)iz+(z»+2i‘)dy+{4j«— 3*)d2 ■* 0 
ExattPLE S s*y<he4-(y*2~**»)dy— y*ii3 »* 0 
EXUOLEtl 

(2«— y3)ie+(2yz— Z2)dy— (z*-*zy+y*)ds 0 

63 Solutions of Pdx+Qdy+Bdz => 0 when dOt iate> 
grable In this case the cuires given by dr/P * dy/Q ~ dz/S 
have no orthogonal eurfaco trajectones However there is no 
lack of curves orthogonal to those of this famOy In fact any 
arbitrarily given surface w completely covered by curves which 
satisfy the given equation and which are orthogonal at every 
intersection to the curves dxfP = dyjQ »= dzjB For let 
/(z y z) = a denote an arbitrary surface Any curve on it 

satisfies the equation ^<iz+— dy+^ds=aO But this forms 
dx 0tf ^ 

with the given equation a pair of aimultaneous equations which 
denote a set of curves of which one goes through each point of 
space Hence the surface / = a ts cwnpletely covered by these 
curves 
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SOLUTIONS WHEN NOT INTEGRABLB 

Example L Find the curves which satisfy the differential 
equation {x-\-z)dx-\-(y+z)dy -)r(z—x—y)dz = 0 and lie on the 
sphere x--\-y^+z^ = r®. The curves must satisfy also the dif- 
ferential equation xdx-{-ydy-{-zdz = 0. Subtracting this from 
the given equation we get zdx +zdy —{x+y)dz = 0 or 

dx-\-dy dz ^ 

x+y 

from wliich x-\-y = kz where k is an arbitrary constant. The 
required curves are the intersections with the sphere of this 
family of planes. 

ExAirPLE 2. Find the curves which satisfy 

{l-\-2a)xdx +y{l—x)dy —zdz = 0 
and lie on the hyperboloid x^-\-y^~z^ = c*. 

Example 3. Find the curves which satisfy 

ydx-{-{y+zx)dy~[z-\-xy)dz = 0 
and lie on the surface x-^y ~ zx. 


EXAJIPLES our CHAPTER X 

1 n dx dy dz, 

■ I’+z* ” ~ 2z' ■ 1 /S— z’’ ~ z^—x^ ~ a:>— j/“' 

3 ^ ^ dz ^ dx _ dy dz 

y ® y')' ’ X ~2y'~ cos[x'^-{-zy 

^ ^ — dy _ dz 
' xy~ y^~ xyz—ix^' 

= 2g + 5g+6z = ®. 

•jdy.^dz 


8— 4-^'.^- ft dx dy 

9 . {y—z)dx ,>r{z—x)dy ->f-{pt—y)dz = 0. 
10. ydy+zdz = V{l-y’- 2 =)d*. 


11. (xdxY,^{ydyy—{zdz)",^2xydxdy = 0. 

12. yz'(x^~yz)dx JrxH{y'^~xz)dy ,^xy>[z^-xy)dz = 0. 

13. 2«(la:+2!/zrfy-(a;S+t,>_jS)rf3 = q. 
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11 yz{y+»)dr+«(*+*)rfy+*s{i+y)* *= 0 

15 lr(y»— »)££ j:+{s*+*— y)ii = 0 

18 (i^— a^«— + — 0 

17 Puidtiiee«tuatioBaofsciiive(luoiigbUiepo0t(I>I I)ortbogonAl 
to the famJy of surfacw (*•— jr*— 2rjfrfy-f ted? <= 0 

18 Find the curves on Uw eOipMid * which eatiafy 

19 Fnjd the cttTPW oa the cylinder* v*+*’ “ r* (f eilntnuy) whxh 

Mtisfy the equAtioa dx — y(tz— l)dy 4-*dz 0 

SO lo order thet the equetioo Pdz +Qdy+77d:-f5dl 0 be m 
tegrabla it most be intcsr>bte when nny one veneble ts eonitant From 
this feet state four eondittons which are necessary for lotegrabilit} and 
show that any one of these may be obtained from the other three The 
fonowingtsagenetal method of mlr^ntisg such rquaticins Take sand 
I constant and aoira the resuhing equation. RepUce the arbitrary con 
stant in the eolation by /(s t). Differentiate the new equation in the 
fonjr WMblee and compare with the ongmal diffeTetitial aquation to 
fiad/<t t) 

*1 Urd-a+Dfc -l-U-f »+*)dy +Il4-a4->>ila +U-l-y+»)dt - 0 

2S (y*+ftrs*)da +(Sry— 6y%*— <*)dy — (<yH— 4**t— 8f*Jd* + 

+(lr*a->'3yr>-f &r)df • 0 
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PARTIAL DIFPEEENTIAL EQUATIONS OF-* 

THE FIRST ORDER 

Partial differential equations are those in which the derivatives 
involved are all partial derivatives. This implies that at lepst 
two variables are independent. This is the case in many physical 
problems in which, for example, the time and one or more 
coordinates of an arbitrary point vary independently. As a con- 
sequence partial differential equations play a dominant role in 
mathematical physics. The subject is a large one, and only a 
brief introduction is in order in an elementary book. The present 
chapter is concerned with equations in which only first partial 
derivatives appear. In order to furnish an idea of the kinds of 
solutions to be e.xpectcd, we show fn’st how such equations may 
be derived from their primitives. 


64. Formation of partial differential equations by the 
elimination of arbitrary constants. Let (f>{x,y,z,a,b) = 0 
be a known relation involving two arbitrary constants a and b 
and two independent variables which may be taken as x and 
y. By partial differentiation we obtain 

^ I = 0 

• dx'dz dx 


d(f> d(f> dz 
By dz dy 


From these two and the given equation the constants a and 6 
may be eliminated, the result being an equation of the first 

ordor,/|a;,y, 2 ,^,^J = Oovf(x,y,z,p,q) = 0, where we employ 


the notation — 
dx 


oz 


Vi, 

In the following examples partial differential equations are 
to be formed corresponding to the given equations, a and b being 
arbitrary constants to be eliminated. ® 

’*13 
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ExAsms 1 (x— o)*+(y— J)*+s* = f* 

By partial differecttatioo mth r^ard to x and y w get 

X— o+r— = 0, y— »= 0, 



or 2*(p*-H®+l) = r* 

Et-ouxe S (x— o)*-r(y-“0)*+{*— <*)* = ^ 
ExAsmj: 3 x = «x+fcy+«3 
Example 4 2 c= (X'(-a}(y4'&) 


65 Forxoatioa of partial differentia] etjnatlDfls by tbe 
ellnunation of arbitrary functions. This metbod has no 
analf^e m tie case of ordinary differential eijnationa and is 
aeoordiagly of pecoliar ngni&caace ta tie loterjwtsboo of 
pastsal ^eRQtial equtioas. The equation * «» *+/(y) yield# 
by diffenotiahoD cz/cx sb 1 so matte; «hat fonetion /(y) u 
Siadlariy, froms = /(z+y)efededneef»/&t » Ps;cy These are 
particolar examples of the foOomng general ease Snppoee « 
and r are t<ro known fimction# of x, y and 2 and that ti and r 
are connected by an aibitraiy fonetic^ relabon which may be 
symbolized in the form dC**!**) ** 0 or c «= /’(«) We aaemne 
only that the fonctioas ^ and t as well as u and r possess first 
partial denratires with respect to all their vanables Then 
from 4l«, c) = 0 we get 




ett\cV ^ct)^ct\c3^^czl 


?e , ev 
=;+*5: 
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which takes the form 
I3u8v dudv\ 

\& 8z dx)^ ~ dy 8x dx dy 

or Pp+Qq = -fi where P, Q, and R, the coefficients in the pre- 
ceding equation, are known functions of x, y, and z. 

The results of the present article suggest that in the general 
solution of partial differential equations arbitrary functions play 
the role which is played by arbitrary constants in the case of 
ordinary equations. This will be borne out in the subsequent 
discussion. 

In the following examples partial differential equations are 
to be formed by eliminating the arbitrary function /. 

Example 1. z =f{xy). 

Example 2. z = f{ax+by). 

Exasiple3. z=f{x-+y-). 

Example 4. = 0. 


66, The linear equation of the first order. Lagrange’s 
solution. By eliminating the arbitrary function ^ from the 
equation ^{u, a) = 0 we obtained an equation of the form 

This is called the linear equation of the first order.f A relation 
w(a;,y,r) = a will be a solution of (1) provided the derivatives 
8zjdz and dzjdy obtained from it reduce (1) to an identity. From 
M — a we have 

& _ Bu jdu 8z 8u jdu ' 

8x Bx/ dz’ By dyf dz’ 

Hence (1) is satisfied if 

( 2 ) 

1 Tho distinction should bo noted between tliis definition and that for 
ordinary differential equations. In this case tho equation is linear only in tho 
partial derivatives of tho dopondcnt variable, not necessarily in the dependent 
variable itself. 
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that 13, if the function « satisfies (2) Again, bom v = a we 


la order that (2) and (3) be true simultaneously we must have 
* * 

e Q B ^ ' 

where the differcnOals refer to the surface tt =« fl Conversely, 
if « =» a IS an integral of (4) then from (4) and (3), (2) fallows 
whence by division by Iv/Bz (I) foUons Hence u *= a will be 
a solution of (I) if and only if it is an integral of the set (4) 
Suppose now that u » a and t> = 6 are two independent 
integrals of (4) IThen ^|«, t) = 0, where ^ is an arbitrary func- 
tion 13 also a solution of (1) For by the ehmioatton of ^ we 
obtained (art 65) the equation 

t duSv m 

ft By ByctjBz \BxH Bs Bx) By By Bx BxBy 
Also, EUiee and tr «« 6 are integrals of (4), we have 

from (t) ft,, a,. ft.. 

4>|+^S=<>' 


SuSv &uSv Bu9v PuBv BuBv Budv 

By Bz Bz By BzBs Bx Bz BxBy By Bx 


Hence (5) is equivalent to (1) 

The foregoing solution of the linear equation u due to 
Lagrange The ordinary equations (4), upon which the soluDon 
depends, are called the subsidiary equations The curves which 
satisfy (4) and furnish the hey to the solution of (1) are called 
the chaTo/tlenatic euri'es or cAoracteruties of (1) For a given 
function ji the solution ^{a, v) 0 (or v F(u)) represents a 
surface which is completely covered by characteristics of which 
one passes through each point <«> the surface With ^ arbitrary 
the solution consists of all surfaces which can be built up from 
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LINEAR EQUATION OF FIRST ORDER 


cliaracteristics. Through an arbitrary curve C wliich, is not 
itself a characteristic there passes a surface which is the locus 
of all the characteristics 'which intersect C. The determination 
of such a surface through an arbitrary curve C is a particular 
case of a more general problem first considered by Cauchy. 

The foregoing method of solution applies as well when the 
number of independent variables is greater than two. Thus if 
Xi, Xo,..., x„ are independent, the form of the line{}r equation is 


The subsidiary equations are 


•+-P/I 


dx„ 


= E. 


dxi 


dx^ _ dxs _ _ 


P,. 


dz 

E' 


(6) 


and if Vi{x,y,z) = a (i ~ 1, 2,...,n) are n independent integrals 
of this set, the general solution of (6) is = 0. 


Example 1. a-+b~ = 1. 

dx dy 

The subsidiary equations are dxja = dyjb = dz. Two inde- ‘ 
pendent integrals of this set are x—az = and y—bz — c^. 
Tims the characteristics are a family of lines parallel to the line 
x — az, y — bz. The solution of the given equation may be 
written y—bz = F{x—az) or <l>(x—az,y—bz) = 0 arid denotes 
all surfaces which may be built up from these parallel lines, 
that is, all cylindrical surfaces with elements parallel to the line 
X = az, y = bz. 

Example 2. {7nz—ny)^+{nx-h)^ = ly-mx. 

OX dy 

The subsidiary' equations are 

dx __ dy dz 

Ttiz—ny nx—lz ly—mx’ 

Two independent integrals of tliis set were found (art. 68) to be 
Ix-\-Tny-\-nz = a, 

= b. 

Thus the characteristics are the intersections of a set of con- 
centric spheres with a set of parallel planes. Tliey are therefore 
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circles vhicli lie m parsUcI planes, at right angles to the Une 
xjt = y/m = zjn, and hare their centres on this line The 
general solation = 0 consists of all 

snriaces which can be boilt op from these circles, that is, 
all sorfaces of revolntioo ahoot the hoe z{l =: y/m « z/n. 

Example 3 (z—o)^+(y-‘b)^^z—e Solre and inter 

pret geometrically 

Sz 

Exampus 4 SolTe and interpret geometn 

cally 

67. Standard forms for which solutions are easily found 
Under this head we consider four types of equations of the first 
order and obtain solntiona each lovolTing two arbitrary con 
stanta 

I Epatvma <ntA only p and q pruent Sneh an eqnabon, 
F(p,q) «a 0, IS obnonsly satisfied by s ^ or-f 6y+e pronded 
a and b are chosen to satisfy the relation Fip, b)^0 Thus, for 
example, j»*+3+ 1 « 0 has the solation * * or— (a*+l)y4-c 
An eqoation may eomebmes be redoced to the form E(p, y) 0 
t^meansofatranafonaabonofraruble^ Tbusmflzp.g) « 0, 
if we act X » logz, we hare zp » &/3X =* P, say, and the 
eqaation beoocies/(P,gJ — 0 

n Equatwnamihordyp.q andzpramt If in i'ls.p, ff) == 0 
we suppose z a foscbon of the single variable u = le-f-ay where 

^ dz Sv tSz dz Stt 

a 13 an arbitrary constant, we have — = ^ ana ?- = t-t- 

■' e* dwhi ^ dnby 

or p = g = a— Hence the ongioal eqoaboti is satisfied 
^ ^ general solution of this ordinary 
eqnabon is a solution of the given equation with two arbitrary 
constants As an examp^, the mdi^ry equabon correspond 
ing top* -f-g* = 42 IS (l-ha*){dz/dw}*s=! 42 which gives the solubon 
0+o*)z == or (1-f fl*)2 »= (x-j-ay+6)* 

in Esuations o/ the/om/(*,p) = ^(y,g} In this equabon 
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the variables are said to be separated. If each side of the equa- 
tion is set equal to an arbitrary constant a and the resulting 
equations solved for p and q, the %'alue of z is then found by 
integrating dz = pdx-]-qdy. Applying this method to the 
equation px- — g-d-y we have p = a/a;®, q = 'whence 

dz = ^dx +^(fl—y)dy and z = — |(a— y)l-f 6. 

X“ X 

IV. Equations of the form z = px-\-qy-\-f(p, q)- This equation 
has a close analogy to Clairaut’s equation (art. 17) which is 
solved by replacing the derivative by an arbitrary constant. 
This suggests as a trial solution for the present case 
z = aa:-f 6y-|-/{a, b). 

From this equation we get by differentiation p — a, q = b, and 
hence the differential equation is satisfied. 

Example 1. p--\-q^ — 4. 

Example 2, pz+q = z®. 

Example 3. z = px-^qy-{-4p—q. 

Example 4. yp — xq. 

Example 5. pq — p+q. 

Example 6. z® = pq. 

Example 7. ^J{l—x^)p = xq{l+y^). 

68. Complete, general, and singular integrals .f We have 
seen that a partial differential equation of the first order may 
be formed from a primitive by the elimination of either two 
arbitrary constants or an arbitrary function. In the converse 
process of solution we have seen how to obtain for the linear 
equation a primitive involving an arbitrary function and for 
the equations of art. 67 primitives each involving two arbitrary’'' 
constants. The latter type of solution, that involving two 
arbitrary constants, has been named by Lagrange a complete 
integral, the former involving an arbitrary function, the general 
integral. It remains to point out that the same equation yields 
both types of integral and that from an integral of either type 

t The discussion of this nrtiolo applies to partial diSorential equations irith 
itro independent r-arinbles. 
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one of the other may be Smnd In the first place it is clear 
that from the general integral a complete integral may be found 
by replacing the arbitrary fonctuHi by a specific function in 
volving two arbitrary constants Infect any number of com 
plete mtegrala may be foond in this way and any complete 
integral appears as a two parameter &tmly of surfaces contained 
in the general integral The name complete integral appears 
from this pomt of view to be a munomer It gets a measure 
of justification however from the fact that the general integral 
may also be obtained from a complete integral This is done in 
the following manner Let jr z a 6) = 0 be a complete 
integral with the arbitrary constants a and b This denotes a 
two parameter family of surfaces Ifweact b = q5(a) where is 
a known function we pick out from the complete integral a one 
parameter ■family of surfaces This one parametei family will 
in general possess a surface envelope whose equation may be 
found by the usual process At any point of this envelope the 
values of X y z p q determine the equation of the tangent 
plane and hence agree with the corresponding values for one of 
the surfaces o t Hence the envelope is it«elf a solution 
of the ditferentialequatioQ If now ^(a) be taken as an aibitcacy 
function we reach the conception of an infinite Kt of envelopes 
of one patatneter famihea of the eotfaces f »= 0 These eon 
Btitute the general integral If it is possible to find for it a single 
explicit equation the latter will contain an arbitrary function 
The elimination process involved in finding an em elope cannot 
generally be earned out when ^(a) is arbitrary Wo may how 
ever choose particular functions ^ in a great variety of ways 
and obtain particular surfaces beloi^ng to the general integral 

A third type of integral exists in case all the surfaces of the 
general integral (which includes any complete uilegtal) possess 
an envelope This envelope may be found from a complete 
integral by elunmatuig a and b from the three equations F 0 

t Since t y t p q deteniune a tangaot pldne we may think of such a a«t 
of values M a plant ekmtnt which we nay wuaUze as a small square or 
circular port oa of a pUoe with the po nt (z y z) at ts centre Plane elements 
are useful lA the geometric interjaetabmi of part al differenUal equal one os 
me Ime (^cnrienlii m ot or&neiy fidSe enliol equal one 
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dFjda = 0, dFjcb = 0. In this process it may happen that loci 
other than envelopes appear. The final test for an envelope is 
that it satisfy the differential equation. 

The three types of integralf are illustrated in the following 
examples: 

Example 1. z = 

Tills equation belongs to type IV, art. 67, and has for a com- 
plete integral 

The envelope of this two-parameter family of planes is found 
b}' eliminating a and b from this equation and 

This elimination, which is left to the student, gives 

The singular integral is therefore a sphere. The general integral 
consists of all surfaces which are enveloped by the planes 
2 = ox-f ^(a)y— A:4l-}-a2-}-{^(c)}-}, 
where ^(o) is an arbitrary function. It is clear that any one- 
parameter family of tangent planes will envelope a developable 
surface tangent to the sphere and that a surface of this type 
exists which touches the sphere along any curve on it. All such 
developable surfaces constitute the general integral which in- 
cludes, thei-efore, tangent cylinders, tangent cones, etc. The 
family of tangent cylinders constitutes a two-parameter family 
of siurfaces which might equally well be taken as a complete 
integral. A tangent plane would then appear as an envelope of 
a one-parameter family of these cylinders. 

Example 2. 2-(l-}-p--t-g-) = r^. 

Tills equation was obtained (art. 6f, ex. 1) from the primitive 
(r— o)--f (y— 6)®-f 2- = r - which is therefore a complete integral. 

t Tliero cxisl equations of particular forms which possess integrals not 
included in these three tyites. Such integrals are termed epectal. CC. Forsyth, 
DiJfcjmlial Equations, 1014, p. 379. 

4313 
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It consists of all spheres of radios r with centres in the zy plane 
The singular int^al is readily fonnd to be s*— r* = 0 which 
consists of tu’o planes parallel to the *y plane For the general 
integral we consider those spheres for which b = ^(a)> that is, 
we unagine a sphere to move so that its centre traces the cnire 
y = ^[x) The surface enTdoped by the mewing sphere is a 
tubular surface All su^ surfaces are included in the general 
integral Again we may pich out fitnn the general integral any 
number of complete mtegrals, the set of right circular cylinders 
of radius r with axes in the xy plane is an example 
Exa>otj 3 5 * = 4? Solve and interpret geometrically 

69 The method of Lagrange and Charplt, We consider 
here a general method of obtaining a complete integral of a 
partial differential equation of the tot order, 

F(*y,sps)«=0 (I) 

At any point of an integral snrlace of (1) we have 

<fe*p(fc+gdy (2) 

If functions p and q in terns of x y, and s are known which 
satisfy (I) and make (2) an integrable equation, then the solu 
tioo of (2) IS also a eolation of (I) Two such functions p and q 
could satisfy not only (1) bat a second independent relation 

“ a (3) 

The function O u subject only to the restriction that (1) and 

(3) be an integrable system and that — Our 

vp vq ^ ^ 

problem now resolves itself into that of obtaining a function 9 
which wiU answer our purpose To this end we differentiate 
(1) and (3) partially with regard to * and y In the following 
equations subscripts mdicate denvatives 

-= 0 . 

t Ttua eondittCA SKp&aa thal no ntatton independent of p »nd 9 exaits 
betwwa F «ndO 
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THE JIETHOD OF LAGRANGE AND CHARPIT 

By elimination of BpjSx from the first pair and of Sqj&y fi-om 
the second pair we get 

F^%-F^%MF^^-~F,<i>^)q+(F,%-F^%)^ = 0. 

Since — = — - = — , the last two equations give by addition 
8x Bx&y by 

Fj,%+F^%+{pF^+qF^)%-{F^+pF,)%~[Fy+qF,)<^^ = 0. 
This is a linear equation of the first order with dependent vari- 
able <I> and independent variables x, y, z. p, q. By the method 
of art. GG a solution is to be found from the system of equations 
dx dy dz —dp _ —dq ... 

F^~T,- pFj,+qF, - F^+pF, - F,+qF: ^ / 

For our present purposes the general solution of this system is 
unnecessary. One solution only is needed. The general method 
may now be summarized as follows. Obtain from the set (4) an 
integral <b{x,y,z,p,q) — a independent of F = 0. Solve F = 0 
and 3> = a for and q in terms of x, y, and z. Substitute these 
values ofp and qindz — pdx -\-q dy and integrate. The resulting 
solution has two arbitrary constants, a and the constant intro- 
duced in the last integration, and is therefore a complete integral. 
Example. j7--fg- = py—qx. 

The set (4) is for tins case 

dx __ —dp dq 

2p—y ~ 2q+x ~ 2p-+2q-—pti+qx ^ q ~ p' 

An integral is obtained from the last two terms, p-+^ = a. 
This and the given equation give for p and q the values 

ay-\-x>j{a[3^+y-)—a~} 

^ s+i? ’ 

_ —ax+y^{a{x^+y^-)—a^} 
x?+y- 
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^ t -‘«+yJo(a:*+y*)-an j^ 

**+y* 


whence g = atan '~+| 
where « = je*+y* or 


eXAUFLES Of CSAPTES XI 

Korz —For the linear e<iu(UMM>s find the general mtegraU for the 
ethers eotnplete integraU 

I fy+*)F+(*+*» - *+y 2 aj>~yj+»*--y»»0 

’ - 0 

7 ^(aoynif 8 1) ■» 5* 

9 10 (jr»-v‘-*»)p+2iy?i-Zr» 

11 ja+sy*-V(**+y’> 

13 *« 2p»+J5y+4rp»+2? 14 U+2-lp+(4«-»hr » 8*»+y 

15. (i’-y*)j>+(v'-t*>3 = 

18 p*+j*— Jjw— Jjy+l = 0 17 2i+F’+9y+3v* » 0 

18 Shaw that if one of P,(^, Ha xero m the equation Pp + Qq » H 
the characteristics are plane eurres id planea parallel to a coordinate 

19 Find tb« surfaces whoee tangent planes all pass through a fixed 
point (i^y^r,) 

SO Find the surface whose tangent plane makes with the coordinate 
planea in the first octant a tctrabedroa of constant TOjtime 

21 Find the surface whose tangent plane cuta off from the lym wtiw te 
axea mtcrcepts whose sum u constant 
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PARTIAL DIFFERENTIAL EQUATIONS OF ORDER 
HIGHER THAN THE FIRST 

70. General remarks regarding formation and solution 
of equations. In a partial differential equation of the second 
order there may occur two fii-st derivatives and three second 
derivatives. These are often denoted by single letters, 

_ Sz Sz 0% 0-2 0-2 

^ Bz’ ^ By’ ^ 8x-’ ^ BxBy' By- 

Similarlj', in equations of liigher order tliere may be four thh’d 
derivatives, etc. As m the case of first-order equations, partial 
differential equations of higher order may bo formed from 
primitives by the elimination of arbitrary constants or of 
arbitrary functions. The general solution of such an equation 
involves arbitrary functions. Many applications of the subject, 
however, require solutions satisfying given initial or boundary 
conditions, and such solutions may often be foxmd without a 
knowledge of the general solution. In some cases the solution 
is assisted by an integration with regard to one independent 
variable. In this case the constant of integration is an arbitrary 
function of the other independent variable (or variables). 

In the first four examples below, / and ^ are arbitrary func- 
tions to be eliminated. 

Example 1. y=f{x—at)+ij>{x-\-at). 

Partial differentiation gives 

I = -aS'{x-at)+aj>'{x+ai), g = a^r(x-at)+a^"{x+ai). 
Hence = a-—. 

This equation has a useful interpretation. Let x and y be the 
coordinates of a moving point and i the time; take a positive. 
When I increases by an amount A and x by a A, the function 
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flz~at) retuma to ita fomer t»Ioc Hence y — /(i— ol) lepre 
sents a curve which as tunc jncreaacs is continually displaced 
to the tight m other words, it represents a wave moving to the 
nght SiBuIariy y = represents a wave moving to 

the left The diffcrenliaJ equation represents this double wave 
motion and is much used ui the study of vibrating stnngs 

Eiamvli: 2 s = /(i)+^lCy) 

Exaupix 3 z 

Exuiple 4 ar-t-/(2) = j+^(y) 

FxASIPLE 5 Solve -^2- ss 2x 
ixey 

Exasiple 6 Solve ^+«*ra=0 (Treat as an ordinary 
Sr 

equation in z and f and inaKe the arbitrary constants fane 
ticms of y ) 


71 Particular solutions found by means of exponentials 
Exponential functions were found to be of great use ui the 
solution of ordinary linear equations with constant coeScienta 
In the cose of (be corresponding partial equstiona these funo- 
lions are often eonvement for finding particular solutions satis 
fying given boundary conditions A more general method of 
dealing with the equations of the following examples is given 
m art 75 and the subsequent articles 


EXiSITLE 1 — 4-~ = 


As a trial solution set s = It will satisfy the equation 

if F-fm* = 0 i and m cannot then be both real except in the 
tnviat case when both vanish Suppose we desire a solution 
satisfying the boundary conditions (J) z = 0 when z = 0 (2) z 
approaches 0 when y becomes mfinite From (2) m must be real 
and negative ray w = ~p Then I = d;ip and a formal sola 
tion IS z = e v*'(c,e'*’'+Cje *>*) la art SI we found, in con 
stdenng the equation iCyld^+P ^ « 0 that the solution 
Cj«^+v,< could be replaced by A cos ^+£sui^ Making 
the same substitution hereweh«Tez= « '’•'(dcoapz+Esinpz) 
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which is easily shown to satisfy the given equation. Since z = 0 
when x= 0, A must be 0. Hence z = sinpx. 

B^z Bz 

Exajeple 2. a-^ = with the same boimdary conditions 
dx- oy 

as in ex. 1, 

3 * 1 / 3 * 1 / 

TC xamp t.-r 3. -4 = a”—?. Find a solution which is finite 
dt- 8x- 

for all values of x and i and such that y = 0 -when x = 0 and 
also when / = 0. 

Bz Bz 

Example 4. — |-a — = 0. Find a solution which takes the 
dx dy 

value c* when y = x. 

72. Homogeneous linear equations with constant coeffi- 
cients, The complementary function. Ordinary linear equa- 
tions with constant coefficients were found to yield to methods 
of solution in which symbolic operators played a useful role. 
We shall find symbolic methods of advantage here also and 
shall denote d/dx by D, djdy by D’ so that, for example, D-D'z 
stands for dPzjd-xdy. The tjqie of partial differential equation 
for which symbolic methods are most useful is the so-called 
homogeneous linear equation with constant coefficients, 

=/(x,y). (1) 

A part of the argument of art. 28 may be adapted to the present 
case and the corresponding results will be merely stated. The 
general solution of (1) consists of any particular integral and 
the complementary function which is the general solution of 

Wo consider first the question of obtaining the complementary 
function. For the equation {D—aD'}z = 0 or p—aq = 0 La- 
grange’s solution gives z = F{y-f ax). Tliis fact suggests that 
F(»/+niix) will satisfy (2) provided wq is a root of the algebraic 
equation aim«-i-{-a. 7 ?i"---t-...+a„ = 0. (3) 

TIus conclusion is easily verified by direct substitution of 
Fiy+niix) in (2). If m-^, vu,..., m„ are the n roots of (3) and 
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are all distmct the geoerftl K^attoa of (2) ts 

' = •f’,(y+»ni*)+/i(y+«i,jr)+ +F„{y+n^f) 


ere the F s are arbitrary functions 

If (3) has two or mote cijual roots the solution (4) contains 
fence than n mdependent functions That it is not then the 
general solution appears from the following example which is 
a case in point 

(ly^iaDD +a n’*)z ^ 0 or {D~aD}*z===0 
By nsing *• =» (,D—aD )2 this equation becomes (D~aP')v = 0 
from which V a» i'j(y+a») Hence 

{D—aiyyt ss f,ly+ar) or p~ay » Pjly+tf*) {®) 
The cbazactensiic equations fenr (S) are 




tfs 


One integral is y+or 

whcnee-sl^—jr * B or - st-t — t 

■fiM) -fiM) .fily+a*) 

general integral of (S) u therefore 


/i(y+<») 

A second integral is goh from 
£ The 


sr sa ^-fy+or) or z « a-F,(y+flz)+ii(y+«) 

riV+®*) 

This result gires the terms in the aolution of (2) corrtsponding 
to a doable root a of (3) This suggests the more general result 
•ft bich ina> be directly renfied that corresponding to an r fold 
root a of (3j the terms in the solution of (*’) are 
j' ^Fi(j/-}-eu)+x' *Ji(y+ajr)+ +*^ i(y4-o^)+'S(!'+®^) 
Suppose now that (3) has* pair of coiyugate imaginary roots 
aii^ The corresponding terms m the solution of (2) ate 
J',(y+or+i^)+iJ{y4-wr— which are in general jmagmary 
If jowerer we set = 4+^ “id F* = where 4 amd 4 
are arbitrary real functions we get 
^ty+ar+«^)+^(jr+i»r-i4r)4- 
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which is real. If, for example, ^[u) = u-, tp{u) = the 
result is 

= 2{y+ax)^—2§-x?—2^x/{(y+ax)-+^x-}. 
Example 1. r = s. 

Exajiple 2. r+G5+9/ = 0. 

Example 3. {lP—2D-D’SDD'~)z = 0. 

Example 4. (D^-D''-D'+DD'^-D'^)z = 0. 


73. Homogeneous linear equations with constant coeffi- 
cients. The particular integral. The next problem is that 
of obtaining a particular integral of equation (I)^ art. 72. Fol- 
lowing the notation used in Chapter VI for ordinary equations, 
we denote this equation by F{D, X)')z — f{x, y) and a particular 

integral bv 7 z , J {x,y). The operators D and D' obey the 

following tiirce laws: 

(1) they are distributive over the terms of a sum, 

(2) they are commutative with constants and with each other, 

(3) they obey the index law. 


From these facts it follows that the meaning of F{D,D’)f{x.y) 
is unchanged when F{D, D’) is resolved into linear factors which 
may be uTitten in any order. The argument of art. 35 may now 
be applied to shou' that the inverse operator ljf{D.D') may bo 
MTitten in factored form or separated into partial fractions. 

This brings us to consider the value of — - — jjJ {x,y) which 
is a solution of the cqimtion 

{D—aD')z—f{x,y) or 2 ^—aq=f{x,y). (1) 

The equations of the cbaractcristics of (1) are 


dx 

One integral is y-foa; = 

<ilS 


_ dy _ dz 
-a ~f(x,y)‘ 

= A . For a second integral we have 

R 
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dz = /(x,A~as)dx wluch gires on integration a ralue of s and 
hence the formula 

(2) 

in which, after integration, A la replaced by y+ax 
ExiMTLE 1 (Ild-jy){B— 3D')s *= ye*. 

No* = j(A-3i)eiIs<~(A+3}er~3ze’ 

which, on the eubetitution of y+3;z for A, becomes 
Again 

(?+;».■ = j M+i+s)!**- M+2)<*+».- 

which, oa the eubeUtutioa of y—x for A, gives (y+2}e* The 
complementary foaction ta 3x) and the general 

solution, i «» &+2)e*+^i0'— *)+^r(y+5*} 
iSceond nulAod 

p+zrjiD-ji?)®** “ 

From (3) we get 

Hence a particular integral is 

I = ^b-i+aj+sj"* - jto+si** - (!(+2>!* 
ExAwrLB 2 {IA+iD2y+4Jy*)s = g?y 
ExiMyLs 3 (i)*-f-i>i)'— 6iy*)» — zcosy 

Examfle 4 (Z)*— iy*)i = ^amy 

The analc^ with ordinary lumat equations would lead us to 
expect that short methods may sometunes be employed to deter- 
mine a particuUr integral Hus is the case when f{x, y) has the 
form ^(oz-fhy) We have 

I>^(ax+hy) »= B^'{ox-i-by), 
i/^toi+hy) =» bi>’{ax+by) 
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So D-(f> ~ a^", DD'<j> = ab<f>', etc. Hence 

F{D,D’)<f,{ax+by) = F(a,b)<f,M{ax+by). 
If F{a, b) ^ 0, this gives by division 

F(D,D')i^^±M = ^M^ax+by) 
±(a,o) 

or, if we replace by 


(3) 


J|... J ^(<,*+65)P(»,+ij,)].. 

The case in which F{a, 6) = 0 is an exception to this formula. 
It is dealt with by means of the following example: 

(D—aD')z — xr<l3{y+ax) or p—aq = af<f>{y-\-cuc). 

A solution of this is found by Lagrange’s method to be 

iC** 

whence we get the formula 


1 af+^ 

^^,xr^(p+ax) = —^{y+ax). 

It follows that 

Example 6. {JF—DD’-)z = 8in(2a;+t/). 


A particular integral is _-^^^sin(2x+y) = ^cos(2a:+y). 

The complementary function is .Pi(y)+J^ 2 (!/+*)+^ 3 (y— «). 
Example 6. (Z>— D')2(jD-}- 2Z)')2 = e=^+v 
The complementary fimction is 

J’i(y— 2a:)+Ej(jf+z)+a:ii’3(y4-a;). 

A particular integral is 


(D-D'f D-\-2D' 


e*+v = 


[D—D’f 3 


21 “3~ 


= IxV+v. 
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Example 7 {Ifi-lrD»D‘—\iDD*)i = Sr+y 
Example 8 (D*- 10Z?2)'+25O^)2 = Iog(3j;-4y) 


74. Partial differential equations o{ mathematical phy- 
sics. Mention will now he made of some of the partial djf 
ferentnl equations which have been much studied on account 
of their use m mathematical phjaics For the deduction of these 
equations from physical hypotheses the student w referred to 
treatises which deal at length with these questions t The 
interest of the physicist is not so much in obtaining the general 
solutions as m obtaining particular solutions which satisfy cer- 
tain boundary conditions For (he cristence of such solutions 
the physical problems which give nse to the equations furnish 
convincing evidence 

ElaalKvtbralxmt We have already met (art 70) the equation 
for the transverse vibrations of a stretched string t 




(I) 


Boundary conditions in this problem have to do with the posi 
turn of the ends of the string and with the shape and velocity 
of the etni^ at a given uutant whieb u generally the beginning 
of thtt motion , the case of a struck piano string winch la initially 
straight and that of a harp string which is plucked and released 
are examples Siimlarly, for a stretched membrane which is in 
equilibrium along the ay plane and la made to vibrate trans- 
versely the equation of motioa is 


«* 

Examples of boundary conditions for this case will readily occur 
to the reader 

The flow oj heat If heat flows along a conducting rod the 
temperature at any point depends upon the position of the point 


t For example A G M eli*t«r, PaHtd D^mnt ai t quatuiM oJ ilalht 
maluoC FSyticJ 

t TKa conatATil a which Egurcs m ttua %x«l lira (oltoss uig equations depenOa 
on physical properties such as tlw demitv and tension of tho stnng or mem 
brsne, the conductivity of tlie conductor etc 
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and upon the time. If a; ho tlic absci.s‘?a and it the temperature 
of an arbitrary point of the rod, then u satisfies the equation 

^ = a-~. ( 2 ) 

e( dx- 

Boundary conditions here liave to do with the position of the 
ends of the rod and the initial distribution of temperature along 
it. If the flow of heat is maintained in such a way that the 
temperature at anj' point remains constant wjth the time, we 
liave ctijdt = 0 and hence dhi/cx- = 0. For a flow of heat over 
a flat plate or througli a three-dimensional conductor the equa- 
tions are generalizations of (2). Thus for a three-dimensional 

6t ^c.r- by cz-j 

If the flow is maintained so that the temperature at any point 
does not change with the time, it is called a steady flow and 
the temperature then satisfies 

d-u . B-u , e-u , . 

?3+5?+ST “ "• <*> 


which is known as Laplace’s equation. If, for a solid conductor 
bounded by a closed surface, the temperature of the surface be 
kept constant with respect to the time but continuously variable 
from point to point it is highly plausible that the temperature 
at any point witliin the solid will not vary with the time. In 
other words, a steady flow of heat will be maintained. The 
corresponding mathematical problem is that of obtaining a 
function which throughout a given volume will satisfy Laplace's 
equation and which will take on a given set of boundary values 
on the bounding surface. This is known as Dirichlet’s problem. 

Other examples of floiv. The same equations as for the flow 
of heat characterize the flow of electricit 3 % the variable u in this 
case being the electric potential. Another application of the 
equations in two and three dimensions is to the so-called irrota- 
tional flow of an incompressible fluid, in which u is the velocitj' 
potential. This last example is, in firct, useful in providing 
a visual analog}' for an}' ca.«e of flow according to the same 
equations. 
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The 2>euiontan poUvtu/J In tt»e ease of the gravitational 
attraction of matter a force function or potential function exists 
with the property that its partial denvatives dit ftr Znjiy dujZz 
are the x y and z components of the force at the pomt {x y z) 
In empty space this Kewtoaiaa potential function satisfies 1a 
place 8 equation TVithin a solid of continuous density p it 
Kitisfies Poisson e equation 




75 The equation ^ Founer senes W® shall 

now illustrate a method of modi use for obtaining particular 
solutions satisfying boundary conditions Consider the problem 
of the vibrating tlnng Suppose the atnng is of length I with 
its ends fixed that it is ortguiaUy at rest and distorted into a 
curve nith the equation y « /(xj Then rre require a solution 
of the follorring equation and boundary conditions 


rt* 


ex’ 


(>l»U, = o WsU-o W »!,, = /(*) 

We shall seek a solution of the form y XT where X 13 a funo 
tion of * alone and T of t alone The equation -will be satisfied if 


or if 


T^, 

«tr*' 


1 d^r o*d*X 
T dP “ X dx* 


Since the left member of this equation is indepcodcnt of x and 
the right member mdepeadent of t each must be constant If 
this constant is A we hare 


iTX Ay 


( 1 ) 


Suppose first that A = 0 We get y = XT’ — (Ciar+Ct)(c,<+«J 
In order that conditions (i) and (u) be satiafied. for all values 
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of t we must have Cj = Co = 0. The solution then becomes the 
tri\nal one y — 0 which cannot satisfy condition (iii). 

Next suppose A is positive. We get then 

y == (cie'“®+C2e“^*)(c3e''^+C4e-'^^). 


. In order that condition (i) be satisfied 

(ci+c^Xc^e-^^+Cie-^^) = 0 
for all values of t. Hence Co = — c^. Also from (ii) 

whence e“ — e “ =0 which cannot be true if A > 0. This 
solution therefore cannot satisfy the boundary conditions. 
Finally take for A a negative constant, —k-. Equations (1) 


become 


— + ~X-Q, 


■+k^T=^0, 


Ik k \ 

whence w = {A.cos-®+jBsin-a:i(C’cosI-<+DsinitO. May the 

\ c a J 

constants now be chosen to satisfy the boundary conditions ? 
From (i) we get A = 0; from (ii) £ sin A-i!/a = 0, and since B^O, 
klja — TTiTT OT k = mair/l where m is an integer. There is now 
no loss in setting B = 1. Then 

. mn moTT^ , „ . manX 

y — sm— xlCoos— ^t+Dsm— p-tl, 

dy . imr ( ^tnaiT . marr. , -.man mav \ 

i = k-xT®— ')• 

Now (iv) gives D = Q. Hence the solution satisfying (i), (ii), and. 

(iv) is y = Csin^^cos^^^^ for any value of C and any in- 
c c 

tegral value of m. The sum of any finite number of terms of 
this kind will also bo a solution satisfying the same three con- 
ditions, as n . . 

2 . i5r lair, 

CfBm—xcos——t, 

i-l 

and it is plausible that the same will be true for an infinite 
series of such terms provided the leqmsite conditions of 
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convergence hold There remains condition (m) When / = 0 
the last solution gnea . 

(2) 

Tlie qtic»tioii is now raised whether an arbitrary function /(a-) 
cm bo oxpics«ed in a sum of trigonometnc terms of this sort, 
or 111 on infinite senes of such terms This problem was pro 
posed b^ Fourier and its answer terms an important province 
of mathematical analysis Bnefiy it is possible means of a 
trigonometric sum (2) to exprees approximately anti by means 
of an infinite series to express exactly a function /(i) in the 
interval 0 < * < I provided /(r) is aubject to certain very mild 
restriction^ Assuming that such » development is possible we 
mar obtain the coefficients of the senes as follows Suppose 

/(!■) = Ci6m’p+c.8m'^'+c,sm^+ (31 

lest us get trr / I then /(e) « /(({/*) F(f) and (3) becomes 

^'(f) - f,sm{+fs*«'2f+c,ein3f+ (4) 

Multiply (4) throughout by smwf and then integrate the two 
Rieinbere between the hmits 0 and w We assume here that the 
integral of the intuute ecccee ts obtained by integrating term by 
term Then 

[ F{^)m\nidS =• e^ j sin/smnf +c, j ein2fsinn^(/f + 

’ * » (5) 

\OW, if ff MI 

j suiiif suimf i J [«>e(w— w)|— co8ln+m)f] d? 

) » 

— ^ r w”(w~*’*)f «n(nd-w)|'l* _ 

5[ «~m n+m J, 

Al» = 

Hence (o) reduces to 

J FC^lsiimf => c, Jw 
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whence ~ ~ J F(i)smnid^ 

0 

or, if the variable of integration is changed from i tax, 




. nvx j 
sin— r- dx. 

V 


Tliis formula gives the coefficients in the series (3) and with 
these coefficients the series y — ^ Cjsin^ajcos^^i is a solution 


of the given equation and four boundary conditions. 

The series (3) or (4) ivith the values found for the coefficients 
is sometimes called a half-range sine series. In like manner a 
function may be developed into a half-range cosine series. The 
further development of this subject will be found in books on 
Fourier series. 


76. Laplace’s aquation in two dimensions. In this article 
and the next brief discussions of Laplace’s equation will be 
given, using polar coordinates. This is intended to illustrate 
further the method of solution used in art. 76 and to suggest 
the coimexion between Laplace’s equation and certain im- 
portant ordinary differential equations of the second order. In 
two dimensions the equation is 


d^u.d^u 


( 1 ) 


For many problems polar coordinates are an advantage, in which 
case a; = rcosff, y — rBinO, or r = ^(x^+y-), 0 = tan-^yfx. 
Then (1) transforms into 


^ 1 a« 

dr-r dr'^r^ 88^ 


( 2 ) 


as the student should verify. Suppose we ^desire a solution 
M(r, 6) which is continuous within and on the unit circle r = 1 
and which takes on a prescribed set of boundary values on 
thisoixclo, ,.(!,«) =/(9), ,3) 

f{d) being periodic of period 27r. Applying the method used in 


4213 
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nrt 75 we set 
6 alone Tten 


or 


a iimction of r alone and © of 


dr*'^r dr^r*dr ' 


1 ^ 

<fr/ " Qdd* 


Since theleft memberis independent of $ and the right mdepen 
dent of r each most be a constant which we nay dencto by A 
Then «PO/ef5’+AO = 0 In order that 0 bo periodic A must be 
positive and m order that it have the period 2ir VA must be an 
integer AVe shall therefore denote A by n* Then 


The equation in J? is of the homogeneous type and gives readily 
Xae:Ar*+Br-" AhoO aeo$n$+^Bian$ Since tt is desired 
to be continuous whenrsO fimnstbeO We may take ^ 1 

and hence u <s r*(acosn5-l-b8mnd) This is a soluUoo for any 
zero or positive integral valoe of n In order to satisfy (3) the 
question is raised of whether a finite or infinite senes 
2 (o»cosn0+b^amn5) 

may be found to represent an arbitrary function f(S) This 
problem sgain belongs to the subject of Fouxiet senes 


77 Laplace’s equation in three dlmensiooB La this equa> 
tion we ehsU denote the independent vanable by V The 
equation IS . ... 

A homogeneous function of sr tr and z which satisfies this eqaa 
iKm is calkd a solid spherical har 
moDic For the study of sphencal 
harmonics it is an advantage to trans 
form the equation to spherical co 
oidmates If these coordinates are r 
{the distance from the centre) 6 (the 
longitude) and ^ (the co-Iatitnde) we have 
x^rcosfiBin^ p = rBm3sui^ z 



= rcc«^ 
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In spherical coordinates a solid spherical harmonic has the form 
where n is constant and « is a function of 6 and only. 

The transformation of the equation may be made in two 
stages by setting v = rsin^. For then x = vcosd, y = usin0. 
Hence, regarding 2 as constant, we have, by (2), art. 7G, 
cPV g=F _ 1 5F 1 aw 

0x2 ~ dv- 8v 002 

and tlie left member of (1) becomes, when V is regarded as 
a function of v, d, and 2, 

027 , 027 , 1 07 , 1 027 

02,-2+ 022+ j; 0y + e2 352- 


It remains to change from the independent variables v and 2 to 
r and Since u = rsin^, 2 = rcos^, we may again use (2), 
art. 76, regarding 9 as constant. Then 

0!;2 + 822 “ 0r2 +r 0r +r2 0^2 ’ ' ^ 

Finally the derivative 07/0U of (2) must be changed to the hide-, 
pendent variables r and Since v = rsin^ and r = 

4> = tan"2y/r, we get 

07 07 dr 07 0<^ _ 07 V 07 I/2 

eu dr 0U+ 0f dv ~~ dr ^{v-+z‘^y d<p l+r^jz- 


= — sm* 
dr ^ 


07 cos^ 
'dj, 


(4) 


Substituting (3) and (4) in (2) we get Laplace’s equation in_, 
spherical coordinates, ' 


027 2 07 1 027 cot^ 07 , csc2^ 027 

g^2+ r2 0^ 


dr- ' r dr r- d<l>'^ r- 002 ~ 

K 7 = r”u is a solid spherical harmonic, then u must satisfy 
the equation 


n(n-l)r"-2«+|7jr«-2M-f-L 

co^ a 

^ r2 0^+ r2 002 ’ 

p + COt^^ + CSc2^^+7j(7, + l)« = 0. (6) 
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A BoIatioQ of (6) 19 called & rar&ee aphencal bamooic If u is 
& fonction of ^ alone (4) becomes 

Thu equation transfonns mto I«gendie'a equation hj means 
of the BuhstitutioQ i 3= cos^ for 

du, ,du <f*« ,<f« . 

and hence jiu 

amV^-2«»^g+«(«+l)« * 0 

or (1— a^+nlnd-l)** = 0. 


which u Legendre's equatian (art 45) On account of this con 
nexion with Laplace's eqnaboa the 
Legendre functions are also known 
as surface zonal hannoaics 
Ihe tranafonnahon of Laplace's 
equation into e^lindneal coordinates 
u also osefnl That is, x rcos9, 
y m rsin z a 2 B/ sppljing (2), art 76, we get the tnns 
formed equation 

1 er 1 g«r 




Let US attempt to flod a aolntion of the form F <= SQZ, these 
letters denoting functions of r, 0, z respechrely By substituting 
and dividing by SQZ we get 

i. ‘**■^ 4 . J- ^ 4 . JL _ 0 

dr'^r'Q d0*'^ Z <fe* 


As the first three terms do not depend on z. while the last term 
IS a funcDon of z alone this term must be constant Setting 
2 — h* we get 2 =» C|e**+f,«^ and also 


Rdt* 


7+-t*^ = 


1 epe 
■© as* 


By the same argument as before, each member of this equation 
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1 

is constant. We shall set — — = m- whence 

0 

0 = c^cosTtiO+CiSinmO 

and = 0. 

ar® dr 

Changing the independent variable to t where kr — f,'we have 
dB_,dR 

dr dt’ dr"- ' df?' 

whence fs ^ 1 ^ = 0. 

dt^ qX 

This is Bessel’s equation (art. 46), Thus the study of Laplace's 
equation makes contact at this point with the theory of Bessel 
functions which from this circumstance are also known as 
cylindrical harmonics. 

EXAMPLES ON OHAPTEB XU 

Find solutions in the form of exponential functions for 1, 2, and 3. 

1. r+j r= p, 

2. r+s+p—q ~ 0. 

3. r— {— 2p+23 == 0. 

Find the general solutions of 4, 6, 6, 7. 

4. r—x = t—y, 

B. 6r+ll«— = 2x— y. 

0. r— «— 6f = xy, 

7. (D^ — 3D’D'+3DD'* — D'^)z = 8in(a; — y). 

8. Verify that the following functions are solid spherical harmonics: 

ax+by-{-cz} (a;’+y*+s’)”*; tany/x; x*— y'; 2s®— x®— y*. 

giy aiy gzy 

9. Find a solution of ^ + ^ + ^ = 0 in the form V = XTZ, 

these letters denoting functions of x, y, and s respectively. 

10. Represent by a Fourier sine series the function x in the interval 
from 0 to «•. 

11. Represent by a Foiuler cosine series the function x in the interval 

from 0 to w, . 



HISTORICAL NOTEf 

Ik the course of the foregoing treatment the names of vanous 
mathematicians have occiured in connexion with particular 
equations ot methods of eolation A bnef outline of the histon 
cal eetting of the subject wAl now be giTCn 
Ihilerential equations made their appearance as soon as the 
calculus had been invented in the latter half of the seventeenth 
century The formal methods of integration which have been 
dealt with in this hook were nearly all developed during the 
next hundred years The two founders of the calculus were Sir 
Isaac Newton (1642 — 1727) and Gottfried Wilhelm Leibmtz 
(164^1716) both of whom contnboted to the study of dif 
ferenbal equations The following quotation &om Ball aHulory 
of ilathmaHca shows the trend ofNewton s ideas The object 
of the second or inverse method of fluxionst u fcom the Simon 
or some lelationa involving it to determine the fluent or moro 
generally an equation being proposed exhibiting the relation 
id the Susons of quantities to Snd the relations of tho o 
quantities or fluents to one another This is equivalent eithc.i 
to integration which Newton termed the method of quadrature 
or to the solution of a differential equation which was called by 
Newton the inverse method of tar^ents Newton classified 
some simple types of equations of the first order and made use 
of power senes m ohtauung solations The contnbutions of 

Lmbnitz included the solution of the equaticmy^ = 2^(x)F(y) 

oy 

by quadratures which amotmts to a separation of its vanables 
the solution of the homogeneous equation of the first order and 
of the Imear equation of the firat order 
Contemporary with Newton end Leibnits end following them 
was a group of matbemsticians who used the new methods of 
t For fuller actounts of tlifl hwtory <if ths subject the reader la referred to 
IV W R Sail, A Short Uutor^ oj JVolAemottef I«ndon F Cajon Htototy 
of ilathematier hevYoA or E !■ Inee (M nary Vi^mrUtal on* 
I«cdon, 1927 p 629 

t Tb Kewton s nomencbtore tbe tenos fluent and fluiion were used for a 
function end ts denvative or diflermtialcooffic ont 
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analysis 'vrith great power. Differential equations were studied 
by these men on account of their intrinsic interest and also 
because of their occurrence in problems of geometry and 
mechanics. I’oremost among this group were the Swiss mathe- 
maticians, James' Bernoulli (1654-1705) and his brother John 
(1667-1748). Among the problems solved by James Bernoulli, 
who was the first to use the term integral, are those of the form 
of a catenarj', the form of an elastic rod fixed at one end and 
acted on by a force at the other, the form of a sail filled with 
wind, the curve in W’hich a pendulum should swing in order that 
its period should be independent of the amplitude, etc. The 


equation P?/ = Qy" was proposed by James Bernoulli and 


solved b 5 ’- liis brother John and by Ijelbnitz; John Bernoulli 
made explicit the method of separation of variables; he solved 
the problem of orthogonal trajectories and the problem of the 
brachistochrone, i.e. the curve along which a particle falling 
under the influence of gravity will travel from a point A to 
a lower point B in the shortest possible time. 

The Italian, Count Jacopo Riccati (1676-1754), was first to 
point out that a second-order equation which does not contain 
the independent variable explicitly may be reduced to a first- 
order equation. The equation which bears his name is 


% = A(x)j,HP(x)y-l-C?(x), 


of which a special case is = bx'^. Daniel Bernoulli 

ax 


(1700-82), the son of John mentioned above, showed that this 
equation is integrable in finite form only if m is of the form 
— 4l7(2h±l) where h is zero or a positive integer. For other 
values of m, Riccati’s equation is the simplest example of one 
which cannot be solved by quadratmres. 

The device of using differentiation to assist in the solution nf 
a differential equation is due to Jean le-Rond d’Alembert 
(1717-83). The equation a:^(y')-}-g^(y')-f xCy') = 0 is known as 
d’Alembert’s equation. To the same period belongs Alexis 
Claude Clairaut (1713-65). The equation known by his name 
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IS a paiticalar case of d Al^b^ s Iq its study Claur&ut 
obtained singular eolntionst and sbowed that they were not 
included in the general solutions 

One of the most powerful of the eighteenth century mathe 
maticians m the use of caicolua was tieonhard Euler ( 1707-&3) 
Among his coninbutions to differential equations was an ex 
plicit discussion of integrating factors of which isolated examples 
had occurred before Euler set up types of equations which can 
be Bolred by integrating factors ofpien forms In his work be 
made &equent use of the change of vanablea Much of the 
theory of linear differential equations is due to Euler Be also 
extended and unproved Newtons method of integration in 
senes 

The penod of development of the elementary methods with 
which this book deals may be said practically to have closed 
with the work of Euler By the middle of the eighteenth century 
the formal methods of solution were known and the attention 
of mathematicians turned to more fundamental questions One 
result of the labours of the early analysts was to show that the 
method of quadratures is insufficie&t for the solution of other 
than specut cases A more modem point of view was initiated by 
Joseph Louis Lagrange (ITSA-ISIS) He introduced the adjoint 
equation which is satisfied by any integrating factor of a given 
equation. He abo developed the method of vanation of para 
meters Problems concermog the existence and nature of func 
tions satisfying differential equations of vanous forms were 
mvestigated in the nineteenth century notably by Augustm 
Louis Cauchy (1789^1857) From the tune of Cauchy on the 
literature of differential equations becomes increasingly volu 
mmous 

The subject of partial differential equations presented more 
difficulties than that of ordmaiy equations and progress was 
slower The impetus for their study came first from physical 
problems The problem of the vibrating stnng was considered 

t The geoiaetn«&l s grqfimnwi adiiticpta 'wu pennted <nit hj 

bat tba funbsr snslysis of Um iubj«ct left to tbs Uttec ti»i£ 
oS the nineteenth century Tbe (leatmeot vow used foQovi the vork ot Arthur 
Ceyley (1S21 
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early by Brook Taylor (1686-1731) and later by Daniel Ber- 
noulli and Euler. A more significant contribution was made by 
d’Alembert who gave the general solution of this equation and 
who encountered partial differential equations in problems on 
the motion of fluids and the theory of winds. The work of Euler 
and d’Alembert in this connexion links up with that of Jean 
Baptiste Joseph Fomier (1768-1830) who gave the equation for 
the flow of heat. The first successful attempt to deal with the 
subject systematically was that of Lagrange to whom is due 
the solution of the linear equation of the first order and the 
further method known by his name and that of Charpit. Pierre 
Simon Laplace (1749-1827), in connexion with the idea of 
potential, introduced the equation known by his name, 

dx^ Sl/2 1- 022 

and also the spherical harmonics which are among the solutions 
of this equation. 

The further development of the subject, though of great 
interest and importance, goes beyond the subject-matter of this 
first coiurse. 


4212 



ANSWEES TO EXAMPLES 

CBAPfEB I 

Alt. 3 pages Ex i It) dyjdx => m (a) f/ xidyldx)+{> 
lw)>*(p*+l)>=f* r*»){*-»)*&»»+I)»r*a+p)* wlierep-rfyAfa 
Art.) p.g.r E« ! (,_S)g+i +^)'_0 

- -i@r- {■+§ns 

£xaEDpIes at end Page 8 

. 3-0 = (,- 4 )‘ -+&’)■ -S- 

4.yds-ixdy=^0 6 (i-sWl+P*) •* <*+n>l* 

« sW fcp 1 (p*-i»V+ip-»o 

a a =• 0 

10 j^ax-l)dy+(li*-3}^-.c.(-3}<ic = 0 

11 (*+yp)(jp-j) “■ ij> 1* tp*+i)(**+y'J* = <(p»-p)* 

IS (H-Sp)(l+p»+ji5)-J}(*+sp) — 0»l«Tej=«iJV/it* 

U *yj+*p*— VP • ^ 

CBAPTBB tl 

Art # page 10 E* 3 **— i*V+4**v*+V* • « 

£x 3 e*+»*+y8a«'{'MeyB< 

Art 7 page 11 £x 2 

Fx 3 a^(l+y)-< Ex « («-3Xy+3)’(y-J)' -*<*+*) 

Art. 8 peg? II Ex 2 logVfa^+V*) ^ ^ ’(y/zl+a- 
Ex 3 x'ti^+Zv*) » e Ex. 4 . **+y* » «*» 

Art S pagsi: Ex 3 <3i+y-6)* - c<2*+y-3) 

Ex i Sz— 10y+lll<>g{lSx-'Sy— 13)« a 
Ex 6 <x+yP+fl*-2y+3iog(»+y+l> = « 

Art 10 page U Ex. 4 S*V*~6®*y ' =•« Ex 5 «;•— I ^txy 
Ex 6 *V+y*=«* Ex 7 Sx^«jr<+**V*=»e 
Art 11 page 16 Ex 3 6ry=.i»+« Ex 3 y^V+**) "= *+*• 

Ex 4 ly-l)a»;»»e. 

Art 12 page 18 Ex 3 y(e— x) =>**“* 

Ex 4 **lv'(l+y*)+rf = y+e 
Ezamplee at end Page 17 

1 y-«=>c(I+2y) 2 logx— <y+ 2 rJ*/(y+*P = « 

5 y •= a* (log*+e) 4 2y’+2*y+** =• w 

6 ar-i ■= eoaP(logcoa9+e) 0 *y+*y+**~y* “ * 

7 2y(*+a)* =. 6{z+<»)*+«. * s--2y 41oe(®r~5y+ 10) e 

9 »y •= ee*(ry+I) 10 (f*+l)8in^ =• c 11 (y— l)tan* «s e— * 
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ANSWERS TO EXAMPLES 

12. tr^l*+y = C. 13. tan»«+icottf — Bin(tf+i)— = c. 

14. 2xy+l = Cl/’. 16. a^+2x^ = cy. 10. Saic’i' = e®>'+c. 

17. Bm(y/x) = cx. 18. 6(a**+6*i/) = atan(a5a;+c). 

19. Vi/logs = (logx)»+c.. 20. log(x*+j/’)— 2/j/ = c. 

21. ylx = tan{c-a/V(x»+ 2 /’)}. 22. j/* = cc**/*-!. 

23. x’+y*+4a:i/+2x— 8i/ = c. 24. xlogy = xe*— e®+c. 

26. y^€!~^l^ = 2 J dx 4"C. 

26. = Boc0-f-tan0—ff+c. 

27. (l+x)(l— x+y) = cc''. 28. x’/’y’+^V^'’ = c. 

29. l()g(./+2)4-2tnn-‘(i/-2) (.r-.6) c. 30. e» = (o— y’)tftnx. 

CHAPTER III 

Art. 14, pngo 20. Ex. 3. (y— c)(y— 2x— c)(y— 4x— c) = 0. 

Ex. 4. y = sin(x+c). Ex. 5. (2y— x’— c)(y— ce*) = 0. 

Ex. C. (y— cx){y’— i’— c) = 0. 

Art. 16, page 21. Ex. 3. 2y = x’+cx*. Ex. 4. xy— c’x+c = 0. 

Art. IG, page 22. Ex. 2. (y+l)’(p’— 1)’ = c, with the given equation. 
Ex. 3. 2x+y’ = cy. 

Art. 17, page 23. Ex. 2. y = cx4-2c/(c’+l). 

Ex. 3. y’ = 3ex+c’/’. Ex. 4. logy = cx+c’ " 

Examples at end. Page 23. 

1. (y— ee’*){x’+4y— c) = 0. 2. c*x’+o’ = 2yc. 

3. y’ = cx+cV4. 4. cWr = with the given equation. 

6. cv = cc*4-c". 6. X = asec(y+c). 

7. x+c — 2f-nn”’jD— 3/p, with the given equotion. 

8. 4yc = (xc+4m)’. 9. (2y— x’— c)(2y+x’+c)(xy— cy— 1) = 0. 

10. y’ = cx’+c’. 

11. log(2P— l)x— 1/(2P— 1) = cwhere P’ = (4y*— x’)/4x’. 

12. ly=+yv/(x’-i-y’)}/xJ+log{yJ-V(x'+y:)}/c.T-' -=• <1. 

13. y’ = C.T+C’. 14. 4y-f-2x’ = 4cx — c’. 16. 2y = x(ce*— c“*e~*). 

16. (y— i-c)(y— ce*)(l— x+y— cc-*) = 0. 

17. r = aBin(0+c), in polar coordinates. 

18. x+p — - CC”’’, with the given equation. 

19. y> = cx’-a*c/(l+c). 20. x’+y’+2V2cx+c’ = 0. 

CHAPTER IV 

Art. 19, page 27. Ex. 1. x’+4y = 0, a singular solution. 

Ex. 2. y(y— x) = 0; y = x, a singular solution; y = 0, not a solution. 
Ex. 3. y’{i’— y’) = 0; x— y = 0 and x+y = 0 are singular solutions; 

y = 0, not n solution. 

Ex. 4. 9y* = 4x’, a singular solution. 


T' 
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Art 20 page 28 Ex 1 **+J^ =» 1 
Ex 2 For ei 3. *rt 19 **+^+2«af+i?/2 =. 0, 

For ex 4 art. 19, 9v— a O 
Examples at end Page 33 

1 y ^ e*+a/# y* 4a* 3 (y— e*)*-pa*e »» 0, 4*y = a* 

3 y— e* = ^(1— e’)~eeor^, jr=* *Bi* 

4 y+c = e*«, y = *+l S 27(y— «r)* •* 4<*, y = ** 

0 (3y+e]' s 2ar*. Oy » ■a 0, eusploeus. 

7 (y+c)‘ •= (*+«)*. *— y =• */t^> *^Sf. eusp locus 

8 x*+y'~3c(x+y]+c' ** 0 0 y B 0, * >9 y, tao locus 

9 y-e* »= «e/(e-l) (ai+y-o)* = 4iy 

10 y* s (x— eP y 0> cusp locus sztd eovelope 

11 y<9e(x-e]< 27y <9 4*< 

13 c<+y''4o*+5e'~l e0.**/6+y*s 1 

13 y*+<V - e. W “ 1 

14 e**+y' •• e/(e— l)».(**+y*'-I)*—4**y* m 0 (conic* touching the 

four ndea of a squveV 
18 y—e» s e*/»+j(*/* m «•/» 

16 9(8y49/«l'«S98«s*.8*+2y«0 

CHAPTER V 

Art 39 page 39 f-% 4 100 000(12}" 

Ex 6 (Iog30)/log3 «• 4 33 hour* 

Act 34 page 37 Ex. 3 80(l>«->'') lb 
Ex 3 *« V(J«/»'W-2)/(3'*A»t+«-lJ Ex 4 74 nanutc* 4e’C 
Art 25 page 39 Ex 2 y* ^ 2aT4-e Ex S'yece'i* 

Ex 4 Ex 6 r * ee*H«'-H e, g r.fcaaie 

Art 29, page 41 Ex 2 y » cx* Ex 3 3’+Uogr)* « e 
Ex 4 log^x'+y’J+tan *y/* = **f^= 

Examples at end Page 44 

1 x*+y’ = ey 2 **+y* = cr 3 (x~y)'-3a(«+y)+a* = 0 

4 *y="A/2 5 r = 8 «*+y* = 2logM 

7 2**+y* = e 8 !c‘-y»=.c 9 (ar-ej'+y* = c»-l 

10 CiKle centre midiray be t ween the given pomte, radius lislf the 

_giveo sum 

11 xV('>‘+4:'}+S'’A' 2 (tAegiven pomta (a,0), (— o 0} the given 

rwjtangle =» t>) 

12 f = «(l-cos^J 13 T-=e^P 14 r«e(l+cos^) 

16 i'/»+y‘'' =■ o*/* 18 19 4 98 per cent 

20 tS", 21 1 22 10/» 23 29’ C 

24 1 09r = 8 + 100(0455J» 25 y = *tana-*'(y/2v»>oc»a 
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ANSWERS TO EXAMPLES 
26. c*’— y* = a-cjil+c) -with fixed points (—0,0) and (a.O). 

•27. Portions of two parabolas. 

CBAPTER VI 

Art. 29, page 47. Ex. 2. y = CiC*+Cje'’*+CjC~‘*. 

Ex. 3. y = 

( V7 . V7 \ 

.4 cos — *+ Bsin — a: j. 

Ex. 3. y = c'*(Acos6a:+Bsin6a;). 

Art. 33, page 61. Ex. 2. y = c"(CiX+Cj)d-e*®(C3®+C4). 

Art. 35, page 64. Ex. 2. y = 6/2+CiC“*/*+Ci®~*®- 
Ex. 3. y = — a:/4+CjC“+Cje~“. 

Art. 36, page 66. Ex. 2. y = e^/T+CiC^I^+Cgtr*^. 

Ex. 3. y = ^=‘/(b^—a^)-i-Cie^-i-e~^J^^gC03~ax+CgBm^atj. 

Art. 37, page 67. Ex. 2. y = e"“(c4a:4-c,— sinr). 

/ V3 V3 \ 

Ex. 3. y t= a^*/3+C4C®+c~*/*lcjC03 — x+Cjsin — *1. 

Ex. 4. y = e“(a:V2-*+l)-fCie*+c,e‘®. 

Art. 38, page 69. Ex. 1. y = ^8ini*+Cie®-f Cje~*+Cje**+C4e~»». 

/ “VS ^3 \ 

Ex. 2. y = sin*+6^^/*(<!iC03— a:+c,sin-^a;j. 

Ex. 3. y = a;(Bin3a:— C08 2a:)4-CiC0s2x+Cj8in2a:+CjC08 3a:+C4Bin3*. 
Art. 39, page 60. Ex. 2. y = j(2r*— 10r^+16r)+CiCO32aT+c,8in2aj. 

Ex. 3. y = 2r’+3a;*+27a;4-Ci+®*®*+®s®*** 

Examples at end. Page GO. 

1. 10+oo3-^a:(Cje*/*+c,e~®/*)+sin^x(Cje*/*+C4e“®/*). 

/ 2 

2. ■^(8sin4a:— 16cos4r)+Ci€r*4-c,xe^. 

3. *e**4-cie”+c,c“”/^ 4. x+l+Cie"®+c,e®+C3a»®. 

6. c“’(a;^— 6r)4-c^(Aco3a:+Esina;). 

-1+VS_ -1-V5 

6. K«+8a:3+48x*+168r+288+Cie»+c,e 2 +c,e~2~®. 

7. ic“®(8in»+cosx)+Ci6“*®+<rt*^c,co3^®-f CjSin^xj. 

8. e“®(Ci+c,a:+C3a:*+C4a:^4-ar‘/24). 

9. (x cos 6x)/20 d- (a: cos a:)/4 -J- Cj cos 6®+ Cj sin 63;+ Cj cos a: d-C4 sin x. 

10. c“r“+V(o+l)(o+2)+e“*(CxX+c,). 

11. Jz’c-*®+Cje-’®+e-»(cj*d-C5). 

12. xcos2r+CjCosx+CjSinx+CjC032xd-C4sin2x. 

13. e*®/’(xd-Ci)— c“‘‘®/*(xd-c,). 14. e®(c4d-c,x— 2cosa:— rsinr). 
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15 c,«to+c,«-**-(l+3j*+3i*} 

16 ^Rinx— co«2«+]/5+C|eiaVl»«+^CMV5f 

18 (zeosx]/3+(2uiix)/9+e,eoa2z+^«in2i! 

19 —x*oosx—thmz~2xeo»z+uax’^eoax+e^ 

SO a^n®+*cosi+ei*m*+e^«)«* 

21 — **co8*+x&maf+(ei+ej»Xc»*“«*+®**“**) 

22 **-l0i>4-48J:*-ia2*+e,+e»«-«+«»«^ 

23 l«-*'{3ieia2*+eQ#44!j+*'“(>l«»»**+®»®2*) 

24 f,e-*«+e,«-“f»-t-*(l»ail*+l^+8) 

26 X*— Sx*+3»izkS+ei+e,z4-e,e(iex-f-C«Bm> 

26 l/6+(coa2c)y3— (co«iz}/6-)-«iCMtV3£+r,iijt2V3ir 


CHAPTER Vtt 

Art 40 p«g9 83 E* 1 y«=e,*+<»^ Ei 2 y «» ar*(ei+eilog*) 
Ex 3 y B z(Aeo«}agz-f-Bclologz) 

Art 41, pags64 E* 2 t/ •m e*lx •— 2*^>+e,»**+e,*r"* 

Ex 3 y ■« **/30+iF/ZO— t+ar*{^+^«odog»+4|ni»log*) 

Art 43 page 67 Ex 3 4«*>ywe(X*+et«4-ri 

Ex. 4 y«e,»+«|*r^/**‘ 

Art 44 page 73 Ex 6 

+‘^»(>+,V+S^6'‘+ST5TTi''+ ) 

Ex 6 

y - a(l 2+3 3;z+3 ix‘+4 6rt+ )+ 

+6<4r(3 8+* 7x+7 9rt+B 11*»+ ) 


+b(i+=: 


-*•+- 


-3-14 6 


*•+- 


- 11468 . 


at end. Page 74 
1 7^6+ei*+e,*» g -(*+I)‘/6+Cj(*+lJ‘+c/*+l) » 

3 (**-4rt)y = i*+ei*+«, 

4 (s'log£)/l2+e,x'+s~'{e,«oe(V3loga}+e|nQ{V31ogz)} 

e -a*/2+e,+e,**+c,*‘ 6 o(*+X^/4) 


- 1 - 


2 M»‘+l)^ 3 * 2 {»+lX»+*> f 

8 a*+6(l-£*-rt/3-**/6- ) 

9 (ox/4+76/12)/a*+Ct,/(ar+»)+i^+6)-* 

10 ®*y*— ay* B cr^-c, II I6y ». -^fclngx+e, A*+f,ar» 
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ANSWERS TO EXAMPLES 

CHAPTER VIII 
c®4-l 

Art. 40, page 70. Ex. 2. y = xjc log{l+c*)+Ci. 

c 

Ex. 3. j/ = — 6m(x+c)+Ci®+c,. Ex. 4. y = a:’+cloga;+Ci. 

Art. 47, page 77. Ex. 1. y— clog( 3 /+c) = x+c^. 

Ex. 2. cy— 1 = c*(a:+Ci)V4. Ex. 3. log(cy+c®) = c(a:+Ci). 

Art. 48, page 77. Ex. 1. y = cBm{az+Ci). Ex, 2, 2y“^ = (ox+c)*. 

Ex. 3. 2e-''/* = c— a:. Ex. 4. y = sec(a:+c). 

Art. 40, page 79. Ex. 2. y = — a:’/4+Ci*e*’+Csa;. 

Ex. 3. y = ar'(Ci+CjC*). Ex. 4. y = — e“y2+Cie^+Cje* J e*+*' rfx. 
Art. 60, page 80. Ex. 2. y = xe*/*{ci+c,®). 

Ex. 3. y = (Ci«‘+Ciar5}/(a:— 2). 

Art. 61, page 81. Ex. 2. y = (e*+e**)log(l+e“®)+Cie®+c,e**. 

Ex. 3. y = a:Han“*x— }x*log{l+x*)+CiX*+Cjar*. 

Art. 62, page 82. Ex. 2. y = (Bin?x)/4+Cisin*a:+CiOsox. 

Ex. 3. y = -x*/o>-12x/a*-24/oHcie“''®+c»e-“''®. 

Examples at end. Page 82, 

1. y+Ci = logBec(»+c). 2. y = e®{c(x+l)’— x/4+Ci}. 

3, tany/2 = ce’'*®. 4. y = — x*+cx| erV^dx +CiX. 

6. y s= c** J ar‘6“*V* dx + 0 ^**. 6, y = (x+a) 712 +c{x+o)‘+Cj. , 

7. y = ear* J dx +Ciar*. 

8. y s= ein* J e®(cso*+ccsc*x) dx -f Cisinx. 

0. y = e’®V>{ciC03 V3x+c,ain V3x). 10. x+c^ = log(e*/''— c). 

(i) y = — Bmx*+cie®*+c,e'-®', (ii) y = {Ci{l— x»)+c,x}/(l+x*). 

12. y = e®+Cie®/x+cJx. 

13. (i) y = Ci(2x+7)+c,e**— i+e*® J e”*®(x+3)log(x+3) dx. 

(ii) y = Cxe®+e® J c~*(x+l)>(x+c,) dx. 

CHAPTER IX 

Art. 63, page 86. Ex, 2. (x— o)*+(y— 6)* = r“. 

Art. 64, page 87. Ex. 2. A parabola. 

Art. 66, page 89. Ex. 3. a = (y/fc*){fc<— l+e-*<). 

Ex, 4. 2 seconds nearly; 19'6 feet per second. 

Examples at end. Page 92. 

1. C3in~\/(v/r)- y/(t/(c -y)' x+c,. 2. 2y = Cie"+cf*e~®®. 

3. y+c, - logiifc^=^>_V{l-(c-x)*}. 

4. y = e®/®+e-®/‘‘. 
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6 2(y+l;)= He*/*+€^'*)jfcufv«t«iebe»y-axiaatongva 

6 y = tz*/83V ’ y = 

8 » •» + 4w/Vn. sr. 16/V17. 

B (-|<ir/^(ss(*-7-iaiog^)) 

10 v= V(s/lK«*>'w>«-l)/{«»^+lJ 

•-«,«= glogW+e^wW+e-"^**') 

n i - i^y^. «> = "^(e^ >2 a*(a*-**| = b* 

CHAPTER X 

Art, 58, page SB E* 3 y = <».*‘ = ^jr 
E* 4 («+o)>— >,(z+e) 

E* B **—*‘*> *1, **+*••» «,y* 

Art 69,pegefi3 Ex 2 y «» <ie»+e,r»,» «,e»-(,r* 

Ex 3 y *» 3eQS*+<|«'“*+«j«*^**. 

Ex 4 x•.S(^!+t9^^e^■»•l•.vm6t+S+$e^0-»’/•+e^e• 

Art 82,pt2«l02 Ex 3 x()r+t) 

Ex 4 x*y— 3M'f3y2* K « Ex. 8 rfy+yjt^e 
Ex 0 *•— ay+y* •» «. 

Art 03, page 103 Ex 2 The i&t«rsectuMi 0 witb the hypnbolMd of the 
cylinders y* » *i(*— i). 

Ex. 3 The isteiscctioas with the given wfare of ihe qaadiics 
iry-* = e 

Examples at end. Page 103. 

1 *»— y* *. e,f, jf*— t*=«e,r 

2 x+y+* =* o, *’+y*+*' “= 0 

3 a»— jt* ,= 0, <*— y)s*P »» h 

4 ;cV— ‘><eec(o+x)+tan(04-O = ^ S y = ar, «— 2 = 6«* 
0. y= I/12-*/2+lIe*/12-3e,a-**/2-,^«-*», 

* = I/36+x/0-«»/a+<fi«-"+e,a~»* 

? y 2 xe<wz— {3xsu>s)/24-(se()e*-(-c,smz, 

2 = 4zsinx-'(&ccos7)/2d-e«eosx4-c,auiz, where 

5f, = -I3(i-c,+I,«f, =.«,_l2e,~I/2 
8 * = « ‘i»st+Ci#“*'+e» y = 

0 (y-*5 (t-x) *= « 10 (*-e?+y+»* = 1 

11 (i*+y»+2»-fXx’+y’— **-«>=»<> 12 z/y+y/i+a/x = c. 
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ANSWERS TO EXAilPLES 
13. *•+!/»+!:* = C3. 14. xifHx-{-z)(y+z) ~ a/(o+z®). 

16. (x^-y)l(y*-z) = a. 16. (a:+3)/y+{y+z)/a; == a. 

17. s»+J/*+z* = 3s, j/ = 2 . 

18. The intersections of the ellipsoid with the spheres = c. 

19. The intersections of the cylinders with the parabolic cylinders 

j/* = ax: 

21. xy->rXz->rXt-\-yz+yt+zl = c. 

22. an/*+3x*2t-2i/’2*-t/t’+42t* = c. 


CHAPTER XI 


Art. 64^ page 106. Ex. 2. (a:+zp)(l+g) = (y+2g)(l+p). 

Ex. 3. 2 = px+qy+pq. Ex. 4. s = pq. 

Art. 65, page 107. Ex. 1. px = qy. Ex. 2. bp = aq. Ex. 3. yp — xq. 
Ex. 4. xp+yq = z. 

Art. 66, page 110. Ex. 3. = 0. Ex. 4. fb{,yfx,z) = 0. 


Art. 67, page 111. Ex. 1. z == aa:+j/^(4— o*)+c. 

Ex. 2. logz— a/z = x+try+fi; Ex. 3. z = flw:+6y+4»— 6. 

Ex. 4. 2 = a(a:’+j/*)+6. Ex. 5. 2 = aB-hoy/(o— l)+c. 

Ex. 6. ®+o’i/+6 = alogz. Ex. 7. z+a7(l— »*)— otan~V = 6. 
Art. 68, page 114, Ex. 3. A complete integral, 0*2 = (®+(jy+b)*. 
Examples at end. Pago 116. 

1. E{{®+y+2)(y-«)’,(y~®)/(2— y)} = 0. 

2, E(®’+y>+ 22 ,*i/) = 0. 3. F{xila—xjb,xja—xjc,xja~z) = 0. 

4. {z—bx—by—c){z—b'‘x~by—c) = 0. 

6. 2= a(®+l)V2+(3/— l)V2a+6. 

6. 2 =’sin{(®+ay)/(l— n)4-b}. 7. 2 = alog»secy+6. 

8. {z—b){z*+d‘—(x+ay+by) = 0. 


9. 2 = aten-»*+Jologi^+6, 10, E{(s’+y*+z*)/z,y/2} = 0. 

11. z,y/®) = 0. 

12. ®j),®i+®j-l-*j,2} = 0. 

13. 2 * = tw;+6y+o*+b. 14 . F{x^~y—z,xy—z’‘) = 0. 

16. F{[x-y)l{y-z),{y-z)j{z-x)} = 0. 


16. 2z+b = a:’+y’+2V(a^+a-l)+(a-l)log{®+V(a^+a-l)}+ 

17. . 1. +W(>''-<.)-«!ogU,+V(!r--«)>. 


19. F{{y yt))l{x x^),{z—Zq)[{x — ®()} = 0. (All cones with vertex 

(z'o.VoiZo)-) 

20. xyz = 2i-/0, where fc is the volume. 

21. {ax-\-by—z)[a+b—ab) = obh, whore h is the given sum. 
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A^S^VEBS TO EXA3IPLES 


CBAPTBS XII 

Art 70, page 118 Ejc 2 * = 0 Ex, 3 =» y^—a^q 

Ei 4. fq^tp =*0 Ex. 5. » « *V+/W+^)- 
Ex. 6 * =/(s)co6ai:+^)saitt* 

Art 71. page 119 Ex t s *= B«-^ttn{fcr;«V 
Ex. 3 y =» BsiBjtfsiaflp* Ex 4. z = 

Art. 72. pag© 121 Ex. 1 * * E,ly)+E^+x) 

Ex. 2 z = xE,{y-3x)+E/y-ax) 

Ex 3 s = E,(y)+E^-x)+E,(y+3r) 

Ex 4 z « ru+r)+4(q+u)+pf(y+u)+4(y-ti)-i4<y~u). 

Art 73 page 122. Ex 2. x = «*»/l2-**/lS4-xr,<y-2*)+E^-Sz) 
Ex. 3 **(sm»)3e+<xtc»v)«+r,(y+2*)4-r^-3*) 

Ex 4 »*l«»*iBy)'2-i-E,ty+x)4-E,l»-x) 

rage 124 Ex. 7 * = (Sx+y)*/S10+E,(yHE^y+3x)+E^-4x) 
Ex. 8 *«. {2(3x-4yWog(lr-4y)-3(3x-4y)>3/21164- 

+*Ei(»+5rJ+Ate+W 

ExBtnplea xt end. Page 133 

1 2 * - •^♦•'wheixw » »(I+JKI-/) 

3 « - ««»+" of * - f**-*-* » 

4 *-*>/«-x*y/S+E,<y4.x)+iV(y-*) 

6 z . 28z*/2Ifl-xV/»+E,(ty+x)+E^3y-7a) 

« »-aVe+*‘-«+^«ty+Sf)+E,(y-2f) 

7 x-|eoafx-y)+/',(y+x)+xE^+x)+*»E^+x) 

® irk>0, 

vhere « letX.f<f iIIktO. 

le,MaV(— if Ir < 0. 
E*s/Cy 0. Z— and i+|+p» « 0 

10 x»2(<mx-(szn2f)/24-(ain3f)/3-_} 

11 *=|-t{{c«i)/i*+{08e3,)/3»+jeaa5x)/5«+-} 
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Adjoint equations, CG. 

Atmospheric pressure, 3S. 

Auxiliary equation, 47. 

Ball, W. W. K., 134. 

Bernoulli, equation of, 16. 

— James, John, Daniel, 135. 

Bessel, 00, 70, 73, 76, 133. 

Boundary conditions, 7, 

Cojori, F., 134. 

Cntenarj', 85. 

Cauchy, 109, 136. 

Cayley, 130. 
e-disoriininant, 27. 

Change of variable, 1(1, 81. 
Chamotoristics, 108. 

Charpit, 114. 

Clairaut, 22, 111, 135. 
Complementary function, 47, 119. 
Complete integral. 111. 

Compound interest low, 36. 
Curvature, 0. 84. 

Cusp locus, 28. 

Cylindrical harmonics, 133. 

d’Alembert, 135. 

Degree of a differential equation, 2. 
Dirichlot, 125. 

Discriminant, 20. 

Elastic vibrations, 124. 

Electricity, flow of, 125. 

Elomcntaiy functions, 1, 67. 
Elimination, 6, 105, 106, 117. 
Envelope, 24, 112. 

Euler, 136. 

Exact equations, 9, 64. 

First integral, 65. 

Forsyth, 113. 

Fourier, 120, 137. 

General integral. 111. 

Harmonic motion, 61, 80. 

Heat, flow of, 124. 

Homogeneous equations, 11, 101. 

— linear equations, ordinary, 62 : 
partial, 119. 

Hypergcomctric equation, 74. 
^J’poeycloid of 4 cusps, 32. 


Ima^ary roots of auxih'ary equation, 
48. 

Ince, E. E., 25, 134. 

Indicial equation, 70. 

Infinite series, 68. 

Intcgrability condition, 99. 
Integrating factors, 12, 66. 

Lagrange, 80, 107, 114, 136. 

Laplace, 125, 129, 130, 137. 

Law of cooling, 37. 

— of gravitation, 87. 

— of growth, 34. 

— of mass action, 37. 

— of motion, 87. 

Legendre, 72, 73, 75, 132. 

Leibnitz, 134. 

Line element, 3, 95. 

Linear dependence, 46. 

— equations, 14, 46, 96. 

Mathematical physics, 124. 

Newton, 37, 87, 134. 

Newtonian potential, 126. 

Node locus, 28. 

Normal form of second-order linear 
equation, 79. 

Normal length, 38, 40. 

Oblique trajectories, 40. 

Operators, D, /(D), 49j 1//(D), 61; 
0, 63. 

Order of a differential equation, 2, 
74. 

Orthogonal trajectories, 40, 41, 100. 

Parametric equations, 21. 

Particular integral, 47, 61, 52, 63, 121. 
p-discrirainnut, 25. 

Peirce, B. O., 74. 

Plano element, 112. 

Poisson, 126. 

Primitive, 0. 

Reduced equation, 46, 78. 

Resonance, 89. 

Riecati, 74, 75, 135. 

Separation of variables, 10, 111. 

Series, solutions in, 67. 

Simultaneous equations, oi, 90. 
Singular solutions, 24, 111. 

Special integrals 113. 
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Spl>enca]luirrao(uca, 130 137 
SuboomaU 38, 40 
Sube diary eqoatuHU 108 
Sabtugent 3S S3 
GiUp«iision bridge, 87 

Tm locus, 28 
Tangent lengdi, 38 39 
Teylot 137 
Trftctnx,39 


T^ajectonee, 40, 100 
TnaMcddentel fiuwtuioa, 48. 

Vaneluoofparamotci^lS 80 
Yeloaty of escape, 38 

Wavematioo. 118 
Webster A O 124 

Zonal hanaoolca, 132 



